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Zusammenfassung 



Das Thema dieser Dissertation ist die Starke Wechselwirkung zwischen Gluonen und Quarks 
mit Schwerpunkt auf den nicht storungstheoretischen Aspekten des Gluonensektors. Es wer- 
den hier Kontinuumsmethoden verwendet, um insbesondere das Phanomen des Farbeinschlus- 
ses zu untersuchen. Der Farbeinschluss, welcher die Detektion der elementaren Quarks und 
Gluonen als freie Teilchen verhindert, verlangt ein Verstandnis der langreichweitigen Wech- 
selwirkungen. In der Storungstheorie konnen nur kurzreichweitige Korrelationen verlasslich 
beschrieben werden. Ein nicht storungstheoretischer Zugang ist durch das Dyson-Schwinger- 
Integralgleichungssystem gegeben, das alle Greenfunktionen miteinander verkniipft. Eine L6- 
sung fiir den Gluonpropagator wird im asymptotisch infraroten und ultravioletten Grenzwert 
erzielt. 

In Kapiteld] werden redundante Freiheitsgrade der Yang-Mills-Eichtheorie durch Fixierung der 
Weyl- und Coulombeichung vor der Quantisierung entfernt. Die Quantisierung mit Zwangs- 
bedingungen unter Verwendung von Dirac-Klammern wird explizit ausgefiihrt. Als Resultat 
erhalt man den Yang-Mills-Hamiltonoperator. 

Der asymptotische Infrarotlimes der Korrelationsfunktionen in Coulombeichung wird in Kapi- 
tel[2]im Rahmen des Gribov-Zwanziger-Szenarios fiir den Farbeinschluss analytisch untersucht. 
Das Coulombpotential zwischen schweren Quarks als Teil des Yang-Mills-Hamiltonoperators 
wird in diesem Limes berechnet. Die Ubertragung der Losungen in Coulombeichung auf den 
Infrarotlimes der Landaueichung wird diskutiert. 

Der hergeleitete Hamiltonoperator ermoglicht die Bestimmung des Vakuumwellenfunktionals 
mithilfe des Variationsprinzips in KapitelO Numerische Losungen der Propagatoren in diesem 
Vakuumzustand werden besprochen, und es wird aufgezeigt, dass der vorhergesagte Infrarotli- 
mes tatsachlich realisiert ist. Die Diskussion wird auf Vertexfunktionen erweitert. Des Weiteren 
wird der Einfluss der Naherungsmethoden auf die Losungen untersucht. 

KapitellHist vornehmlich dem Ultraviolettverhalten der Propagatoren gewidmet. Die Behand- 
lung erfolgt sowohl in der Coulomb- als auch in der Landaueichung. Ein nicht storungstheore- 
tischer laufender Kopplungsparameter wird definiert und berechnet. Der ultraviolette Teil der 
Variationslosungen aus Kapitel [3] wird mit den Forderungen der Storungstheorie verglichen. 

In KapitelOwird die Riickkopplung des Gluonensektors auf die Anwesenheit aufierer Ladungen 
(schwerer Quarks) behandelt. Zu diesem Zweck werden koharente Gluonenanregungen vorge- 
schlagen, welche der Anwesenheit aufierer Ladungen Rechnung tragen. Weitere Alternativen 
zur Behandlung dieser Problemstellung werden besprochen. 



Abstract 



The subject of this thesis is the theory of strong interactions of quarks and gluons, with par- 
ticular emphasis on nonperturbative aspects of the gluon sector. Continuum methods are used 
to investigate in particular the confinement phenomenon. Confinement — which states that the 
elementary quarks and gluons cannot be detected as free particles — requires an understanding 
of large-scale correlations. In perturbation theory, only short-range correlations can be reli- 
ably described. A nonperturbative approach is given by the set of integral Dyson-Schwinger 
equations involving all Green functions of the theory. A solution for the gluon propagator is 
obtained in the infrared and ultraviolet asymptotic limits. 

In chapter [H redundant degrees of freedom of the Yang-Mills gauge theory are removed by 
fixing the Weyl and Coulomb gauge prior to quantization. The constrained quantization in 
the Dirac bracket formalism is then performed explicitly to produce the quantized Yang-Mills 
Hamiltonian. 

The asymptotic infrared limits of Coulomb gauge correlation functions are studied analytically 
in chapter [2] in the framework of the Gribov-Zwanziger confinement scenario. The Coulomb 
potential between heavy quarks as part of the Yang-Mills Hamiltonian is calculated in this 
limit. A connection between the infrared limits of Coulomb and Landau gauge is established. 

The Hamiltonian derived paves the way in chapter [3] for finding the Coulomb gauge vacuum 
wave functional by means of the variational principle. Numerical solutions for the propagators 
in this vacuum state are discussed and seen to reproduce the anticipated infrared limit. The 
discussion is extended to the vertex functions. The effect of the approximations on the results 
is examined. 

Chapter H] is mainly devoted to the ultraviolet behavior of the propagators. The discussion 
is issued in both Coulomb and Landau gauge. A nonperturbative running coupling is defined 
and calculated. The ultraviolet tails of the variational solutions from chapter [3] are compared 
to the behavior demanded by perturbation theory. 

In chapter EJ the back reaction of the gluon sector on the presence of external charges (heavy 
quarks) is explored. To this end, coherent excitations of gluonic modes are suggested to account 
for the presence of quarks. Further alternatives for the discussion of this issue are put forward. 
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1. 



Yang-Mills theory as a constrained dynamical system 



1 Yang-Mills theory as a constrained dynamical 
system 



In the development of quantum electrodynamics (QED) [T], local gauge invariance of the matter 
fields served as a guiding principle and it successfully led to a highly accurate description of 
electromagnetic forces pj. The fundamental requirement of locality necessitates the existence 
of bosonic gauge fields which mediate the forces between the fermions. An impediment to 
the canonical quantization of such a dynamical system is the existence of redundant degrees 
of freedom. This obstacle can be overcome for the electromagnetic theory, thanks to the 
linear structure of its equations of motion. Quantum chromo dynamics (QCD), on the other 
hand, bears non-linearities that make the treatment far more complicated. Since QCD is the 
non-abelian generalization of QED, the general structure of a gauge theory is inherited and 
along with it its difficulties. It will be shown in this chapter how the redundant degrees of 
freedom can be related to constraints in the classical phase space, and how quantization of 
such a dynamical system can be pursued. After a brief definition of Yang-Mills (YM) theories, 
amongst which QCD is the SU (3) version, we first examine constrained systems in the context 
of classical mechanics to pave a way for the terminology. Second, constrained quantization is 
exhibited for electrodynamics as an illuminating example. Finally, we demonstrate how to 
systematically follow the steps of constrained quantization in the difficult case of Yang-Mills 
theory in Weyl and Coulomb gauge. The explicit form of the quantized Hamiltonian operator 
will thus be derived in a quantization approach that differs from the canonical method. 

1.1 Definition of Yang— Mills theory 

After a first sporadic attempt by Klein as early as 1939 to generalize the electromagnetic theory 
imposing non-abelian group structure [3], it took another 15 years until Yang and Mills were 
able to provide the full account of SU{2) non-abelian gauge theory [4]. Its generalization to 
SU (Nc) (with Nc = 3 realized in QCD) is straightforward, a clear introduction to the currently 
accepted definition of YM theory can be found, e.g., in Ref. [5]. There are some good text 
books on this topic, e.g. Ref. [6], yet the main elements of YM theory are outlined here to set 
up a notation for this thesis. 

The basic idea, following electromagnetism as a guiding example, is to cancel the change of 
the Dirac Lagrangian 



under a local symmetry transformation of its fermion fields ipa{x) in the intrinsic spaci 



Cd = il; {i^d^, -m)il) 



(1.1) 



Tpaix) Ipa (x) = Uabix)lpbix) 



(1.2) 



^The roman color indices are used as subscript or superscript interchangeably. 
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1. Yang-Mills theory as a constrained dynamical system 



by the introduction of a gauge field ^JJ(x). To achieve this, the derivative 3^ = g|jr in the 
Lagrangian (jl.lj) is replaced by 

Df, := + gAf, (1.3) 

which transforms homogeneously, UDyU\ whereas the matrix- valued gauge fields 

transform inhomogeneously, according to 

A^{x) = A1{x)T^ ^ A^ (x) = U{x)A^{x)U\x) + -U{x)d^uHx) ■ (1.4) 

The subsequent interaction of gauge and matter fields is controlled by the coupling constant g. 
Originally, this procedure coined the term "minimal substitution". Nowadays, with a deeper 
understanding of the geometrical aspects which allow close analogies to general relativity 
we call the "covariant derivative", defined in Eq. p.3p as acting on a field in the fundamental 
representation, such as the quark field iIj°'{x). 

The unitary group elements U{x) G SU{Nc), 

U{x) = exp (-a"(x)r'') , a e M , (1.5) 

are continuously generated by N'^ — 1 antihermitian matrices Ta = —T^ that obey the Lie 
algebra 

[r",r^] = r''^r^ (1.6) 

with the structure constants For the quarks to transform as 3- vectors in the fundamental 

representation of SU{3), one uses the well-known 3x3 Gell-Mann matrices for iT". The 
gauge fields A^, on the other hand, live in a (A'^^ — l)-dimensional vector space spanned by the 
generators with the additional operation of commutation (|1.6p . By an infinitesimal expan- 
sion of (11.41) in a using (|1.6|) . one can see that must transform in the adjoint representation 
of SU{N,), 

{A^y^ = A''^ + {f'rA'l.a' + ^d^a'' + ©(a^) , (1.7) 

where we defined 

^fbyc ._ jabc _ (J 

Any linear combination of the matrices with field operators shall be denoted by a caret, e.g. 
A^ = A'^T'^, as opposed to the matrices A^ introduced in Eq. p.4p . The covariant derivative 
D'j^ acting on fields in the adjoint representation is defined by 

Df = 6-% + gAf. (1.9) 

Knowing the transformation properties of the gauge field p.4p . one may construct from the 
gauge covariant field strength tensor 

F^,, := ^[Df„D,] = d^A, - d^A^^g\A^,A,\ , F^, ^ UF^.U^ (1.10) 
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1.1. Definition of Yang-Mills theory 



a gauge invariant contribution to the Lagrangian density that is purely gluonic. It is the first 
term in the Yang-Mills Lagrangian 

Lym = ltTF^,{x)F^^''{x)+gr^{x)A^Ax) , (1.11) 

whereas the second term represents the couphng to the quark Nother current of the color 
symmetry, 

j^ix) =^{xh^iT^^{x) , (1.12) 

here treated as a classical external field. This thesis shall mainly be concerned with the 
dynamics of the Yang-Mills Lagrangian (jl.lip and disregard the influence of dynamical quarks. 
It is hypothesized that unquenching has a small effect on the gauge field sectoral The color 
trace in Eq. (11.111) with the convention tr(T"r^) = — ^5"* guarantees that the gluonic energy 
is positive definite, cf. (jl.lSp . 

With the action principle, the Lagrangian formalism provides a possibility to derive equations 
of motion where Lorentz invariance is manifest. Applying a functional derivative w.r.t. the 
gauge field to the Yang-Mills action 

S = J d^xCvMix) (1.13) 

one findj§ 

DfF^ = -gfa (1-14) 
and, using the Bianchi identity [D^j, F^p] + cyd.perm. = , 



^f^yp'^fjabpa^ = . (1.15) 



In order to derive a Hamiltonian, one introduces conjugate momenta 
and readily notes that due to the antisymmetry of F^^ 

ng(x) = o. (1.17) 



^The nonperturbative confinement phenomenon is known to exist in the absence of flavors, Nf = 0, from 
lattice calculations [8j. In the perturbative regime, a small number of quark flavors {Nf < 16) will not remove 
the property of asymptotic freedom. Recent continuum studies [9J as well as lattice calculations [lOJ indicate 
that the effect of unquenching on hadronic observables is small. 

^Note that F^'' transforms in the adjoint representation, therefore it is natural that the covariant derivative 
Du is found here. 
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1. Yang-Mills theory as a constrained dynamical system 



Eq. ()1.17|1 is understood as a constraint on the phase space variables and will be shown to 
complicate the quantization process. The Legendre transform gives the Hamiltonian H, 

H = j d\ {lild^A^, - Cym) 

= I d^x (-^n^n^ + ^i^^F^^' - Am'^t - 9A^oPt.i - 9Atj^) 

= ]^jd\{lll + Bl)-jd\Aig- + gjd\A--i\ (1.18) 

with the abbreviations for the external color charge density Pg^t = Jo > color magnetic field 
Bt = -hijkFt, and 

g^{x) = bfixHix)+gpt^,{x). (1.19) 

In the Lagrangian equations of motion we recognize that the Gauss law can be written as 
Q"" = 0, see Eq. (|1.14ll for u = However, in the Hamiltonian approach, the Gauss law 
cannot be found as an equation of motion but must be implemented as a constraint, as will 
be discussed below. 

The quantization of a gauge theory is a delicate task, already in the abelian case. Since the 
gauge transformation indicates that the field variables A'^{x) contain unphysical degrees of 
freedom, it is not clear a priori whether the canonical quantization leads to the correct results, 
i.e. whether it contains the classical theory in the limit ^ — > 0. It is discussed in some detail 
in the following sections how a gauge theory can be understood as a constrained system and 
how to go about with its quantization. 

Although in the Lagrangian formalism Lorentz invariance can be made manifest, it is some- 
times favorable to use the Hamiltonian formalism instead. The path integral method provides 
an elegant way of quantizing a theory covariantly in the Lagrangian formalism. However, 
it is rather unsuitable to calculate, e.g., the Balmer formula for hydrogen, one of the first 
exercises in quantum mechanics [Hj. The Hamiltonian approach to the quantum theory of 
electromagnetism and YM theory shall be the main focus of this work. Once the conjugate 
momenta are defined and one passes over to the Hamiltonian, Lorentz covariance is lost since 
a reference frame has been chosen. Nevertheless, Lorentz invariance is maintained. This can 
be understood, e.g., by explicit calculation in Coulomb gauge for electromagnetism [12] or by 
criteria for commutation relations in Coulomb gauge quantum Yang-Mills theory [13] . 

1.2 Constrained dynamics 

In this section we shall discuss the importance of constraints of a dynamical system for its time 
evolution and for its quantization procedure. The basic idea was developed by Dirac [HI [15] 

■^In electrodynamics, one defines the electric field as E'' = -IT* and thus gains the familiar form of the 
Maxwell equation V • E = ep. 
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1.2. Constrained dynamics 



and very pedagogically presented in his Yeshiva lectures |T6j. Dirac and many of his successors 
tried to maintain a high level of generality in favor of the applicability to any given system, 
say electrodynamics or even general relativity [TT]. However, it turned out that exceptions to 
the quite general claims can be found. Therefore, only the concepts relevant to Yang-Mills 
theory will be introduced here. For additional information, see e.g. [T8l IT91 120] . 

Consider a classical mechanical system described by a Lagrangian L{qi,qj). The Legendre 
transformation to a Hamiltonian H{qi,pj) with conjugate momenta pk can be achieved if we 
can solve the equation 

dL 

«:=^ (1.20) 

for the velocities qk, which is locally guaranteed by a regular Hessian 

det (^^^) / . (1.21) 



dqidqj 

The phase space T is then defined as a set of independent elements {qi,Pj}- In the case of 
a singular Lagrangian where Eq. p.21|) does not hold, some of the momenta pk are depen- 
dent. Moreover, additional dependences among the qk and pk may arise from the Hamiltonian 
equations of motion. Generally, one may formulate all dependences as a number of constraints 

^m{q,p)-0- (1-22) 

by which the dimension of phase space is reduced. Below, the symbol is explained. The 
constraints (ll.22p then hold in a subset C F of the original phase space. The Hamiltonian 
H can be supplemented by the constraints (ll.22p with some arbitrary Lagrange multipliers 

Ht = H + VmVm (1-23) 

We refer to Ht as the "total Hamiltonian". The common variational calculus with constraints 
gives rise to Hamiltonian equations of motion in F/j, 

qk = ^ + vm^ ^ {q,HT} (1.24a) 

Opk OPk 

Pk = -^-Vm^ ^ {P,Ht} (1.24b) 
dqk dqk 

Above, we have made use of the Poisson bracket for any functions f{q,p) and g{q,p), 

df dg df dg 
dqk opk opk dqk 

Note that there is a 'ss' symbol in Eq. (|1.24l) . a notation due to Dirac which emphasizes that 
one is to first calculate the Poisson brackets and then set the constraints (|1.22|) to zero. These 
kinds of relations are termed "weak equations" and are specific to the use of Poisson brackets. 
It may be stressed here that the canonical quantization with such weak equations leads to the 
circumstance that some quantum constraints cannot be imposed as operator identities but only 
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1. Yang-Mills theory as a constrained dynamical system 



as projections on states. Below, we shall introduce generalized Poisson brackets that make the 
usage of weak equations redundant. 

The time evolution of any given phase space function g{q,p) that does not explicitly depend 
on time is given by the Poisson bracket in a concise way: 

g^{g,HT}. (1.26) 

The total Hamiltonian Ht obviously serves as the generator of time evolution. However, the 
multiplier functions Vm contained in it leave an arbitrariness. Given some initial conditions, 
say g{0) = go, g{t) at some finite time t is not unique since we can choose the Vm to be one 
value or another. After an infinitesimally small time 6t, the function g will yield 

g{dt) ^go + {go,H}6t + v^igo, <fm}dt . (1.27) 

Imagine we take two different sets of multiplier functions v!^ and v'^. Comparing these two 
arbitrary choices to calculate g'{St) and g"{6t), using (|1.27p . leads to 

g"{6t) « g'{6t) + em{go, ^m} (1.28) 

with em = {v'm — v'm)5t. Evidently, the two physically equivalent values g'{5t) and g"{6t) 
are related by a term generated by the constraints ipm- This suggests that the physical state 
remains unchanged under transformations of the kind 

9(.q,p) 9^(.q,p) = 9(.q,p) + ^m{g, v'm} (i-29) 

However, we have to be careful with the interpretation, as shown below. 

The constraints (|1.22p need to be checked for consistency, i.e. they should hold at all times. 
In view of Eq. (I1.26P this requirement gives rise to the conditions 

'^k = {^k,H} + {Lpk,ipm]Vm'^^ (1-30) 

which can actually fix the multiplier functions Vm so long as the matrix 

Ckm ■= Wk, Vm] (1-31) 

is invertible. We now distinguish first-class and second-class constraints |16j. By definition, a 
first-class constraint has vanishing Poisson brackets with all other constraints. With first-class 
constraints present, the matrix {Ckm) in Eq. p.31|) becomes singular, det {Ckm) ~ 0. In the 
absence of first-class constraints, all constraints are termed second-class, and one can show |21j 
that 

det(Cfc„) ^ . (1.32) 

If Eq. (I1.32P holds and only second-class constraints are present, we can fix the multiplier 
functions Vm using the consistency conditions (ll.30p by 

Vk = -C^mym,H} (1.33) 
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1.2. Constrained dynamics 



and derive that g{t) with second-class constraints satisfies the equation of motion 

g = {g,H} - {g,ipk}Cj:^W.^,H} ^ {g,H}D • (1.34) 
Here, the Dirac bracket was introduced, 

{f,9}D ■■={f,9}-{f,Vk}C^^,W^,9} ■ (1.35) 

It is clear from Eq. (|1.34l) that the time evolution of the initial state go is unique. The Dirac 
bracket provides a formalism that yields equations of motion with (second-class) constraints. 
Furthermore, the constraints (|1.22p which are regarded as weak equations w.r.t. Poisson brack- 
ets, can be used as strong (ordinary) equations if one uses the Dirac bracket instead, since for 
any function / 

Wk, f}D = WkJ} - Wk, (pj}C~^{ipm, /} = . (1.36) 

These ideas led Dirac to propose a quantization technique for constrained dynamical systems, 
known as constrained quantization^, making use of the Dirac brackets in the fundamental 
commutation relations of operators Qi and j5j, 

ih{qi,pj}D [qi,Pj] ■ (1.37) 

While the Dirac brackets has most of the features of the Poisson bracket — e.g. antisymmetry, 
linearity, Jacobi identity, product rule — the fundamental Dirac brackets are generally not as 
simple as the Poisson brackets, {qi,Pj}D / ^ij- In linear theories, such as QED (discussed in 
the next section) this does not pose a problem to the quantization prescription p.37|) since 
the fundamental Dirac brackets do not involve field variables. However, in the non-linear YM 
theory, constrained quantization is more difficult. This will be dealt with in section [L4l 

Let us come back to the claim that constraints generate transformations among physically 
equivalent states, see Eq. (|1.29|) . For second-class constraints, this can certainly not be correct, 
since all multiplier functions are fixed, as we have seen in Eq. (|1.33p . To understand this caveat, 
note that the transformations (I1.29P generated by second-class constraints kick the system out 
of the reduced phase space F/j since for g = (pk 

V>k = ^'k + Wk,^m}em 7^ . (1.38) 

That is, the transformations (I1.29P generated by second-class constraints are not symmetries 
of the theory. 

On the other hand, if we consider first-class constraints only, then the Poisson brackets all 
vanish, {ifi, ifj} = 0, and the transformations p.29|) are indeed such that the system remains in 
the reduced phase space Tr where the constraints are satisfied. They give rise to an equivalence 
clasfl 

[(9,^)] = {(9^P') er^l ('Z,!?) ~ ('7^/) , e€R} (1.39) 
'it is also associated with the lengthy expression "reduced phase space quantization". 

®The equivalence relation ~ is here represented by the transformations infinitesimally defined in Eq. Ill.29p . 
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1. Yang-Mills theory as a constrained dynamical system 




Figure 1.1: Left: Ambiguity in time evolution of an initial point {qo,Po) in 
the reduced phase space Tr. Different choices (v, v' , v") of the arbitrary 
multiplier functions lead to different trajectories in Tr. Right: At a fixed 
time, the gauge orbit (solid line) represents physically equivalent states. 
The gauge fixing condition x is admissible, whereas x' ^Siils to be unique 
and x" is not attainable by a gauge transformation. 



of phase space variables that all correspond to the same physical state. Any set of {qi,Pj) 
uniquely determines the state but the reverse is not true, see Fig. II. 1[ These considerations 
infer that the first class constraints are the generators of symmetries of the theory. 

It should be noted here that the above statement, known as "Dirac's conjecture" is not a 
rigorous theorem. Counter examples can be constructed [22], and it takes further classifications 
of constraints ("primary" and "secondary") to identify the exact set of generators of symmetry 
transformations in general. However, all physical applications agree with Dirac's conjecture 
and the issue is somewhat academic. In particular, the gauge transformations of QED and 
YM theory are generated by first-class constraints (see below). In this context, the equivalence 
class (|1.39ll is also referred to as the gauge orbit. 

The quantization in the presence of first-class constraints is performed either in the path- 
integral formalism [23], or by a projection on the physical Hilbert subspace [5l[2l]. Alterna- 
tively, one can fix the gauge on the classical level and then quantize with the Dirac prescription 
(|1.37|1 . When fixing the gauge, one effectively picks out a single representative from the equiv- 
alence class (|1.39ll generated by the first-class constraints. In practice, this is achieved by 
imposing supplementary constraints ("gauge conditions") 

Xn ~ (1.40) 

that turn the first-class constraints ipm into second-class ones and hence obey 

det{{xn,^m})^0. (1.41) 

The constraints ifm that formerly were first-class thus cease to cause an ambiguous time 
evolution of the system. With all constraints being second-class, all multiplier functions are 
fixed and one can use the Dirac bracket to proceed with the quantization. 
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1.3. Quantum electrodynamics 



In order to arrive at a "physical gauge" where all unphysical degrees of freedom are eliminated, 
the number of gauge conditions must equal the number of first-class constraints, hence the 
matrix that appears in Eq. (|1.4ip is quadratic. In principle, the gauge conditions (|1.40p are 
quite arbitrary, except for the requirements of uniqueness and attainability. Uniqueness can 
be established locally as follows: Consider in Fig. 11.11 (right panel) the gauge orbit at a fixed 
time tf and assume without loss of generality that the point {q^ ,p^) satisfies Xnil^ :P^) = 0. 
The transformation 

Xn^Xn=Xn + (im{Xn,H^m] (1-42) 

within the gauge orbit is not to yield any further solutions Xn ~ unless it is the identity 
transformation, e = 0. Otherwise, the gauge fixing condition specifies (at least) two elements 
of the gauge orbit, see x' in Fig. 11.11 With the condition p.41|) on the Xn, one can see from Eq. 
(I1.42P immediately that e = is the only solution to Xn ~ ^nd therefore the point {q^ ,p^) is 
unique0 Moreover, a sensible gauge condition has to be attainable by a gauge transformation. 
Otherwise, we would have the situation in which the gauge fixing condition does not intersect 
with the gauge orbit at all, see x" in Fig. 11.11 

In summary, the general procedure of quantizing classical mechanics in the presence of con- 
straints has been discussed. It has not been shown how to find the entire set (ll.22p of con- 
straints starting from the Lagrangian. This can be achieved with the so-called Dirac-Bergmann 
algorithm [lH[15l[25]. It incorporates the stationarity of constraints that are found in the equa- 
tion of motion or in the mere definition of conjugate momenta. 



1.3 Quantum electrodynamics 

An illuminating example of the formalism exhibited above is the quantization of electrody- 
namics in Weyl and Coulomb gauge. There are no difficulties in promoting the ideas from 
point mechanics to a field theory, the discrete indices of the generalized coordinates are merely 
replaced by the continuous spacetime dependence. 

The electromagnetic theory is the abelian version of the general set of YM theories described 
in section 11.11 recovered by setting f"'^'^ = in the Lie algebra (II. 6|) . In the absence of external 
charges, one finds from (II. lip the well-known Lagrangian density 

C = -^F^uF^-" , F^, = d^A, - d,A^ . (1.43) 

As quoted in Eq. p.l7|l . the conjugate momentum of the Aq field vanishes and we can write 
down a first constraint fi aj§ 

Mx) ■■= no(x) = « . (1.44) 



^The global issues of uniqueness will be discussed below in the context of YM theory. 
*In field theory, an infinite number of constraints is specified, one at each spacetime point. 
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1. Yang-Mills theory as a constrained dynamical system 



The Hamiltonian function depending on the fields ^0(2;), Ak{x) and the conjugate momenta 
n^(a;) reads 

H = d^x (n2(x) + b2(x)) - j d^x AQ{x)dkYl^{x) . (1.45) 

From now on, we write Hk{x) for the components of the contravariant tensor n'^(x) = F^^{x). 
This change in notation will become useful in the nonabelian case where the field operators 
also have a color index. The Hamiltonian (|1.45l) determines the time evolution via the Poisson 
brackets. The fundamental Poisson brackets are given by 

{A^{x,t),Ii,{y,t)} = 5^,5\x,y) . (1.46) 

Since the constraint (ll.44p is stationary, 

992(x) :={no(x),ifT} = {no(x),if + J d^yvi{y)Uo{y)} = dkUk{x)^0 (1.47) 

has to hold weakly where vi (x) is an arbitrary multiplier function. We recognize Eq. (|1.47p as 
the Gauss law. Following the Dirac-Bergmann algorithm, we can check if further constraints 
arise from requiring that (p2 is constant in time. One finds 

!^M^),H + J d^y{vi{y)Uoiy)+v2{y)dkUkiy))^ = (1.48) 

and thus the complete matrix of constraints Cij, see Eq. (|1.31ll . is given by (pi and 932 which 
are obviously first-class, {y?i,(/?2} = 0. 

As we have shown in the preceding section, the first-class constraints serve as generators of 
gauge transformations. Here, it can be explicitly verified that the transformations A^ 
A^ + 5A^ with 

5^^Aq{x) = I (i^x'{Ao(x),v3i(x')}ai(x') = ai(x) , 5^^Ak{x) = (1.49a) 



5^^Ak{x) = j d^x' {Ak{x),ip2{x')]a2{x') = -dka2{x) , S^^Ao{x) = (1.49b) 

are symmetries of the Lagrangian (I1.43P . The general form (11.4(1 of gauge transformations in 
the Lagrangian formalism, with generator T = zl of the group U{1) and gauge coupling g = e, 
reads 

^^(x) ^ ^^(x) = A^{x) + ^d^a{x) . (1.50) 

If we wish to identity an a(x) in Eq. (jl.SOp that corresponds to the transformations p.49|) . 
it is recognized that Eq. p.49ap corresponds to spatially independent gauge transformations 
a{t) with ai{t) = doa{t)/e whereas Eq. (I1.49bp corresponds to time-independent gauge trans- 
formations with a{x) = —a2{x)/e. 

In the quantization procedure, the problem arises that due to the gauge invariance (|1.49|) there 
are more variables than physical degrees of freedom. The constraints (|1.44p and (I1.47P then 
should be taken into account in the Dirac formalism. Before we do so, it is briefly sketched 
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1.3. Quantum electrodynamics 



here how the canonical quantization is usually pursued with Poisson brackets. With the choice 
of Weyl gauge, Aq{x) = 0, the first constraint, no(x) ^ 0, does not cause any difficulties for 
the quantization of the time-components of the fields. The canonical Poisson brackets are 
used to promote the spatial classical fields to operator-valued fields that obey the equal-time 
commutation relations 

[Ai{x),Uj{y)]=i6ij5^{x,y) . (1.51) 

Since the Gauss law (|1.47p is an equation that holds only weakly on the classical level, it is 
not a surprise that it leads to a contradiction with the commutation relations, 

= [Ak{x),dJUj{y)] = dy[Ak{x),U,{y)] = idl5\x,y) + . (1.52) 

The requirement that the Gauss law holds as an operator identity is therefore abandoned and 
one usually restricts the Hilbert space to the kernel of the Gauss law operator Q, 

9^nfc(y) 1^-) = 0, (1.53) 



Let us mention here that the quantum theory is described in the Schrodinger picture where 
the field operators are independent of time. 

Constrained quantization has the generic feature that the constraints can be imposed as strong 
equations once the Poisson brackets are replaced by the Dirac brackets. In the so quantized 
theory, the constraints hold as operator identities and can be used to eliminate unphysical 
degrees of freedom. Before the Dirac brackets can be introduced, all first-class constraints 
have to be degraded to second class by imposing gauge fixing conditions. We now choose Weyl 
and Coulomb gauge to render the constraints (|1.44p and p.47|) second class. 



Xi(x) =^o(a;) «0 
X2{x) = dkAk{x) « 



(1.54a) 
(1.54b) 



The Weyl or "temporal" gauge condition p.54ap fixes the gauge transformations (ll.49a|l gen- 
erated by the first-class constraint Ho ~ 0, whereas the Coulomb gauge condition (ll.54bp fixes 
the gauge transformations (|1.49bj) generated by the first-class constraint dkllk{x) which 
is the Gauss law. Thus, the gauge is fixed completeljlf] on the classical level and we will refer 
to the joint conditions ()1.54ap and (|1.54bj) as the temporal Coulomb gauge. 

All (second-class) constraints are now collected in a vector (p = {(pi, ^2-,Xi-,X2), and we redefine 
the matrix p.31|) by Cij{x, y) = {(^j(a;), (j)j{y)} which can be evaluated for a fixed time t using 
the fundamental Poisson brackets (I1.46P to yield 



{Cij){x,y) 



( 



b^ix.y) 
\ 










-d'^5^{x,y) 



ix,y) \ 

d'^6''^{x,y) 



/ 



(1.55) 



Up to global gauge transformations. 
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1. Yang-Mills theory as a constrained dynamical system 



Since Cij is regular by construction, we can also write down its inverse, 



V 
























1 

4TT\x—y\ 






(1.56) 



ATT\x~y\ 

The fundamental Dirac brackets can be calculated straightforwardly from the definition p.35|) . 
{A{x),Uj{y)}D = {A{x),Uj{y)} 

- 1 d^x'dy{A{x),ct)^{x')}C-,i{x',y'){My'),^jiy)} 



5,,5\x,y)- [ d'x'd'y'df6{x,x')—^ ndf5iy',y) 

J 4:Tt\x' — y'\ ■' 



5i,5\x,y)-dfd^- 



1 



tij{x)6^{x,y) =: tij{x,y) (1.57) 



^ Att\x — y\ 

where the last line defines the transverse projector tij{x). The definition of the matrix- valued 
distribution tij{x,y) will be useful as well. 

Quantization is achieved by the prescriptior0 



{Ai{x),Uj{y)}D =tij{x,y) 



[Ai{x),Uj{y)] = itij{x,y) 



(1.58) 



which provides information on how the momentum operator 11^(03) acts on the wave functional 
(^1 = "^[A] in the coordinate representation. Yet, not all components of 11^(3;) are deter- 
mined by Eq. p.58|l since the transverse projector has zero modes, dlj.tij{x,y) = dltij{x,y) = 
0. To be specific, the longitudinal components of the field operators are left unconstrained 
by Eq. (I1.58p . We split the field operators into longitudinal and transversal components, 
Ai = Aj- + a| with A-^{x) = tij{x)Aj{x) (and for IT accordingly), and find from Eq. (|1.58p 
that 



[Ai{x),Iif{y)]=itij{x,y) 



(1.59) 



The additional information needed for Aj^ and n|J, comes from the constraints themselves, 
Eqs. (fOTll and (ll.54bh . that now hold as operator identities. The longitudinal projections by 



• ■- 



5ij — tij of both field operators are identically zero. 



Al{x)=eijix)Aj{x) = 
uj{x)=£ij{x)Uj{x) = 

Here, no contradiction to the commutation relations (I1.58P arises, 

= [Ai,{x),dJUj{y)] = dy[Ak{x),Uj{y)] = idjtkj{x,y) = , 



(1.60) 
(1.61) 



(1.62) 



^"The field operators are denoted by the same symbols as their classical counterparts. It should be clear 
from the context which objects the symbols refer to. 
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1.4. Constrained quantization of Yang-Mills theory 



unlike in the quantization procedure with Poisson brackets, cf. Eq. ()1.52|) . The transverse 
components of the momentum operator is found in view of Eq. (I1.59P to be 




(1.63) 



Quantization is thus complete and one may find from (| 1.451) the gauge-fixed quantum Hamil- 
tonian 

H[A^,U^] = ljd^x (U^^ + B'^[A^]^ . (1.64) 

After all, we have arrived at a Hamiltonian expressed in terms of gauge-invariant variables. 
Note that A-^ is left unchanged by gauge transformations (ll.SOp and Hk = F^o transforms 
homogeneously in view of (ll.lOp and hence transforms trivially in the abelian theory. 



1.4 Constrained quantization of Yang-Mills theory 

The transfer of the techniques used in the previous section to non-abelian theories is aggravated 
by non-linearities in the equations of motions. Nevertheless, constrained quantization of YM 
theory in the temporal Coulomb gauge is feasible. This section is devoted to presenting the 
steps of calculation that lead to the gauge-fixed Hamiltonian operator of YM theory with static 
external quark fields. 

Usually, the canonical quantization is performed in the Weyl gauge where the spatial compo- 
nents of the gauge field can be treated as a set of Cartesian coordinates, and after quantization 
a coordinate transformation to curvilinear coordinates in the Coulomb gauge is undertaken 
[211 ESl [23, [2E]- This corresponds to an incomplete gauge fixing before quantization and is 
known to bring about a number of difficulties [29]. For instance, if only the Weyl gauge, 
Aq = 0, is enforced, the Gauss law cannot be used as an operator equation and has to be 
imposed as a projection onto a physical Hilbert space [30l[5], see also [31]. The physical states 
are not normalizable as an artefact of the residual gauge invariance [32]. Moreover, the com- 
mutator [^",^^] may be anomalously broken which obscures the projection on states [331 [34] . 
Another difficulty in fixing only the Weyl gauge is that the perturbative gauge field propagator 
is plagued by unphysical poles [35| . 

In the present thesis, the gauge shall be completely fixecil^ already at the classical level by 
imposing the temporal Coulomb gauge conditions 

Ag(x) w (1.65a) 
lDf{x)Al{x) = dkAUx) ^ (1.65b) 

and only afterwards, the quantization will be issued with the Dirac brackets. This choice was 
also discussed in [36l[37j. Complete gauge fixing, such as the temporal Coulomb gauge (|1.65p . 

"See footnote El 
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1. Yang-Mills theory as a constrained dynamical system 



has the advantage that it can be chosen such as to remove all unphysical degrees of freedom. 
The Gauss law then holds as an operator identity. In exchange, the algebra of the quantization 
procedure is slightly more involved. 

Analogously to electrodynamics, the first constraint to momentum phase-space is given by the 
time component of the momentum operator 

The Hamiltonian (jl.lSp with j" = for simplicity, 

H = lld'x (UlUt + B^.Bt) - I d'xA^oG" (1-67) 

serves as the generator of time translations The stationarity of the constraint (|1.66p gives a 
further constraint, 

:= WUx),Ht} = Q%x) = bl\x)Il\{x)+gpl,,{x) ^ , (1.68) 
where the equal-time fundamental Poisson brackets were employed, 

[A1{x,t),Iil{y,t)]=5^J'^H\x,y) . (1.69) 
The constraint ip2 is found to be stationary using the well-known identity [31] 

{g%x,t),g\y,t)] = gr''G'{x,t)5\x,y) « . (1.70) 



The above quantity weakly vanishea^^l in view of Eq. (I1.68P and therefore the constraints ipi 
and ip2 in Eqs. p.66p and p.68p form the complete set of constraints. The latter are first-class 
and are turned into second-class by the gauge fixing condition p.65p . 

Before proceeding to the calculation of the Dirac brackets, let us introduce a convenient short- 
hand notation. The coordinate dependence is absorbed into the color index which is thus 
understood as a collective index. Vectors in Lorentz space keep their Lorentz indices to avoid 
an over-estranged notation. All fields are considered at a fixed time. Repeated indices are 
summed over, i.e. for the color indices this implies an integration over the implicit coordinate 
dependence. 

For instance, 

Ai{x)^Ai, ^^(^) = ^^ (1.71) 

In addition, derivatives are written as matrices in coordinate and color space, 

d^5'^{x,y)6''^ = df , Df{x)6^{x,y) = Of . (1.72) 

^^In the same way as in the previous section, we write n^(a::), meaning n'°"(a;). 

^^A theorem due to Dirac [14] ensures that a Poisson bracket of two first-class constraints is again first-class. 
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1.4. Constrained quantization of Yang-Mills theory 



Note that these matrices are antisymmetric under the exchange of the collective indices, 

5f = -a^ , bf = -bl^ . (1.73) 
Also note the symmetric matrices 

5'''5\x,y)^5''\ 5''%{x,y)^tf^ , <5%(^, ^ . (1.74) 

The constraints along with the gauge- fixing conditions are collected in a vector (j) with elements 



4>l = 
<^l = 

to give a regular constraint matrix C^^ 
calculation yield 

safe 



^2 
X2 



T-ra 
ga 



3 



(1.75a) 

(1.75b) 

(1.75c) 

■fc -fc - (1.75d) 

{{(f^m, 4>X\) ■ The relevant Poisson brackets for its 



(1.76a) 
(1.76b) 
(1.76c) 
(1.76d) 

[G-^^ (l-76e) 



Dl 



dj'^{Aj, n^} — uj. 



In the last line we have defined the non-local matrix which is symmetric in view of the 
antisymmetric matrices in Eq. (|1.73l) . Furthermore, due to the Coulomb gauge condition 
(ll.65bp . it obeys the relation 

[^-Ijafe 



and its inverse G is defined by 

All (weakly) non-zero matrix elements (C"*)., 
that the constraint matrix reads 



-liacQcb 



rafe 



(1.77) 



(1.78) 



and its inverse is given by 






^ab 





[G- 







-llab 



are now given by Eqs. (I1.76al) and (ll.76ep so 



-5" 








[Q-l^ab 






(1.79) 



[G- 



/ 








^ab 





\ 













Qab 






_Sab 













V 





_Qab 








/ 



(1.80) 
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1. Yang-Mills theory as a constrained dynamical system 



At this point, let us emphasize that the matrix G, which will be the classical counterpart of 
the ghost operator, can only be made explicit by an infinite series expansion in the coupling 
constant This has led some authors to refrain from the constrained quantization with 

such a constraint matrix, see Refs. |T9l EBl EH] ■ In Refs. (lOlIll], a Hamiltonian operator is 
obtained but one is left with operator ordering ambiguities. In other gauges, where Aq ^ 0, the 
expression for C""^ becomes even more complicated [20]. Our approach aims at a comparison 
to the Christ-Lee Hamiltonian [24] in the temporal Coulomb gauge. Therefore, let us press 
on with the formal expression for G defined by (|1.78ll and with the inverse constraint matrix 

arm . 

The fundamental Dirac brackets are defined by 

{AlUl}n = {AlUt}-{Air:.}{[C-Y''UnWtK} (1-81) 
and, using Eq. p.76|) . are calculated to yield 

{AluDn = {AIuI}d = {AIK}d = (1.82) 
{Af, n^jo = SijS"'' - (-£>r)G"'^6»f = dijS"^ + D^^G^'^df =: T^f (1.83) 

It is clear that the Dirac brackets (ll.82p must vanish since, by construction, the Dirac bracket 
of a constraint with any phase space function vanishes, recall Eq. (jl.36ll . In particular, the 
Gauss law functional = cp2 has vanishing Dirac bracket with any function ©[^j H], 

{a",O% = 0, (1.84) 

as one may easily convince oneself, cf. Eq. (I1.36p . All second-class constraints equations hold 
strongly. The non-symmetric matrix T^J' occurring in Eq. (|1.83l) can be understood in some 
sense as a non-abelian generalization of the transverse projector ijjEl It obviously has the 
projector property = T, and the projector orthogonal to it is given by L := 1 — T. We 
thus have 

= -b^^G'^'^df , T^f = 6f^ + D^^G^'^df , TL = LT = 0. (1.85) 
Just like for tij, it is important to note that T^!' has zero modes, 

dfTj;^ = T^^bf = . (1.86) 

We are now ready to pass to the quantum theory with the prescription 

i{Al,nl]n^[Al,Ill]. (1.87) 

It is understood that all fields are operators henceforth. The second-class constraints (|1.75p 
now hold strongly and can be imposed as operator identities, 

= = e?"^ = 4" = . (1.88) 

^''On top of that, zero modes of the matrix infer that G exists only in subspace of r^. This will be 
discussed in chapter [2] 

^'Letting D ^ d, one regains the familiar expression l|1.57p . Note the sign in G = —{Dd)^^. 
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1.4. Constrained quantization of Yang-Mills theory 



In order to arrive at a representation for the momentum operator 11^ that is to obey 

[Ai'^,U'^]=zT^^ , (1.89) 
it is helpful to split it into the orthogonal projections (from the right side) of T and L, 

^^ = ^^ + ^^^ uI'' = u'^t^^, uj:'' = u'^L%. (1.90) 

The commutation relations (I1.89P are thus decomposed into two separate ones, 

= (1.91) 

[A^,nf] = iT,f (1.92) 

and a situation similar to the abelian case, cf. section [L3l occurs. Recall that the longitudinal 
component ITH obeyed a trivial commutation relation equivalent to in (|1.91|) and had 
to be determined by the (strong) Gauss law (|1.6ip . The transverse component H-*- obeyed 
[A^,nj-] = tij which has the easily found solution (ll.63p . Here, the "non-abelian projector" 
T p.85|l in [A-*-,!!^] = iT contains field operators A-^. Nevertheless, a representation ofH^"' 
can be found to be 



as verified by plugging Eq. (I1.93P into Eq. p.92|) . Note here that T is not symmetric. 
The component is determined by the Gauss law p.68|l which yields 

Dful = DfTf^nf^ + Dfll'-L% = nf = -gpl,, . (1.94) 

where the contribution of vanishes in view of Eq. p.86|) . To solve Eq. (|1.94p for 11^, one 
may use the Helmholtz theorem to write 11^ as the gradient field of some field operator <I> since 

(a X n^)^ = e,,fc9f nf = -e^^udfnlDtC'^df = . (1.95) 

We can then write li^" = -df^'' and solve the Gauss law (fL94l) by = -^G^Vext- The 
operator then is explicitly given as a functional of field operators A-^, 

Il^,'^[A^]=gdtG'^[A^]pl^,. (1.96) 

One may subsequently verify the commutation relation (ll.9ip . 

Before we go on to derive the Hamiltonian operator from the classical one (|1.67p by promoting 
the fields and to their operator-valued counterparts, operator ordering is briefly dis- 
cussed, cf. Ref. [21]. Gauge fixing gives rise to non-trivial factors in the Hamiltonian, as can 
be seen in the path integral formalism. Using the the Faddeev-Popov trick [i2], the kinetic 
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1. Yang-Mills theory as a constrained dynamical system 



energy Ei^ yieldfH 

\ I vA^j[A^]^*[A^] (n^'^^ + n^'^^) j[A^]{ul^ + n^)) ^a^] , 



(1.98) 



where the state norm (^| ^) was chosen such as to cancel the group volume, and 

J[A^]=Det{G-^)[A^] (1.99) 

is the Faddeev-Popov determinant. The last line in Eq. (I1.98P is obtained by taking the 
hermitian adjoint of the product ^7(11^ + H^) with respect to the inner product in functional 
space when bringing the ^* to the left. One can read off from Eq. (|1.98|1 the form of the kinetic 
part of the Hamiltonian operator 

Hk[A^,u^,u^] = ^J-\A^] (n^'^^ + n^'^^) j[A^] (n^" + n^«) . (i.ioo) 

The decomposition of the momentum operator IT = 11^ + with the expressions (|1.93p 
and (|1.96|) can now be inserted into (ll.lOOj) . giving rise to four terms. Let us consider these 
separately. One term reads 

= ^Ipl^^F^-pl^^ (1.101) 

and is recognized as the Coulomb interaction of external charges via the Coulomb Green 
function F := G{-d'^)G. Using the relation f^l^T^^df = -f^^Lf^df as well as the identification 



^^The result is reminiscent of the Laplace-Beltrami operator A with a non-trivial metric g^i, and its deter- 
minant g — det{gfj,„), see Ref. [43], 

A = --^a^Vff3'"'a. , (1.97) 

as pointed out in Refs. |24l I44| . The Weyl gauge operators are then in Cartesian coordinates, whereas the 
Coulomb gauge operators correspond to curvilinear ones. In contrast to Ref. |24| . we here have no Cartesian 
coordinates to start with, since the temporal Coulomb gauge was already fixed on the classical level. 
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1.4. Constrained quantization of Yang-Mills theory 



of the dynamical gluonic charge pdyn h 

Dfui" = ni'bf = gAt'ui' =: gp'^^^^ (1.102) 



we further find 



= ^-pt.tG'^'{-dt^df)G'^''''Di''ui^ 

= ^ pl^^F^'^' pI^ (1.103) 



and 



2 J 



m ^mk^ 90k ^ Pext 



5 1 TT±f) Abb' /~ih' a' pM,' a pfiuc /-icd d 

— 2 'j '^k '-'k ^ Pext 

= y ;^Pdyn^i^%fxt , (1-104) 

two contributions that account for the Coulomb interaction of dynamical and external charges. 
Using the identity CA^Xn' = Ct' +Cn.L'r^kLXn' we get two further contributions 
to Hk, 

ttTT _ 1 1 rr^at ^YfTa 

^k - 2 'j ^ 

2. _J_TT-Lfe^^ba rjrpCajrJ-C 

2 jT' mk'J -'^ nk'-^n 

11 TT-La 'T'TT-'-a i ^ ^ J>' ^r<b'a ( '^a'af>aa"\ /~ia" c'c' 

— 2 J ^ ~2 'jP'^y^^ \~^k (^k Pdyn 

= 1 In^ jn^« + ^ jF^'^V^;„ , (1.105) 

namely the kinetic energy of transverse gluons as well as the Coulomb interaction of dynam- 
ical charges. The magnetic potential operator Hp can be obtained straightforwardly from its 
classical counterpart in Eq. p.67|) since it only comprises the coordinates Af^. It simply yields 



Hp = -Bl[A^]Bt[A^]. (1.106) 



1 

2^ 

''With the conventions used here, pi = "''n^ * = -A-^"'' ■ n*. Thus, one arrives at dkU-l = -^(Pext + 



Pa 
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1. Yang-Mills theory as a constrained dynamical system 



Adding the contributions (ll.lOip . (ll.KK^Il . (I1.1()4D . (11.10511 and (ILlOfip . we finally arrive at the 
Hamiltonian operator in temporal Coulomb gauge. In explicit notation, it reads 



+ l_ I d\d'yjp'^{x)JF'^\x,y)p\y) . 



(1.107) 



Here, we have defined the total charge density 

p'^{x) := pl^^{x) + pl,,{x) = At{x)Ill{x) + i;\x)tT'^i^{x) . (1.108) 
For future reference, let us write down the explicit form of the Coulomb operator F, 

F''\x, y)= j d^x' ^^[x, x') (-^2) G^\x', y) (1.109) 

and of the inverse ghost operator 

[G-T{x.y) = -5'^'dH\x,y)-gAt{x)dt5\x,y) . (1.110) 



Christ and Lee ^4] derived the same result as in Eq. (|1.107p using a different approach. They 
quantized the Cartesian coordinates in Weyl gauge and only subsequently transformed into 
the Coulomb gauge, keeping the Gauss law as a constraint on the wave functional. Here, 
we were able to reinforce the Gauss law as a operator identity, having used Dirac's concept 
of constrained quantization. It is reassuring to see that the two methods produce the same 
result for the Yang-Mills energy spectrum, keeping in mind that in general such an equivalence 
cannot be established [M B5l BGlHTl HH] . 



30 



2. Infrared ghost dominance 



2 Infrared ghost dominance 

The complex structure of the Yang-Mills Hamiltonian (|1.107p derived in the previous section 
does not allow for a rigorous calculation of the energy spectrum or the Green functions. Per- 
turbation theory in the gauge coupling g provides a possibility to describe interactions with 
high momentum transfer in QCD where asymptotic freedom guarantees that g can be treated 
as a small parameter. The infrared sector, however, is not accessible by perturbative methods 
since we know from the renormalization group that g increases as we go to lower energies. An 
understanding of long-range interactions of quarks and gluons thus calls for nonperturbative 
methods. Color confinement, the experimental evidence that quarks and gluons have never 
been detected as asymptotic states and must be confined into color singlets, is to date one of 
the most challenging topics in theoretical physics. 

By means of lattice calculations, it has been possible to penetrate the infrared nonperturbative 
sector of QCD and uncover a confining potential between (static) quarks [Ml [8]. At present, 
however, available lattice sizes do not suffice to describe the Green functions in the deep 
infrared [50l EI]- Despite the ever-growing computational power, a numerical simulation of 
confinement on the lattice will always fall short of proving it. Nevertheless, the achievements 
of the lattice community with calculations that are particularly reliable in the intermediate 
momentum regime, are indispensable for the mutual dialogue with the continuum approach. 
A joint theoretical investigation of QCD phenomena is a promising project for progressive 
research. 

The continuum approach has the intriguing feature that the asymptotic infrared limit can be 
studied analytically. In the last decade, a new understanding of infrared QCD has arisen from 
studying continuum Yang-Mills (YM) theory via Dyson-Schwinger equations. The Landau 
gauge has the advantage of being covariant, allowing for straightforward perturbative calcu- 
lations, and it therefore encouraged many to intensive investigation of the infrared properties 
of YM theory [52l [53l [S^. In Coulomb gauge, non-covariance brings about severe technical 
difficulties which are only recently on the verge of being overcome [SSllSHllS^. Nevertheless, 
the Coulomb gauge might be the more efficient choice to identify the nonabelian degrees of 
freedom. It is well-known that screening and anti-screening contributions to the interquark 
potential are neatly separated in Coulomb gauge perturbation theory |[58]. As for the in- 
frared domain, the Gribov-Zwanziger scenario serves as a transparent confinement mechanism 
[591 EQ]. 

In this chapter, the Gribov-Zwanziger scenario shall be introduced and analyzed in the tempo- 
ral Coulomb gauge. It will be thus possible to provide information on the asymptotic infrared 
behavior of the Green functions and recover a linearly rising potential between static color 
charges, i.e. heavy quark confinement. A full numerical calculation in chapter [3] will be seen to 
reproduce the infrared asymptotics given here analytically. Moreover, the infrared asymptotics 
of Landau gauge can be retrieved from a generalization of the results in the temporal Coulomb 
gauge. 
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2.1 Gribov-Zwanziger scenario of confinement 

Having outlined the issues of redundant degrees of freedom in section [L2l a discussion of gauge 
fixing in the context of Yang-Mills theory may follow naturally. In particular, we have touched 
upon the issue of uniqueness of a gauge fixing condition condition 

X"[^] = 0. (2.1) 

Uniqueness was shown to be locally guaranteed by the criterion (|1.41|) which reads for the 
classical YM theory 

Det({x^an) 7^ 0. (2.2) 

Locally means here that only a small neighborhood of the point A in configuration space that 
satisfies (|2.1|) is considered. An infinitesimal gauge transformation according to Eq. p.29|l . 

xn^]-x"[^] + {x^an , e\ (2.3) 

A=A 



shows that A is the unique solution to (|2.1|) only if Eq. (12. 2p holds. Let us now turn to the 
quantum theory in the Schrodinger picture. Gribov f59] discovered that if the Coulomb gauge 
condition x"'{x) = dkA1^{x) = is chosen, uniqueness of its solution A is lost if the matrix 
with components 



[G-Yi^,y) = [-d'6'''-gAf>d^j6Hx,y) (2.4) 

develops zero modes. From the explicit form of this is seen to happen for large magnitudes 
of gA, i.e. in the nonperturbative domain. The field configuration A can then shown to be 
connected by a gauge transformation to one or more A^ , called Gribov copies, that also obey 
the Coulomb gauge condition. Gribov copies lead to incorrect results in the calculation of 
expectation values due to over-counting of physically equivalent field configurations. The 
criterion for avoiding Gribov copies is the quantum analogue of Eq. (12. 2p . noting that G~^ = 
{dkAk,G}, and reads 

J[A] = Det (G-^) / . (2.5) 

It is illuminating to see how the above uniqueness condition follows from a direct calculation. 
Consider a gauge transformation A^ where U{x) = e~"^'^^ 1 — a{x) is very close to unity. 
The field A'^ is a Gribov copy of A if 

= dkA'i! = dkUAkW + -dkUdkU^ 

9 

= u(- (dka) Ak + Ak (dka) - - (dka) (O^a) + - {d^a) + - {dta) (dka)) 
V 9 9 9 J 

= ^U{{d^a)+g[Ak,idka)])U^ 

= i[/r"((rt2 + 5if4)a^);7t. (2.6) 
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A 








1 J 




^ x = o 



Figure 2.1: Configuration space and its restriction to the gauge fixing 
condition X = 0, the (first) Gribov region J7 and the fundamental mod- 
ular region A. 



Since the matrices U are regular and are linearly independent, Eq. (|2.6|1 can be written as 

J d^y[G-^y'{x,y)a\y)=0. (2.7) 

Thus, if the condition (|2.5|1 is fulfilled, then a°'{x) = and = A, and the solution A of the 
gauge fixing condition (|2.ip is unique. 

The configuration space is therefore decomposed into the so-called Gribov regions, alternating 
in sign of the Faddeev-Popov determinant J. By definition, the first Gribov region Q contains 
the perturbative vacuum, gA = 0, and has a positive definite Faddeev-Popov determinant, 

n = {A\ dkAk = , J[A] > 0} . (2.8) 

In order to avoid Gribov copies, it is necessary to restrict the configuration space to the first 
Gribov region 0, bounded by the Gribov horizon d^l, see Fig. 12.11 Practically, it is not 
straightforward to implement this restriction. One might naturally ask whether a gauge fixing 
condition exists, different from the Coulomb gauge, for which the Faddeev-Popov determinant 
is manifestly positive and the Gribov problem does not arise. However, a theorem due to Singer 
[6T] states that (on a three- or four-dimensional compact manifold) there exists no gauge fixing 
condition that does not have the Gribov problem. 

A restriction to can be realized by an action principle, as noted by Polyako\El- It requires 
that the norm of the gauge field Ai:{x) be minimal along its orbit, parametrised by the 

^See footnote in Ref. |62j . 



33 



2. Infrared ghost dominance 



gauge transformations U{x) = e 

Fa[U] := ld^xtT{A'i{x)y 

= Fa[1]-- [ d^x tr {a{x)dkAk{x)) + 
9 J 

^ j d^xd^ytr(a^{x)G-^[A]{x,y)a{y)^ . (2.9) 

The stationarity condition sets A = A and the minimum condition yields J'[A] > 0. Action 
principles are a common technique to implement gauge fixing on the lattice [63] , 



However, the restriction to O is not sufficient to exclude gauge copies. Various examples 
[Ml [65] indicate that there are Gribov copies inside $7. A proof of existence was given by Ref. 
|66| . see also Ref. [67]. These remnant Gribov copies show that the minima of -Fa[C^] in Eq. 
(|2.9p are degenerate. Lattice calculations confirm that there are indeed Gribov copies inside 
[Ml EHl EO]- A further restriction of configuration space is needed to completely exclude 
Gribov copies. The fundamental modular region A is defined to be the set of uniquqj absolute 
minima of the functional (12. 9|) among gauge orbits, 

A := : Fa[1] < Fa[U] VC/} . (2.10) 



It was shown that the interior of A is indeed free of Gribov copies whereas the boundary OA 
still contains Gribov copies [71]. The gauge condition specified by Eq. (I2.10p is also termed 
"minimal Coulomb gauge" [72]. Although conceptually vital, the fundamental modular region 
as yet lacks the utility for explicit calculations. Little is known about the boundary dA. 
Fortunately, it was shown in Ref. [73] by means of stochastic quantization that despite the 
Gribov copies inside O, expectation values are computed correctly if the integration domain 
is rj instead of A. This is due to the finding [^ that the dominant field configurations lie on 
dQ n dA, the common boundary of O and A|f| We shall therefore restrict the configuration 
space to Q in the following. Gauge fields A are always transverse and the superscript _L is 
abandoned for brevity. 

The restriction to a compact region in configuration space, be it or A, has important conse- 
quences for the infrared sector of the theory. Due to asymptotic freedom, the ultraviolet is not 
affected by the horizon since when the coupling becomes small, all relevant configurations are 
in the vicinity of gA = 0. It was shown that this point always has a finite distance from the 
horizon [77]. The nonperturbative infrared sector, on the other hand, can in principle be gov- 
erned by any region within From statistical mechanics, we know that in a compact sphere 
with radius r of high dimension A^, the probability distribution is concentrated at the bound- 
ary, due to the "entropy factor" r^~^dr. If Yang-Mills theory is regularized on the lattice, 

^Up to global gauge transformations. 

^In this context, let us mention that the field configurations on dflCidA can be identified as center vortices. 
These field configurations were found to drive the confining mechanism on the lattice in Coulomb gauge [75] 
as well as in Landau gauge [76j. 
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Zwanziger argued |60] that the situation is comparable to the latter example, the dimensional- 
ity N even diverges in the continuum limit. In the asymptotic infrared, the momentum space 
ghost propagator 

Dcik) := ^^2^5"'/ d^x {G[Ar'{x,y)e-'''<'--y^ , k = \k\ (2.11) 

is therefore expected to diverge strongly, 

lim[k'^DG{k)]~^ = , (2.12) 

k — ^0 

due to divergent eigenvalues of the operator G on the Gribov horizon. Eq. (|2.12|) is known 
as the horizon condition. Gribov discussed perturbatively in Ref. [59] that only for A; = the 
ghost propagator Dcik) can have a pole. Zwanziger's horizon condition (|2.12|) . on the other 
hand, is a nonperturbative statement. Since the Coulomb potential Vc{k) between external 
color charges is proportional to 

J d^x (G(-a2)G>"^ {x, y) e-^H^-y) , (2.13) 

the infrared enhancement of DG{k) from effects on the Gribov horizon might be the driving 
mechanism for Vcik) to diverge as k~'^ and for a linear confining potential between heavy 
quarks to emerge. This notion is known as the Gribov-Zwanziger scenario of confinement in 
the Coulomb gauge. 



2.2 Stochastic vacuum in Coulomb gauge 

Vacuum expectation values that are solely dominated by entropy in configuration space are 
intriguingly simple. According to the Gribov-Zwanziger scenario, infrared correlation func- 
tions of Coulomb gauge YM theory fall under this category and hence can asymptotically 
be described by a vacuum wave functional ^o[^] stochastically distributed in configuration 
space. Whilst the full structure of a vacuum expectation value, in particular the perturbative 
regime, is sensitive to a non-trivial probability distribution the infrared limit is fully 

determined by the geometry in phase space, and one may write ^'o[^] = 1 to find correct 
results |78| . This choice of the wave functional is here referred to as the stochastic vacuum and 
suggests the following schematic prescription for the asymptotic evaluation of Coulomb gauge 
expectation values of any operator 0{k), depending on momentum k, 

m Oik) 1^) = J VAJ^*0{k)^ ''^^ J ^^^^ = (^ol O \^o) . (2.14) 

n n 

We introduced a factor Vol{i^) = f^'DAJ'[A] to normalize the wave functional, (^'o| ^o) = 
1, which is possible in the compact region Q. A legitimate question to ask is whether the 
factor of the Faddeev-Popov determinant J' does not spoil the argument that the probability 
distribution is strongly enhanced at the Gribov horizon. Recall that the Gribov horizon di} 
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was defined as the boundary of 0, where the Faddeev-Popov matrix G~ picks up vanishing 
eigenvalues. Hence, the factor J' = DetG~^ will suppress the probability distribution on dQ. 
Nevertheless, it is not clear a priori whether this suppression is strong enough to spoil the 
enhancement of entropy and we shall put forth as a working hypothesis that it does not|l| 
On the lattice, the density of eigenmodes of the Faddeev-Popov operator is under current 
investigation |80| . 

Although the prescription (|2.14|) for infrared correlations can be motivated from the Gribov- 
Zwanziger scenario, it cannot be taken for granted that the solutions thus obtained are mathe- 
matically unique on the one hand, or physical on the other. To further contemplate Eq. (|2.14l) . 
note that ^'o[^] = 1 is the correct wave functional in 1 + 1 dimensional Yang-Mills theory 
(EH [82]. However, the entropy argument is lacking in 1 + 1 dimensions where the gauge field 
becomes a compact quantum mechanical variable. In 3 + 1 dimensions, the dynamics might 
be substantially different. After all, we will try to find the vacuum wave functional in 3 + 1 
dimensions that describes not only the infrared limit but also ultraviolet correlations. In the 
Hamiltonian formalism, the Ritz-Rayleigh variational method gives the possibility to introduce 
a variational parameter into the wave functional and determine it by minimizing the energy. 
This method will be discussed in chapter [3] to discuss the full momentum dependence of the 
Green functions. It will turn out that the full solutions agree with the solution in the stochas- 
tic vacuum for asymptotically small momenta [83j [HH [85j [86] . Turning the argument around, 
these findings support the Gribov-Zwanziger scenario and the notion of a stochastically dis- 
tributed probability distribution in the infrared. The investigation of the stochastic 
vacuum pursued in this chapter is quite different from a variational approach in that it lacks 
equations of motion arising from minimizing the energy. However, we may extract information 
from the geometry in configuration space. In the following section, it will be shown how the 
boundary conditions of the Gribov region give rise to Dyson-Schwinger integral equations that 
determine the properties of the Green functions. 



2.3 Dyson-Schwinger equations 



A convenient tool for the evaluation of expectation values are generating functionals. In par- 
ticular, the Dyson-Schwinger equations that intercorrelate all Green functions can be straight- 
forwardly derived from the path integral representation of a generating functional. The clas- 
sification of Green functions into full, connected and proper Green functions is common and 
the same applies for generating functionals [6]. Full Green functions corresponding to the 
stochastic vacuum ^'o can be found by 



{0[A]) = (M/qI 0[A] \^ 



1 



0/ 



Z[j] 



o 



6_ 



(2.15) 



j=0 



''As a simple example [79], consider a d-dimensional sphere with < r < 1 and let the probability distribution 
r) be damped at the boundary by a factor (1 — r)" with order n{d). The entropy factor from the integral 



measure, drr"^ ^, enhances the probability at the boundary. For large dimensions d, 
then sharply peaked at r = 1 — provided ^ ^ 1, despite the damping factor. 



-1(1 
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if we define the generating functional Z[j] of full Green functions by 

Z[j] = J VAJ[A] e^-^ , j-A^ j SxA'i{x)Uk{x)jl{x) . (2.16) 

It is easily verified that Eq. (|2.15|) indeed reproduces the expectation values in the stochastic 
vacuum, see Eq. (|2.14l) . Dyson-Schwinger equations follow from a quite simple statement, 
namely that the (path) integral of a total derivative can be expressed by boundary terms. 
This idea is applied to the path integration in configuration space, bounded by the Gribov 
horizon. The gauge-fixing procedure endowed us the Faddeev-Popov determinant J7 in the 
path integral, see Eq. (|2.4p . One is thus led to the identity 

= / ^^ ^^[^] -"-^ = I ^^^[^] (^^^ + tMy)jLiy?j (2.17) 

since at the Gribov horizon the Faddeev-Popov determinant vanishes by definition, ^'[AJIgQ = 
0. Applying to Eq. (I2.17P a functional derivative 6/5jf{x), dividing by Z[j] and setting sources 
to zero, j = 0, gives 

t#(.,y) = -(^?M^^). (2.18) 

In order to make sense of the above expression, first note thajfl [83] 



"°^M-T^.^'"°:''-^'--TvMr°->)), (2.19) 



SA^{x) SAf{x) 
where T^'°'{x) is the tree-level ghost-gluon vertex, here matrix- valued with components 



{Ty{x))'\y,z) = -^-^^^^^r^ = 9{f^)\,{^)5H^,y)dj!S\y,z) . (2.20) 



Second, plugging (|2.19|) into (|2.18p . the object {AG) can be identified as the dressed connected 
ghost-gluon vertex that may be decomposed b}l§ 

{A1[x)G[A]) = fd'y{An^)A]{y))l^^^ 

= J [A1{x)A]{y)) (G[A]Tf {y)G[A]) 

= j d^yDl'{x,y)DGT'j{y)DG (2.21) 
into the proper ghost-gluon vertex F^ with propagators attached, namely the gluon propagator 

DtHx,y):=(A^{x)A'{y)) , (2.22) 



'it is convenient to denote by O a matrix with elements 0°'''{x, y) — {x, aj O \y,h) in color and (continuous) 
coordinate space. The trace "Tr" sums up diagonal elements in both the latter spaces. For notation, see also 
appendix [Al 

^If the generating functional is Gaussian, Wick's theorem infers the factorization in the first line of Eq. 
I|2.2ip . Without the use of Wick's theorem, as for the stochastic vacuum, one may use functional methods |73| . 
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and the (here matrix- valued) ghost propagator Dq-, with matrix elements 

D''^{x,y) = lG[Ar\x,y)) . (2.23) 



The present calculations circumvent the common introduction of ghost fields by interpreting 
J[A\ as a Gaussian integral over Grassmann variables. In fact, the introduction of Grassmann 
fields and their sources spoils the argument that boundary terms from the Gribov horizon 
vanish [73]. It is therefore advisory to stick to the functional determinant J[A\ and not to 
introduce ghost sources. Nevertheless, we refer to Dq as the "ghost" propagator and as the 
"ghost-gluon" vertex. 

With the above definitions, the Dyson-Schwinger equation (12.18(1 becomes 
tf^{x,y) = TTUA'i{x)G[A])Tf{y) 



= Jd'zDt^{x,z)Tr(^DGTl{z)DGr'.'\y)) . (2.24) 

It would be pleasant to invert the gluon propagator in Eq. (|2.24|) . but the transverse projector 
tij is singular. Yet, without loss of generality, we can define the matrix Da by D^^{x^y) = 
tij{x)D'^{x,y) and its inverse, D^^Da = 1. The translationally invariant function 

[DT^T'i^^y) ■■=t^,i^) [DATi^^y) (2-25) 

is referred to as the inverse gluon propagator0 One can now deduce from Eq. (|2.24|) the 
coordinate space gluon DSE in its most useful form, 

[DT.'Y\x,y) = Tr (DGr?(^)I?Gr°'^y)) . (2.27) 



In Fig. 12.21 the gluon DSE (|2.27|) is interpreted diagrammatically. By the loop diagram on the 
r.h.s., the so-called ghost loop, the gluon propagator Dij is correlated with the ghost propagator 
Dg- The latter needs to be determined by its own DSE. To derive it [83|, first note that 

= - I d'xAtix)T'^\x) . (2.28) 

Multiplying the above equation with Go from the left and G from the right gives 

G = Go + Go y d^x AI{x)tI'''{x)G . (2.29) 

Taking the expectation value of (|2.29|) and using the identity (|2.2ip for the connected ghost- 
gluon vertex yields 

Dg = Go + Go^Dg (2.30) 



''Note that D^^^ is inverse to Dij only in the transverse subspace of the full Lorentz space, 

I d'z [D-,'] {X, z) [D^,f {z, y) = f^ix, y) . (2.26) 
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Figure 2.2: Tie set of Dyson-Schwinger equations in the stochastic vac- 
uum. Dashed lines represent interacting connected ghost propagators 
whereas the dotted line is the tree-level ghost propagator. Curly lines 
represent the gluon propagator. Empty blobs stand for proper ghost- 
gluon vertices and dots represent tree-level ghost-gluon vertices. 



A multiplication with Gq from the left and Dq from the right gives the final form of the 
ghost DSE in coordinate space, 



shown in Fig. 12. 2[ Let us emphasize that Eq. (|2.32l) was not derived from the path integral, in 
contrast to the gluon DSE (I2.27p . The "operator method" using Eq. (I2.29P offers an alternative 
to the conventional path integral method, without the introduction of ghost sources. Restric- 
tion to the Gribov region Q is ensured by cutting off the path integral at dQ, the operator 
method formally also accounts for this restriction since the operator G in Eq. (|2.29p is not 
defined on dQ. 

Having at hand the gluon DSE (|2.27p and the ghost DSE (|2.32|) shown in Fig. 12.21 one may 
now try and solve for the two-point Green functions ("propagators") Dq and Dij. However, 
the proper ghost-gluon vertex that appears in both DSEs couples the propagators to higher 
n-point functions, thus giving rise to the typical infinite tower of coupled integral equations. 
After establishing a connection to the Landau gauge, an approximation for the ghost-gluon 
vertex will be motivated in section [231 which makes a solution for Dq and Dij feasible. 

2.4 Connection to Landau gauge 

Linear covariant gauges are a natural choice for perturbative calculations in the Lagrangian 
formalism. With Lorentz invariance being manifest and after Wick rotation to Euclidean 
space, perturbative loop integrals can be computed with standard methods jSTj. However, this 
is at the cost of losing a positive definite metric in Hilbert space. One cannot simultaneously 



with the ghost self energy 




d^ix'y'uv] [r°'"(tx; 



,)] {x, x')Dff{u, v)Di'' {x', y') [V^^iv)] "'"{y', y) . (2.31) 



[D-^T^^^y) = [G,Ti^^y)-^''\^^y) 



(2.32) 
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maintain locality, positivity and manifest Lorentz invariance [8H|- A Hamiltonian approach 
is thus obscured and the potential between static quarks as it occurs in the Coulomb gauge 
Hamiltonian is not accessible in the Landau gauge. The confinement problem must therefore 
be studied by alternative means, e.g. the Kugo-Ojima confinement criterion [89l [90] . 

Usually, the class of linear covariant gauge conditions x"(^) = d'^A'Jj^{x) — B"-{x) 0, is 
implemented in the path integral by means of the Faddeev-Popov method. The arbitrary 
function is integrated out with a Gaussian weighting of width ^, which effectively smears 

of the Lorenz gauge condition, 

5{d^'A1) ^Q-U^^"^^^^ (2.33) 

controlled by an arbitrary parameter ^. Apart from the Faddeev-Popov determinant J = 
DetTW, with 

5{^^'A^l{x)) 

the YM Lagrangian is thus supplemented by a gauge fixing term Cgj = {9^^^) ■ The 
Landau gauge is then obtained in the limit ^ — > 0. From gauge invariance, a Slavnov-Taylor 
identity can be derived [91] that states that the longitudinal part of the gluon propagator 
remains unchanged by radiative corrections, 

t^''{k)D^,{k)=i. (2.35) 

It is important to realize that during calculations, ^ is kept non-zero, so that the gluon prop- 
agator has an inverse. The fields ^JJ(2;) are therefore not transverse, only the Landau gauge 
correlation functions are. In particular, the ghost-gluon vertex in Landau gauge does not have 
a transverse projector attached to its gluon leg, since the transversality is imposed "off-shell". 
This is in contrast to the constrained quantization in temporal Coulomb gauge, where the 
transversality condition is enforced "on-shell" and the operators A^{x) are manifestly trans- 
verselfl 

The Gribov problem affects the configuration space in Landau gauge just like it does for 
Coulomb gauge. In fact, the Faddeev-Popov gauge fixing procedure in Landau gauge gives 
rise to the Faddeev-Popov determinant 



J[A] = Det (^-5"" 9^,5^ - gA''J>^^'J (2.36) 

and one cannot miss to acknowledge the similarity to the one in Coulomb gauge. It was pointed 
out by Gribov in his seminal paper [52] that if Wick rotation is employed, the Landau gauge 
condition becomes equivalent to the Coulomb gauge condition. The dimension D = 3 -|- 1 
of Minkowski spacetime infers that d = 4 after Wick rotation, as opposed to d = 3 in the 
temporal Coulomb gauge. 



The terminology using "on-shell" and "off-shell" to describe the gauge-fixing technique is encountered, e.g.. 
in Ref. [73]. 
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Recall the Gribov-Zwanziger scenario in the Coulomb gauge. The entropy argument for in- 
frared dominant field configurations at the Gribov horizon along with the horizon condition 
(|2.12p for the ghost propagator applies to the Landau gauge in Euclidean space just the same 
One may therefore expect that the infrared limit in Landau gauge yields the same set of 
integral equations as for the Coulomb gauge, with the dimension shifted to d = 4. To see how 
it comes about, consider the expression for a Landau gauge-fixed expectation value. 



where the weight of field configurations is determined by the Euclidean Yang-Mills action. 
Shifting the field variable hy A ^ gA, the weight acquires the form 



and one recognizes that to zeroth order of the strong coupling expansion in ^ <C 1, the 
evaluation of expectation values depends only on the Faddeev-Popov determinant [73] and is 
therefore equivalent to the stochastic vacuum (I2.14p . with the difference that the Faddeev- 
Popov determinant in (|2.36ll is to be computed in d = 4 dimensional Euclidean Lorentz space. 

The full structure of Landau gauge Dyson-Schwinger equations was studied intensively in the 
past decade [92l|5lj. It turned out that the infrared sector of the solutions is asymptotically 
determined by the set of equations shown in Fig. 12.21 Apparently, the strong coupling expansion 
in the Landau gauge Lagrangian formalism accounts for the correct infrared physics, while in 
the Coulomb gauge Hamiltonian formalism the stochastic vacuum is a simple but sufficient 
state to describe infrared correlations. The approximation needed to arrive at a solution to the 
coupled Dyson-Schwinger equations shown in Fig. 12.21 concerns the ghost-gluon vertex. This 
will be the topic of the next section. 

2.5 Ghost-gluon vertex 

In gauges where the gluon propagation is transverse, such as Coulomb and Landau gauge, the 
ghost-gluon vertex has the feature of nonrenormalization. In this context, one can motivate 
that the vertex is approximately tree-level for all momentum configurations. The solutions to 
all Green functions will qualitatively depend on the simple structure of the ghost-gluon vertex 
and therefore this section is devoted to assess the tree-level approximation. 

As far as renormalization is concerned, the ghost-gluon vertex in Coulomb or Landau gauge 
is quite an unusual Green function. Whereas one usually encounters ultraviolet divergences 
in a perturbative expansion of a Green function, the ghost-gluon vertex appears to be a fi- 
nite function, to all orders of perturbation theory. This can be understood by considering 
a Dyson-Schwinger equation for this vertex. In the Landau gauge, the ghost-gluon vertex 
Dyson-Schwinger equation was derived in Ref. [53] and can be represented in momentum 
space as shown in Fig. 12.31 The proper ghost-gluon vertex r^(A;; q,p) is related to its tree-level 




(2.37) 



^- f d'^XBCYAlixE) 



■\ J d'^XECvMixE) 



(2.38) 
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Figure 2.3: Dyson-Schwinger equation for the ghost-gluon vertex. The 
l.h.s. represents the proper ghost-gluon vertex, with dehnitions of the 
momenta k, q and p = q + k. The filled blob on the r.h.s. corresponds to 
a four-point function, see text. 

counterpart T^{q) and a connected four-point function rfj,y{k,i;q,p) that will be explicitly 
defined in section 13. 51 Here, it shall suffice to recognize the generic structure of the loop 
diagram, 

T^{k; q,p) = ro (g) + I (fer^'^ip - i)D,,{l)DG{l - v)Tu^.{k, -i; £-p,q). (2.39) 

Taylor anticipated in Ref. [9^ the structure of the Dyson-Schwinger equation (|2.39p and 
perceived that the contribution of T^{p — i) to the loop integration can be factored out due 
to the transversality of the gluon propagator. In momentum space, with the color structure 
suppressed, the contraction of the gluon propagator D^^, with the tree-level ghost-gluon vertex, 
= igQfj,, yields 

r°(p - £)D,,ii) = T°,{p)D,^i£) . (2.40) 

Hence, the degree of divergence of the ^-integration is diminished. Landau gauge perturbation 
theory explicitly confirms that the divergent part of the contribution to the vertex is pro- 
portional to ^, i.e. it vanishes in the Landau gauge limit, see e.g. [95]. Taylor further argued 
that if the loop diagram does not develop any poles in the infrared, one can take the infrared 
limit of the incoming ghost momentum p — > and find from Eq. (|2.39p 

lim r^{k;q,p)=rl{q) . (2.41) 

For the labelling of momenta, see Fig. 12.31 Since one particular momentum configuration, given 
by Eq. (|2.4ip . is found where the ghost-gluon vertex is finite, its multiplicative renormalization 
constant Zi, being momentum independent, has to be finite, Zi < oo. In an appropriate 
renormalization scheme, Zi can defined to be trivial, Zi = 1. Multiplicative renormalization 
of the bare ghost-gluon vertex F^, understood here as a vector, is defined by 

rB{k,q,p;A) = Zi{p,g,i,K)T{k,q,p;p) (2.42) 
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Figure 2.4: Projection of the ghost-gluon vertex in 3 + 1 dimensional 
Landau gauge as a function of one momentum, taken from Ref. [93 ] (for 
details, see there). Both the DSE result and the lattice data show that 
the ghost-gluon vertex is basically at tree-level. 

and relates Tb to the renormalized vertex T. Since in the Landau gauge hmit, ^ 0, one 
has Zi < oo, we can drop all dependences on the regulator A. The dependence of Zi on the 
dimensionful renormalization scale /x must drop out since g and Zi are both dimensionless. 
The number Zi is then defined by a renormalization prescription, i.e. a certain momentum 
configuration P := is chosen for which T{kr,qr,Pr) = r'^(A*). In order to avoid 

cuts in the complex plane, P is usually chosen in the perturbative regime. Nevertheless, 
the renormalization prescription with P = {/^^,/Li^,0} is a convenient one since Eq. (|2.41ll 
then infers that Zi = 1. Marciano and Pagels [96] pointed out that for the symmetric point 
P = with fj? / 0, one generally has Z\^\ (cf. [97]). The symmetric point with 

[i} = 0, however, does not necessarily give Z\ = \ since the infrared limits of momenta are not 
interchangeable, as argued below in section [331 In Ref. [98], the choice P = {0, /li^, /li^} was 
advocated. In the following, we will set Z\ = 1, but come back to it in chapter [H 

Due to the nonrenormalization of the ghost-gluon vertex and the related property (|2.41|) . recent 
Landau gauge DSE studies used a tree-level ghost-gluon vertex for all momentum configura- 
tions [92j[5l]. Of course, this approximation needs to be assessed further. It was discussed in 
Ref. [93] that with the tree-level vertex as a starting point, the iteration of the vertex DSE 
(|2.39p will not lead too far away from the tree-level value. The changes stay in the range of 
10% — 20%. Furthermore, lattice calculations [981 [991 riOO] confirmed that the vertex is basi- 
cally tree-level. In Fig. 12.41 the DSE results along with the lattice results for the vertex are 
shown. Thus, strong support is available that the perturbative arguments of Taylor, stated 
above, hold in the nonperturbative regime. 

In Coulomb gauge, the argumentation for a tree-level ghost-gluon vertex is similar. Although 
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the perturbative calculations in Coulomb gauge are far from trivial [55l [56l EI] , the DSE for 
the ghost-gluon vertex can be shown to reproduce the Landau gauge struture depicted in 
Fig. 12.31 155| . The gluon propagator is transverse as well and therefore Taylor's reasoning for 
nonrenormalization can also be applied to the Coulomb gauge. Within a set of interpolating 
gauges, the nonrenormalization of the ghost-gluon vertex can be established as well, considering 
the Coulomb gauge limit |101| . The Landau gauge calculations shown in Fig. 12.41 were also 
carried out in d = 3 spatial dimensions and the results were qualitatively the same [93j . 
Therefore, it seems to be an appropriate approximation to the Coulomb gauge DSEs to use a 
tree-level vertex instead of a proper vertex, analogously to the Landau gauge, 

T,{x) « Tlix) . (2.43) 



2.6 Analytical solution for propagators 

The approximation scheme (|2.43|) exhibited in the previous section enables us to find a solution 
to the coupled set of Dyson-Schwinger equations. With the Coulomb and Landau gauge 
differing in the number d of spatial dimensions only, it is convenient to keep d as an unspecified 
parameter. The gluon DSE J2.27|) can be transformed into momentum space defining the scalar 
gluon propagator Da(^) b}0 

6'^{k-p)6''^DA{k) = jd'^[xy\Df{x,y)e-'^''+'P-y. (2.44) 

Translational invariance allows us to write DA{k) as a function of one variable which makes 
momentum space calculations convenient. The color structure of propagators is chosen to be 
diagonal, i.e. Z)g(x,y) = DG{x,y)5°'^ . In the approximation scheme with a tree-level ghost- 
gluon vertex, it is verified straightforwardly in a perturbative expansion to all orders ("rainbow 
expansion") that color-diagonal propagators solve the set of integral equations. Using Eq. (|2.1ip 
for the definition of the momentum space ghost propagator Dcik), Fourier transformation of 
the gluon DSE (12.27P expressed in terms of Dr^{k) := tij{k)D~^^ {k) yields 

Dr^{k) = -N,j d''ir^,{k,e)T^j{k,e - k)DG{i)DG{£ - k) . (2.45) 

Here, the tree-level ghost-gluon vertex in momentum space is derived by Fourier transform 
from Eq. K^Hh . 

I d^ixyz] (r°'"(x))'"=(y, Z) ^ik■.-^g■y-,^p■^ ^ ^fa^c^d^^ _ ^ _ ^)pO(^^ ^) 

^Tl{k,q)=igtk,{k)qj . (2.46) 

The color trace in (|2.45p was computed using tr (f"f^) = -iVc5''^ see appendixEl Eq. (ITiBll 
is seen to be manifestly transverse. This follows from the fact that in the derivation of the 



®Note that in momentum space, the matrix notation of coordinate space is no longer used. Moreover, we 
simplify the notation henceforth by abandoning the bold-face typesetting for spatial vectors. Latin Lorentz 
indices now range from 1 to d. 
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DSEs in section 12.31 the Coulomb gauge condition was imposed "on-shell". We shall see in 
section [23] that the "off-shell" gauge condition leads to a (slightly) different gluon DSE. In Eq. 
(|2.45p . we can take the trace in Lorentz space to find, using tii{k) = d — 1 and Eq. (|2.46p and 
writing k = k/\k\, 

D-/{k) = j (l - {l-kf) DciDDGii - k) . (2.47) 

The ghost DSE (|2.32p can be Fourier transformed in the same manner. With a tree-level 
ghost-gluon vertex, it reads in momentum space 

D-\k) = A;2 _ g^Nck'^ j d'^e (l - {£ ■ kf'^ Da{1)Dg{£ - k) . (2.48) 

Eq. (|2.48|) can be put in a form that explicitly satisfies the horizon condition (|2.12|1 by subtract- 
mg k~'^D^^{k) at zero momentum. This subtraction does not introduce an extra parameter 
since by the horizon condition ()2.12p . {k"^ DG{k))~^\k=Q = 0, and we get 

{k'DGik))'^ = g^N.J d'^lDAii) (l - (i-kf) {DG{i) - DG{i - k)) (2.49) 

for the ghost DSE with an integral that converges more strongly in the ultraviolet than the one 
in Eq. (|2.48ll . In the ultraviolet sector, the horizon condition can be used to remove ultraviolet 
divergences in the ghost DSE. 

The set of equations (|2.47p and (|2.49p was first solved without approximations in Refs. [97lll02| 
by a power law ansatz. The derivation and the solution is reviewed here, and some subtleties 
are pointed out. 

With power law ansatze for the propagators, 

^^(^) = (P)^' ^^(^) = (FF^ (^-^^^ 

the integrals in (|2.47|) and (|2.49|) turn into linear combinations of a particularly simple kind of 
integrals we refer to as two-point integrals, 

/(£ ■ k)"^ 
(^2)a((^ _ fc)2)/3 ' a,/3GM,mEN. (2.51) 

In appendix [BJ explicit solutions for the two-point integrals (I2.5ip are given and their conver- 
gence properties are discussed. The gluon DSE (|2.47p may be rewritten in terms of two-point 
integrals as 

(^2)i+"A ^ Jl^f ' (so(aG, 1 + ac) - ^H2(l + «G, 1 + Qg)) (2.52) 

and it converges in the ultraviolet (as £ oo) for l-|-2aG > f ■ Infrared divergences are absent 
so long as ac < f • 

In the ghost DSE (I2.49p . the horizon condition allowed for a subtraction which in turn improves 
the UV convergence. The two-point integral occurring in the "unrenormalized" DSE (|2.48|) is 
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divergent, but it can be unambiguously split into a momentum-independent divergent part and 
a regular part, see appendix [B] The divergent part is subtracted using the horizon condition in 
(|2.49|1 and the formulae (IB.6P for yield the correct values for the regular part. Thus, 

= -g^NcAB^ (^So(l + az, 1 + "g) - ^^2(2 + a^, 1 + og)) (2.53) 

where the remark "reg." indicates that the regular part of Hm(a,/3) is taken. The integrals in 
(|2.53ll are then UV-convergent if ac < 1. In the infrared, poles are avoided for — | < ac < f • 
Only values for aA and ac are considered for which the integrals converge. 

Solving a well-defined (convergent) set of integral equations, we can rescale the integration 
variable in (12311 by ^ ^ Xi, and readily find that H„(a,/3) = /„(a, /3)(A;2)'^/2+™-"-/3 with 
some dimensionless functions Im{o:,(3). These power law solutions can be plugged into the 
DSEs (|2.47|) and (|2.49p and a relation between the exponents oa and ac is gained, 

aA + 2aG = . (2.54) 

Eq. (|2.54l) is referred to as the "sum rule (for power law exponents)". 

Let us emphasize that the sum rule (|2.54p can be directly traced back to the nonrenormalization 
of the ghost-gluon vertex [103| . Following Ref. |104| . a quite general ansatz for the proper 
ghost-gluon vertex is 

r,{k;q,p)=igtkj{k)qjJ2Ca{l)'^ (^)"" (^)"'' , (2.55) 

a 

where the constraint Zj + mj + nj = 0, Vz, guarantees the independence of the renormalization 
scale cj, i.e. nonrenormalization of the vertex. It is readily shown that the sum rule (|2.54p is 
not affected by a dressing of the ghost-gluon vertex such as (|2.55l) . since it turns into 

aA + 2aG = ^^^ + ^{li + mi + ni) . (2.56) 

i 

Any consequence of the sum rule is therefore understood as due to the nonrenormalization of 
the ghost-gluon vertex. 

By virtue of the sum rule (I2.54p . one of the exponents a a and ac can be eliminated. We 
choose to define 

K := ac (2.57) 

and express oa in terms of k. 

The integrals in (|2.47|) and (|2.49|) can now be written down concisely, using the formulae in 

jEi, 

(^2^d/2-i-2K ^ g'^N^AB^lA{K){k'^f/^-^-^'' (2.58a) 

(A;2)« = g^N,AB^lG{K){k^)'' (2.58b) 
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where the dimensionless functions Ia and Iq were introduced, 



1 r(f + 2k) 

2(47r)'^/2 T{d-2K)T{l + Kf ' 



4-(d-l) T{^,-K)n-K)T{\ + K) 

^""^^^ ~ (4vr)^/2+i/2 r(f-2K)r(i + f + k) • ^^-^^^^ 

Eq. (I2.58P requires that 

g^N.AB^lA = 1 = g^N.AB^Ic . (2.60) 
The infrared exponent K{d) is therefore imphcitly defined by 



(2.61) 



Above, convergence of the integrals in the DSEs was discussed. In terms of the exponent k, the 
conditions for convergence can be re-expressed. However, recall that the ultraviolet behavior 
of -Da(^) and Ddk) is not governed by the stochastic vacuum but by perturbation theory. 
The power law behavior (I2.50p is just regarded as a possible infrared limit that will actually 
be realized (see chapter [3]). Hence, the restriction on the infrared exponents by analyzing 
the UV convergence is spurious. Discarding UV criteria for both lA{k) and Icik), we find 
from the condition that no divergences appear in the infrared integration domain the relation 
— ^ < K < ^. Moreover, the horizon condition demands n > 0. Altogether, we have 

< K < ^ (2.62) 

and in this range solutions to (I2.61|) are investigated^ Manipulating the gamma functions in 
(I2.59p . one finds the compact expression 

^ ' ' sin(7rK) T{l + K + d/2)r{d/2-K) ^ ' 

for finding K{d). The solutions are shown in Fig. 12.51 These are exhaustive for k in the range 
(|2.62p and complete a recent calculation [105| . For a given Euclidean dimension d, several 
solutions for k are available. One of them can be found by algebraic manipulations, 

Kaid) = ^ - 1 , (2.64) 

as verified by plugging K,a{d) into Eq. (|2.63p . The other (irrational) solutions Kb{d) and Kc{d) 
were determined numerically. 



There are some exceptional points in the (d, k) diagram that are only solutions to Eq. (|2.63ll 
if approached in a certain direction. One of them is the point (^Z^) that solves Eq. (|2.63|) if 



^"According to the principle of locality, the Green functions should be tempered distributions which restricts 
K further. We do not focus on this aspect here. 
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2 2.5 3 3.5 4 4.5 5 5.5 

d 

Figure 2.5: Solutions for K{d). The box indicates an exceptional point, 
see text. The cross marks the solution favorable for reasons given below. 
The dotted lines are extra solutions that occur for d > 4. JVote that k is 
restricted to the range Ii2.62\) . 

approached along the hne (|2.64|) . This pecuharity appears for even dimensions d and integer 
values for k. Setting 

d = 2m-2Xe (2.65) 
K = n + e (2.66) 

with m, n E N, we approach in the limit e — > the values for even d and integer k. The 
parameter A G R controls the slope of the line on which the point (^^) is approached in 
a {d, K)-diagram. Choosing n < m guarantees the integrability of infrared poles in the ghost 
loop integral, see Eq. (|2.62p . Treating the g expressions that occur in (|2.63|1 carefully, we find 

hm U2m - 2Ae, n + e) = (-l)i+-+"(m - n)^_„ ~ ^^^^^^ (2 + A) (2.67) 

where (a)„ = a{a + 1) . . . (a + n — 1) is the Pochhammer symbol. Setting A = corresponds to 
keeping a fixed dimension throughout the calculation. From Eq. (I2.67P it can be deduced that 
for no values of m, n the equation lid = 1 (I2.63|) is solved if A = 0. Hence the failure of trying 
to find the point (^Z^) as a Landau gauge solution fTHB], with arguments in Refs. [1071 1102] 
that this solution should exist. For the set of points with n = m — 1, i.e. 




(2.68) 



we find from (|2.67p the simple result 

limIid(2m-2Ae,m-l + e) = 2 + A. (2.69) 

e— >0 
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Figure 2.6: Two solution branches of K{d) and their behavior for large 
dimensions d. The upper solution becomes complex at d ^ 20.15, the 
lower one at d ^ 26.18. 

Tuning A = —1, the points (I2.68P are solutions to (I2.63p . at least in principle. One may adopt 
the point of view that the dimension d is fixed and the points in Fig. l2.5l can only be approached 
along vertical lines. In fact, rendering d a variable is here not understood in the context of 
dimensional regularization where divergences are thus regularized. Rather, d is here a fixed 
parameter for which (|2.63|1 needs to be solved. The exceptional points (|2.68p are therefore 
discarded. 

As seen in Fig. 12.51 there are several solutions for k as a function of d. How many solutions are 
found for a given dimension? For even dimensions d, Eq. (I2.63P can be turned into a polynomial 
equation in k of degree d+1 which has the f + 1 integer solutions k G {0, 0, 1, 2, . . . , ^ — 1} for 
A given by Eq. (|2.67|) . The double root at k = is excluded by the horizon condition. Thus, 
these integer solutions lie below (or on) the branch (|2.64|) looking somewhat odd in Fig. 12.51 
There are further non-integer solutions for even dimensions, | at the most being real. 

Apart from the branch Ka{d) in Eq. ()2.64p . there are two other solutions Kb{d) and udd) in 
the range in d G (2,4], see Fig. 12.51 Let us investigate the continuity of these three solution 
branches as functions of d, cf. [75]. Here, we do not consider the limit d — > 1 since the 1 + 1 
dimensional theory does not have any degrees of freedom and power law solutions might not 
be appropriate. Instead, we discuss the high-d range. The Gribov-Zwanziger scenario should 
work particularly well for higher-dimensional YM theory where entropy is even more enhanced 
at the boundary of configuration space. First of all, the solution K,a{d) holds for all d. The 
two other solutions in Ubid) and Kc{d) in the interval 2 < d < 4 connect to the curves shown 
in Fig. l2.6l or higher d. The upper solution approximately yields = | + | for 2 < d < 8 and 
oscillates around k = | for 8 < d < 20. It ceases to be real for d « 20.15 . The lower solution 
matches Kf, = ^ — | for 2 < d < 6, in agreement with Ref. [78], but oscillates around «^ = f — ^ 
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Table 2.1: Power law solutions for propagators in the stochastic vacuum. 



for 6 < d < 26 and becomes complex at d « 26.18. The continuity of the solutions for k as 
functions of d is therefore strongest for the solution Ka{d) which holds for all d. 

A restriction on the set of solutions in Fig. l2.5l comes from dealing with only power laws for the 
propagators, due to the ansatz (|2.50p . The two-point integrals in (|2.5ip cease to be power laws 
if| + m — a — /? = 0. Actually, they cease to exist, due to a logarithmic ultraviolet divergence. 
We have argued that these ultraviolet divergences are subtracted by some renormalization 
scheme in the perturbative sector. One might therefore suggest to consider instead the integrals 
P) subtracted at ^ and then approach |+m — a — /3 = 0. For instance, if m = 0, then 

hm (Eo{a,P)-Eo{a,P)\^^ = lim r(6) [(fc^)- - (^2)-.] 

— ^-/3 ^ ^ (47r)2 r(f) ^-*o 



2 



^1-^- (2-70) 



(4vr)f r(f) 

Evidently, after subtraction of the ultraviolet divergence, the solution of the two-point integral 
(with exponent zero) is a logarithm. These functions call for a separate treatment. The allowed 
values of k (I2.62p are therefore further restricted by 

-/^-^ (2.71) 



to avoid a ghost loop that gives a logarithm. The lower solution branch Kb{d) in Fig 
oscillates around the forbidden values (|2.71ll for high d, exactly yielding (I2.7ip for d = 
6, 10, 14, 18, ... . A discussion on a possible extension to solutions with logarithms will fol 
low in section 



To summarize this section, we solved the momentum space Dyson-Schwinger equations of the 
stochastic vacuum with a tree-level ghost-gluon vertex by power law ansatze and thoroughly 
discussed the validity of the solutions. For the propagators in the temporal Coulomb and the 
Landau gauge, with d = 3 and d = 4, resp., the solutions are shown in Table EH] In 2 + 1 
dimensional YM theory, the temporal Coulomb gauge corresponds to d = 2 and the Landau 
gauge to d = 3. 
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2.7 Linear confining potential 



The Hamiltonian approach in the temporal Coulomb gauge allows for a direct calculation of 
the potential energy between two static quarks. Let the external charge density be given by 
two point charges, separated by the distance r, 

= ,5"3 (^(^(3)(x-r/2) -(5(3)(x + r/2)) . (2.72) 

The color vector Pextix) was here chosen in the 3-direction of SU{2), for convenience. From 
the Hamiltonian (|1.107|1 the |x — y| = r dependent color Coulomb potential Vc{r) can be 
identified by 

Vc{r) = ^1 d'ixyKAx) {F^'{x,y)) p^y) . (2.73) 
With Eq. (I2.72|) . the Coulomb potential Vc{r) can be written as 

Vc{r)=g^F^'{0)-F'^{r)) , (2.74) 
with the self-energy of the quarks, (^^3(0)) = {F^^{x,x)). 

An approximation is needed to express Vc{r) in terms of the known Green function Dcik). 
In arriving at the solutions for the propagators in the previous section, the ghost-gluon vertex 
was approximated to be at tree-level. An approximation for writing Vc{r) in terms of Dq that 
is equivalent to the vertex approximation would be desirable. Let us recall that factorizing the 
expectation value (^GT^G) (G) (G) implies F ^ F", see Eq. (|2.21|1 . In the same manner, 
we factorize 



F''''{x,y)^ = I Sz(G'"^{x,z){-d'^)G'^^{z,y)^ 

(g-\x, z)) i-d^) (g''\z, y)) . (2.75) 



Checking the validity of this approximation is an issue dealt with in section [3?6l The expression 
(ITTill thus turns into 

Voir) =g' I d^k k'DUk) (l - e^'^ '-) . (2.76) 
Plugging in Dcik) = B/{k'^Y^'^ from (|2.50|1 permits an explicit evaluation of Vr (r) 



Vc{r) 




sin (A; 



kr 



1 1 \ 1 1 f sinx 

+ -r-7-; 7T / dx —r—r 



2^2 



4k -1 4k 4«;(4ft; - 1) y e^'^-i 4k(4k - 1) 
-F(-4K)sin(2K7r)r^''-^ (2.77) 



^^In the first line, the pole at = is regularized by integrating only from < e ^ 1 to cxd and letting 
e — > at the end. In the second line, partial integration is used repeatedly, and in the last line an integral 



representation of the gamma function in Ref. [108) . formula 3.761.4 was used 
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for values | < k < | where this integral exists. In the latter interval, two solutions of the 
DSEs exist, see Table ITTl One of them, found on the branch Ka{d) in Fig. 12.51 yields Ha = \ 
and thus a heavy quark potential Vc{r) that rises hnearly. Prom Eq. (|2.77|1 one gets for k — 

^c"(0 = • (2.78) 

By the above linear potential, heavy quarks are held together by a constant force given by the 
Coulomb string tension ac =9^^- Thus, heavy quarks are confined. 

Let us reflect on how we obtained such a heavy quark potential. It is known that Voir) 
cannot rise stronger than linearly |109^ I110| . but it could saturate for large r (see further 
discussions in the chapters below). The linearity of the potential agrees exactly with lattice 
calculations (see e.g. Ref. [H]). Hence, support is given for the validity of the result (|2.78l) . 
It is enlightening to survey the origin of this outcome. The Gribov-Zwanziger scenario was 
exhibited in section 12.11 Gribov copies necessitate the restriction of configuration space and 
Zwanziger's entropy arguments demand the horizon condition. The simplest possible state, 
the stochastic vacuum, is deemed sufficient to describe the infrared properties of the theory. 
Dyson-Schwinger equations comprise the information on the Green functions and with a few 
approximations, we thus arrived at the linear potential (|2.78|) . One might be astonished by the 
simplicity of these ideas. In particular, discarding the entire Yang-Mills action Sym (setting 
= 1) and still producing sensible results seems puzzling. The only physical content in 
the generating functional 

that drives the (infrared) dynamics comes from the gauge fixing procedurel^ This notion goes 
by the name infrared ghost dominance. Merely the Faddeev-Popov kernel M. (|2.34p and its 
boundary conditions on the Gribov horizon dVt are sufficient to give rise to a linear confinement 
potential. The crucial infrared properties are thus governed by entropy in configuration space. 

A caveat in the above reasoning for quark confinement in the Coulomb gauge is that the 
Coulomb potential Vc{r) is only an upper bound to the gauge-invariant potential Vw{r) defined 
by the Wilson loop The Coulomb string tension is therefore larger or equal to the string 

tension from the Wilson loop, ac > o"VF- Hence, it cannot serve as an order parameter for 
the deconfinement transition. The mechanism that lowers the Coulomb string tension to its 
physical value is expected to be due to constituent gluons along the color flux tube. So far, 
the vacuum Green functions are unaware of the presence of external charges. It will be the 
subject of chapter El to investigate the effect of the back reaction of the external color charges 
onto the gauge field sector. 

Some further critical remarks are at order. 

• Perturbative sector. If the solutions in the stochastic vacuum are independent of Sym, 
how do they connect to the ultraviolet regime? 

^^Note that the Faddeev-Popov determinant is gauge invariant. 
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• Approximations. How can we estimate the effect of the approximations made in deriving 
the Green functions? 

• Uniqueness. Are the power law solutions found unique? Which one of them is realized 
in physics? 

In the following it will be attempted to give answers to the questions listed above. The latter 
item is dedicated to the next section. 

2.8 On the uniqueness of the solution 

Whether the solutions for the propagators found in section 12.61 can be regarded as unique 
solutions to the infrared sector of YM theory is an important issue addressed in this section. 
A thorough discussion of uniqueness would probably have to begin with the perturbative sector 
and any approximations should be avoided. Here, we shall turn to a much simpler task. In 
view of the approximated DSEs (|2.47|) and (|2.49|) . it is investigated if any of the solutions listed 
in Table EH] can be considered unique. 

A free parameter 

First of all, let us start with a generalization of the DSEs. In the context of renormalization 
in the perturbative regime, a local action is indispensable and the transversality condition 
on the gauge field, be it Coulomb gauge or Landau gauge, needs to be taken "off-shell". A 
crucial difference to the "on-shell" formulation is that the tree-level Green functions are slightly 
different. In particular, the ghost-gluon vertex does not come with a transversal projector 
attached to its gluon leg. Instead of T^{k,q), the local formulation yields ^1{q) with 

rliq)=igqk ^ Tl{k,q)=igtkj{k)qj . (2.80) 

Consequently, in the local formulation not all Feynman graphs are transverse in the gluon 
momenta. For instance, the 0{g'^) gluon propagator in Landau gauge is known to be comprised 
by a gluon loop, a ghost loop and a tadpole term. Only the sum of these diagrams yields a 
transverse gluon propagator. This subtlety slightly alters the DSEs derived in section 12.61 We 
do not repeat the derivation here, but simply state the results, 

D^\k) = -N, j dHTUt)T%i-k)DG{(i)DG{t-k) (2.81a) 

D^^{k) = k"^ -g'^Nck'^ jd'^l (l-{i-kf^DAii)DG{i-k) , (2.81b) 

as opposed to Eqs. (|2.45l) and (12.48(1 . Here, the proper ghost-gluon vertices were again replaced 
by the tree-level ones. The ghost DSE stays the same. In the gluon DSE (I2.81a|) . on the other 
hand, note the tiny difference that the ghost-gluon vertices ^^il) ones defined in 
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Eq. (|2.80p . Although the nonperturbative ghost loop on the r.h.s. of Eq. (|2.81ap should be 
transverse with a proper ghost-gluon vertex^fl it is generally not transverse with the tree-level 
ghost-gluon vertices (|2.8Qp . Transversality is lost due to the approximation of the vertex. To 
control the longitudinal contributions, contract Eq. (|2.81ap with the tensor 

Rl-ik) := - Q^i,{k) . (2.82) 

The factor (iij{k) would give no contribution if both left and right hand side of (|2.81al) were 
transverse. The sensitivity of the solution with respect to C thus measures the violation of 
transversality and the quality of the approximation used. 

Before we use the tensor R^j, let us calculate the inverse gluon propagator in Eq. (|2.81ap 
explicitly with the ansatz (|2.50p for the ghost propagator £>rj(fc)P^I 



B-^Dr^{k) = g^N.Jd^e 



T{d/2 - K)T{d/2 - 1 - K)r{2K - d/2 + 2) 



2 ^ 1 .,2^d/2-l-2. / Tid/2-K)^T{2K + l-d/2) 

^ '{47ry/^^^ 1 2r{l + Kmd-2K) 



r(i + K)2r(d-i-2K) 



(k) 



r{d/2 + 1 - K)r{d/2 - 1 - k)T(2k - d/2 + 2) 

+ r(i + Ac)2r(d-2K) ' 

= lA[Sij- {4k -d + 2)£,,{k)]{ky/^~-^-^\ (2.83) 

with the function I a given by Eq. (|2.59ap . The Lorentz structure of D^j^{k) is most easily 
recognized by the last line in (I2.83p . It reproduces a familiar result from perturbation theory 
(k = 0): the ghost loop has longitudinal components. Furthermore, one can see that quite 
generally the ghost loop is not transverse. Transversality is maintained if and only if the 
infrared exponent k satisfies 

d-1 , , 

AC = ^ , (2.84) 

as seen directly in Eq. (I2.83p . For other values of k, projections of the gluon DSE (I2.81al) with 
the tensor R\j{k) in Eq. (I2.82P will depend on the value of Q. In this case, the computation 
of the DSEs can be repeated to give a condition /^^(k, C) = 1, analogously to Eq. (|2.6ip . The 
latter condition yields a solution K{d, Q) that is not unique in the sense that it depends on C,. 
This dependence is due to the approximation of the ghost-gluon vertex and must be spurious. 
In principle, C can be regarded as a free parameter. A good approximation will give a weak 
C-dependence. 

We now focus on d = 3 where one can derive that [113] 

I.j^ C) = 32(k-1)accos(2^k) ^ 

l^d{^^,^) (i + 2K)(3 + 2K)(2(l-K)-C(l-2K))sin(7rK)2 



^^If the ghost loop dominates in the infrared, it must itself be transverse. 

^■*The tensor integral may be expanded into scalar integrals by means of the Passarino-Veltman algorithm 

fml. 
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Figure 2.7: The function K{d = 3, C)- -Due to the approxinaation of the 
ghost-gluon vertex, the infrared exponent k depends on the parameter 
C, except for Ha = 1/2. 

determines the value of spuriously dependent on C,. For C, = 1, the two solutions Ha = \ 
and Kb = 0.398 from Table ET] are recovered. In Fig. 12.71 it is shown that other solutions for 
K can be found if C / 1- The important feature to notice here is that while n^, depends on C, 
Ka is C-independent. This can be seen immediately from (|2.84l) by noting that Ka = \ yields 
a transverse ghost loop. Although favored in Refs. |113l 1105] for 3-dimensional Landau gauge, 
the solution k^, is quite sensitive to approximations in the off-shell formulations, and hence 
might be unphysical. For the solution Ka = \, the approximations made are more trustworthy 
and it allows any value of C- Further arguments in favor of Ka will be given in section 13.51 
in the context of well-defined infrared limits of the ghost-gluon vertex. Here, let us finish the 
discussion by noting that it is Ka that naturally yields a linearly confining potential in the 
Coulomb gauge. 

In the Landau gauge (d = 4), the solutions listed in Table ETT] are all sensitive to the value of C, 
and none of the solutions is characteristic in this respect. The accepted value for the infrared 
exponent, is k = 0.595, as reviewed in Ref. [54] . 

Logarithmic degeneracy 

The infrared power law ansatze made so far yielded integrals that can be generically represented 
by 

I(a; p2) - j dH {l^T^{l -p)= I{a){p^f+^ , p^O (2.86) 

for some suitable function (/>(^^). The solutions to the integral in Eq. (|2.86p with 0(^^) being 
a power law were discussed, and it is clear how to find the exponent ip that depends on the 
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dimension d of the integral and the function (p. The dimensionless numbers I{a) can be 
computed with the methods given in Ref. [103| . Each of the DSEs are of the form 

ipy = Iia;p^). (2.87) 

Matching both sides of Eq. (I2.87p . we get the sum rule P = a + ip on the one hand and the 
condition 1(a) = 1 on the other hand. This determines a as well as /3. 

In order to incorporate logarithms into the infrared ansatze such as (£^)° ln"'(^^) we note that 
simple powers can serve as generating functions of the logarithm to some power n G Z, that 
is, 

an 

-(p2)" n>0 



(/)"ln"/ 



5a" ' 

da(p2)" n<0 . (2-88) 



|n| 

As long as the generators of the logarithm, i.e. the derivatives or the a-integrals, may be 
interchanged with the loop integration of d'^l we can compute 

J(a, n; p^) — [ d'^£ (f)'^ In" f(f){i - p) , p^O (2.89) 



by applying these generators to Eq. (I2.86|) . E.g., for n = 1 we simply have to write 

d_ 

da 



J(a,l;p2) = I dH{l^T\nf(t){l-p)= I dH-^{l^T(t){l-p) 



— /(a;p - 



+ 1(a) In/^ (p2)"+^ 



da \ da 

\np^ I{a;p'^) (2.90) 

In the limit p — > 0, the logarithm dominates any constant and thus only the above term 
remains. For any n G N, the result yields 



J{a,n;p') = (p^)"+^ ^ (^) In-'^P^ (2.91) 



f n \ d^I{a) 

k=0 

Equivalently to the case n = 1, one gets the infrared limit 



lim J(a,n;p2) =ln"//(a;/) ,nGN. (2.92) 
p— >o 

Let us turn to n < 0. From the prescription (I2.88P one can derive the formula for J by repeated 
partial integration. To prove the result thus obtained, it is easier to turn the integral form of 
Eq. (|2.88|) into a differential equation, 

()\"-\ 

Ha;P^) = jri;jJia,n;p^) ,n<0. (2.93) 
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For n = — 1, we get 

with a constant C{p'^) that does not depend on a. 
Proof: 

= (p2)"+^/(a) = /(a;p2). (2.95) 

Note in Eq. (|2.94p that for p — > only the term with /c = is relevant. The arbitrary integration 
constants C(p^) can always be chosen such that they are subleading in the infrared. The result 
(|2.94p can now be generalized to any n < 0. By induction, one can show that 

J{a,n;p^) = ^, T] Y , + 0(a'"-^') ,n<0 (2.96) 

lnHp2li^^^ a«El"'fc. Ini:l"'^»p2 

satisfies the differential equation (|2.93|1 . The infrared limit yields 

lim J(a,n;p2) = ^i^, n € Z" . (2.97) 

In summary, we have shown that in the infrared limit for any n G Z the logarithms in the 
integrals J(a,n;p^) in Eq. (|2.89p can be removed from under the integral and placed in front 
of the integral, evaluated at the external momentum p^, i.e. 

lim J(Q,n;/) =ln"p2 J(Q,.^2^) ^ nGZ. (2.98) 

Therefore, if an integral equation of the DSE type, as given in Eq. (|2.87|) . is solved by plain 
power laws in the infrared, an ansatz of the kind (£^)"ln"(^^) with n G Z will also be a 
solution. To see that, note that the logarithm will appear on both sides of Eq. (I2.87P and the 
determining equation /(a) = 1 for the infrared exponent a is left unchanged. In this sense, 
there is a degeneracy in the infrared power law solutions that are insensitive to logarithms. 

However, it is not clear how to rigorously generalize these results to n € M. There are methods 
to define fractional derivatives and fractional integration [114] and a careful treatment might 
recover the result (|2.98p for arbitrary n E M. This is something that would be worth looking 
into. We conjecture here that in the infrared limit. 
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Figure 2.8: Left; Logarithmic dressings of a power law with a = 1 in the 
infrared. The curves confirm that in the infrared Emit, the logarithm 
simply factors out, according to Eq. Ii2.99\) . Right: Here, a = 0.75. Tie 
upper curve with m = 0.5 yields In^ k"^, the middle curve with m = 
yields Ink^ and the lower curve with m = —0.5 gives Inln/c^ in the 
infrared. 



where Ho(a,/3) are the functions defined in (I2.5ip . basically the l.h.s. of (|2.99p without the 
logarithms. Hence, Eq. (|2.99l) states that in the infrared limit the logarithms factor out. 
For the time being, it shall suffice to corroborate this conjecture by a few modest numerical 
investigations. In Fig. 12.81 we have implemented the integrals 



J{k) :-- 
with the simple functions 



(2.100) 




for A; < 1 
for k > 1 



(2.101) 



for several values of a and m. In the left panel of Fig. 12.81 the infrared limit indicated in Eq. 
(fZTOOll with I{a) = r2(3/2-a)r(2a-3/2)/ (87r3/2r2(a)r(3 - 2a)) following from Eq. (iRBajl 
is found exactly. This supports the hypothesis (I2.99p . 



In Ref. |115| . the possibility of logarithmic corrections to the power law solutions was discussed 
using the angular approximation. It is easily verified that the angular approximation is unable 
to reproduce reliable results for the class of two-point integrals we are dealing with. Moreover, 



the hypothesis (|2.99|) needs a more refined treatment than the angular approximation can 
provide. The denial of logarithmic corrections in Ref. [115| is therefore likely to be due to the 
angular approximation. 

A further point is made in the right panel of Fig. 12.81 If the power law exponent a in the 
integrand of Eq. (|2.100p is such that the integral is dimensionless, i.e. a = | with m = 0, 
the results will be a logarithm, rather than k^ = const. If furthermore m 7^ 0, the result 
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is In^"''^'" /c^ . This can be seen in the two upper curves of the plot. The lower curve with 
the choice m = — | renders the exponent of the logarithm zero. The asymptotic result found 
numerically is the function InlnA:^. With zero exponents of power laws incrementing the 
logarithm power, and with zero logarithm power incrementing the Inln/c^ power, it is intuitive 
to expect this logic to continue to In In In /c^. In In In In A;^ and so on. Let us have these general 
ideas be followed up with a specific example. Consider an iteration of the DSEs in d = 3 
with K = i as a starting point, i.e. d'^q {k) ~ The ghost loop then produces, 

after UV subtraction, a logarithm in the IR, cf. Eq. (|2.70l) . Consequently, the IR gluon 
propagator yields D^^ {k) ^ l/{h\k'^). Plugged into the ghost DSE, the log sum rule (|2.99|) 
yields that D^q {k) ~ (In /i;^)/(A;^)^/^. Every iteration increments the logarithm power of the 
ghost propagator, so after the nth iteration, Dg^(A;) ~ {Yn^ k"^) / {k"^)^/^ . The choice = \ 
must therefore be excluded, see Eq. (I2.7ip . 

However, as long as such subtleties are avoided, the solutions for the ghost and gluon propa- 
gators can be appended by 

DA{k) ^ 



{k^y+'^A InT (A;V/i2 
B 

with the conditions (I2.54p and (|2.6ip from the power law analysis and the log sum rule 



Dcik) = (2-102) 



7 + 2(5 = (2.103) 

from Eq. (|2.99|) . Due to the factoring out of the logarithms in the infrared limit, no changes are 
expected for the value of At = ac- Neither are any qualitative changes expected if logarithms 
are present in the infrared, for they are always subdominant to power laws. 

To summarize this section, the power solutions found in sections [23] are not necessarily unique. 
First of all, a free parameter can be dialled in the off-shell formulation to yield different values 
of the infrared exponents. Only the one solution, = ^ for d = 3 (see Fig. 12.5(1 . is insensitive 
to this parameter. Secondly, the amendment of logarithms to the power law solutions give one 
possibility to generalize the solutions. Logarithms will be of particular interest in the discussion 
of the ultraviolet behavior of the Green functions. That will be the topic of chapter ID 
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3 Variational solution for the vacuum state 

While in the previous chapter the stochastic vacuum state was investigated to derive the in- 
frared properties of the theory, this chapter is devoted to determine the vacuum wave functional 
in the temporal Coulomb gauge by variational methods. A Gaussian type of ansatz wave func- 
tional is put forward to calculate the energy and minimize it with respect to the parameters 
of the ansatz. The solutions for the Green functions thus obtained turn out to approach the 
power law solutions for d = 3 of the previous chapter in the infrared. In the ultraviolet, the 
propagators are those of free particles, up to logarithmic corrections. 

By virtue of the (time-independent) Yang-Mills Schrodinger equation, additional information 
is available from the spectral properties of the Hamiltonian, as opposed to the exclusive consid- 
eration of Dyson-Schwinger equations. In this respect, the Hamiltonian approach is superior 
to the Lagrangian approach, utilizing the principle of minimal energy. In quantum field theory, 
variational methods have enjoyed considerable attention |116l 1117] . However, critical remarks 
were issued by Feynman. In the 1987 Wangerooge conference, Feynman claimed that varia- 
tional methods were "no damn good at all" in quantum field theory |117| . essentially for the 
following reason. The only wave functionals that can be used reasonably, are Gaussians. Any 
corrections to the Gaussians can only be calculated numerically, involving intolerable errors. 
The inaccuracies in the ultraviolet sector spoil the reliability of infrared quantities, especially 
in non-linear theories such as QCD where the UV and IR modes mix. Therefore, if Gaussian 
wave functionals are not good enough, there is no chance in finding reasonable results with the 
variational method, according to Feynman. The calculations presented in this chapter are, of 
course, overshadowed by the latter remarks. Nevertheless, on a qualitative level, the results 
turn out immune to the criticism. First of all, a generalization of Gaussian wave functionals 
is feasible and will be exhibited in the following section. Secondly, the Gribov-Zwanziger sce- 
nario infers that the gauge-fixed configuration space itself governs the infrared, regardless of 
the ultraviolet accuracy of the wave functional. The crucial infrared properties do not follow 
from the variational minimization of the vacuum energy, but can be extracted directly from the 
path integration inside the Gribov horizon, along with the horizon condition. We are therefore 
confident that the infrared power laws do not suffer from any other approximations than those 
discussed in the previous chapter. 

It is well-known from nuclear physics that the ground state energy is quite insensitive to 
errors in the wave function. The accuracy of the wave function is tested by the calculation 
of transition amplitudes, not by the vacuum state itself. In the same sense, one may expect 
that the Yang-Mills vacuum energy density is not too far off using the crude approximation 
of a Gaussian vacuum wave functional. Other quantities, such as the vertices, may be more 
sensitive to changes in the wave functional. These assertions will be discussed in the course of 
this chapter. 
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3.1 Gaussian types of wave functionals 



Gaussian types of Coulomb gauge wave functionals are put forward in this section. They were 
first proposed in Refs. [85l [83] . 



\ I d\xy]A'i{xMx,y)A'i{y) 



(3.1) 



where u;{x,y) is a variational kernel and M normalizes the wave functional, (^'| = 1. [Note 
that henceforth, we refrain from using bold-faced symbols for 3-vectors, unless necessary for 
clarity.] The factor of the Faddeev-Popov determinant with the real exponent A allows for a 
further enhancement (if A > 0) of the probability density near the Gribov horizon. One specific 
choice, A = ^, resembles the usage of a "radial wave function", familiar from the calculation of 
the hydrogen atom. In this case, the evaluation of expectation values becomes straightforward 
because the generating functional of full Green functions. 



= (^'lexp 



d'xf,{x)AUx) 



j VAJ^-^\A]e^Y> - j d^[xy\Al{x)u{x,y)Al{y) + j d^xf,{x)At{x) 



(3.2) 



is actually a Gaussian path integral for A = ^. Before we write down the result, let us try to 
understand the meaning of the parameter A. 

A Dyson-Schwinger equation (DSE) for the gluon propagator can be written down for the 
path integral in Eq. (13. 2p in order to relate uj and A to the gluon propagator Dij, 







M"^ / VA 



-J 



1-2A 



6Af{x) 

d'[x'y']Al{x'Mx',y')Al{y')+ I d'x'f,ix')AUx') 



exp 



Af' I VAj'-^^i{l-2X)j^^-2 I dhu{x,z)At{z)+Um{x)jm{x) 



exp 



d'[x'y']Al{x'Mx',y')AUy')+ I a(Vj^(x')^^(x') 



^3^' ^Cf^l\ AC f^l^ 



(3.3) 



Taking a derivative w.r.t. jj{y) and setting sources to zero yields 

= (l-2A)(^|ij|y^5(y)^^-2| d'zco{x,z){At{z)A';{y))^ + tt^{x,y) . 



(3.4) 



Setting A = ^ in the above equation, the gluon propagator is seen to be determined by the 
variational kernel u;(x,y). 



Af(x)4(y))^ = i5%(x)a;-i(x,y) 



(3.5) 
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For A / ^, we rewrite the first term in Eq. (|3.4p by partial integration, 



A 



exp 

We now introduce an approximation for the evaluation of Green functions, referred to as the 
loo'p expansion. Whereas its common definition concerns an expansion in h, we here use the 
term loop expansion for ordering diagrams by the number of loops, integrating nonperturbative 
ghost or gluon propagators. Applied to Eq. (13. 6p . it can be realized that the first term comprises 
one loop, and the second term comprises two loops0 In the one-loop expansion, we may write 
Eq. (1331 as 

d'zu;ix,z)[AUz)A'^{y))^ + t^^{x,y) . (3.7) 
We recognize in Eq. (13. 7p the ghost loop, here denoted by Xij with 

Since by the quantity Xij the non-trivial metric of the gauge- fixed variables is expressed, the 
momentum space quantity 

1 

id - 1) 



Xik) = i)(V2_i) %(^)^°' J '^'^ y) e-^'-^''-'^ (3.9) 



is also referred to as the curvature ^83], with the spatial dimension d = 3. Introducing for the 
gluon propagator the function 0(x,y) by 

(yl^(x)A5(y))^ = h,,{x)n-\x,y) , (3.10) 

and employing the one-loop approximation, Eq. (|3.7|) can be found to yield 

n{x, y) = u{x, y) + (1 - 2\)x{x, y) . (3.11) 

Recall from the last chapter that the DSEs, following directly from the boundary conditions 
on the path integration within the Gribov horizon, fully determine the Green functions of the 



^Each functional trace "Tr", as found in In = TrlnG ^, gives rise to a loop. 
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theory. In the variational approach, on the other hand, the Green functions have a parametric 
dependence on the kernel oj and on the parameter A. The DSE (|3.1ip merely tells us how the 
gluon propagator is related to the variational kernel w. The equation that actually determines 
uj is the Schrodinger equation, H \^) = E \^) , that yields for the vacuum state a variational 
equation of the Rayleigh-Ritz type, 



mm 



(3.12) 



The Schrodinger equation is not strictly a DSE, regarding its origin. Nevertheless, we will 
often refer to the Schrodinger equation as a DSE since, after all, it has the same form of a 
non-linear integral equation. 



One infers from Eq. (|3.11|) that for the evaluation of the gluon propagator to one-loop order, 
the Faddeev-Popov determinant in (13. 2|) can be replaced by a Gaussian, 



J 



exp 



d'[xy]A'^{x)^^{x,y)A]{y) 



(3.13) 



and it suffices to use the following generating functional 



Zx[j\ 



I VAew - I d'[xy]Al{xMx,y)At{y) + J <fxjl{x)Al{x) 



exp 



d\xy]j^{x)t^n{x)^{x,y)j^{y) 



Prom the relation 



{0[A])^ = O 



5_ 



Zx[j] 



(3.14) 



(3.15) 



one can verify that the gluon propagator yields Eq. (|3.1ip . 

It was shown in Ref. [85] that for the calculation of expectation values in the Schrodinger 
equation to one-loop order, the Faddeev-Popov determinant can be replaced according to Eq. 
(I3.13p . Thus, the "Gaussian type" of wave functionals (13.11) are indeed Gaussian to the level 
of approximation. The functional form of Wick's theorem [118] thus becomes applicable and 
results in the concise prescription (|3.15ll to evaluate expectation values. With u! replaced by 
in Eq. (I3.14p . the Schrodinger equation minimizes the energy with respect to the gluon 
propagator ()3.10|1 and the value of A is therefore irrelevant [HS| • We will make the choice A = ^ 
to solve the Schrodinger equation following Ref. [83]. Let us point out here that the restriction 
to the Gribov region is abandoned in the evaluation of the Gaussian path integrals, see Eq. 
(|3.14p . It is therefore of particular interest to compare the results to lattice calculations. 

With a Gaussian wave functional at hand, the gauge fields can be transformed into the particle 
representation and interpreted as particles with energy modes uj{k). Since the function u;{k) 
turns out be a non-trivial dispersion relation, gluons are interpreted as quasi-particles. This 
will be further discussed in chapter [5l 
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3.2 Minimizing the energy density 

The vacuum energy density of Yang— Mills theory with vanishing external charges, p'^^i (x) = 0, 
was calculated in the state (|3.1I) with A = ^ to two-loop order (in the energy) in Ref. [83j . 
Subsequently, the Yang-Mills Schrodinger equation was used to determine the vacuum state 
by minimizing the energy density w.r.t. the kernel oj{k). This leads to the equatior^ 

£(i?) = (3.16) 

which is basically the only equation that follows from minimizing the energy. Since it gives rise 
to a dispersion relation with a mass gap, we refer to Eq. (|3.16|) as the gap equation. The gap 
equation has quite a complex structure, involving the ghost propagator in particular, and we 
have to resort to auxiliary equations in order to calculate the expectation values. One of these 
equations is the ghost DSE derived in chapter [2] and another one accounts for the factorization 
of the Coulomb potential, see below. The entire calculation of Ref. [83j shall not be repeated 
here, but only the essential steps are outlined and some subtleties are pointed out. 

For the expectation value E = {H[A,Il\) with the Yang-Mills Hamiltonian //[^d, 11] given 
by Eq. (|1.107|) . to be calculable by the prescription (|3.15|1 . the momentum operators 11 are 
eliminated by action on the wave functional ^[A\ = {'^\A). We introduce by 

^UU.,^(^nAi) - (-4^-li||^)^e-i/- 

= iQl{x)(^nA]) (3.17) 



the quantity 

Ql{x) = I d'yuj{x,y)AUy) - ^tfc,(x)Tr (^"/(x)) . (3.18) 
One can then show that for the given wave functional ^'i [A], 

2 

E={H[A,U]) = {H[AQ[A]]) (3.19) 

and Wick's theorem (I3.15P becomes applicable. The vacuum energy E = + Ep + Ec can 
be computed [83] and it yields the kinetic energy 



Ek = ^^{2ttY5\Q) / d'k'-^^-^^^ , (3.20) 



the magnetic potential 



, M|^,^,2,)3,3(„) j AA'-^^ (3 - ,i . tf-) , ,3.21, 



^In this section, all expectation values are taken in the vacuum state l|3.ip with A = |. 
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and the Coulomb potential 



{k - k'Y 



{[u:{k)-xm-[^{k')-x{k ')]) 
oj{k)oj{k') 



2 



(3.22) 



These terms arise from the kinetic, magnetic, and Coulomb parts of the Hamiltonian, listed in 
that order in Eq. (|1.107l) . The overall volume factor of (27r)^(5^(0) is due to the translational 
invariance in the absence of localized external charge distributions. In the above expressions, 
two form factors were introduced. The ghost form factor d{k) measures the deviation of the 
ghost propagator Dcik) from tree-levelp 

Dcik) = ^ , (3.23) 

and the Coulomb form factor f{k) measures the deviation from the factorization in the expec- 
tation value of the Coulomb operator F (cf. Eq. (|2.75p ). 

(n^) ^ . (3.24) 

From the definition (I3.9|l . the curvature x{k) is found to define the following integral, 

The gap equation can be derived from Eq. (I3.16P using 

to yield in momentum space |83| 

oj\k) = k^ + x\k) + Uk) + l!^ (3.27) 



with the abbreviations 



I. = jdH[?,-{k-lY)-^^ (3.28a) 



[u{l)-x{i)+x{k)?-u:{kf 
In addition, the ghost from factor d{k) obeys the ghost DSE (I2.48p . 



(3.28b) 



' -1-9'^ fdH(l-{i.k)A (3-29) 



d{k) ^ 2 J \ ' ' J {i- kYto{(.) 



'Note that this definition is different from Ref. |83) . Here, at tree-level d{k) = 1. 
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and the Coulomb form factor obeys another integral equation that will be discussed further in 
section 13.61 

The equations (|3.27|) , (|3.29l) and (I3.30p are here collectively called the Dyson-Schwinger equa- 
tions (DSEs) of our approach to Coulomb gauge Yang-Mills theory. The approximations used 
in deriving the DSEs include the tree-level approximation for the ghost-gluon vertex and the 
one-loop expansion in the gap equation. 



3.3 Renormalization 

Each one of the Dyson-Schwinger equations, the one for the ghost form factor d{k), the 
Coulomb form factor DSE for f{k), the gap equation for uj{k) and also the curvature integral 
xik) require subtraction of the UV divergences. This can be seen by plugging in the ultraviolet 
behavior of the form factors known from perturbation theory. As discussed in more detail in 
chapter HI the ultraviolet asymptotic propagators should attain tree-level up to logarithmic 
corrections. The tree-level values for the ghost form factor as well as the Coulomb form factor 
are simply unity. Whereas the latter are genuinely instantaneous in our approach, the gluon 
propagator may be regarded as the equal-time part of the tree-level propagator in Euclidean 
spacetime, 

oo 

DUk,At = 0)= /^^o|^ = ^ (3.31) 

— oo 

and a first gues^ for uj{k) in the UV therefore is uj(k oo) = k, yielding the dispersion 
relation of a free massless particle. By simple power counting in the ultraviolet, the integrals 
in (|3.29|) and (|3.3Q|1 are found to be logarithmically divergent. The gap equation (|3.27p and 
the curvature integral (I3.25p . on the other hand, contain power divergences. 

Renormalization, the concept of redefining the parameters of a theory such that expectation 
values are free of infinities, can be introduced multiplicatively in perturbation theory [119| . 
That is, the local field operators are multiplied by renormalization constants such as to remove 
the divergences occurring in the perturbative expansion. In the Lagrangian language, multi- 
plicative renormalization may be described by appropriate counter terms in the Lagrangian 
that respect its symmetries. The arbitrariness in defining the finite part of a divergent quantity 
can be controlled by the renormalization group. Nonperturbatively, a systematic concept of 
removing divergences is lacking. This also applies to the present nonperturbative approach to 
Coulomb gauge YM theory. In Ref. |120| . counter terms to the Hamiltonian were introduced 
which are seen — a posteriori — to remove the divergences in the DSEs. Below, we will follow 

■^Recall that Dij{k) = ^Uj{k)uj-^{k). 
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the renormalization procedure given in Ref. jT20] . Thus, some additional parameters are in- 
troduced into the theory. It will be seen that, at least for the asymptotic IR and UV behavior, 
the solutions are independent of these parameters. 



Renormalization of the Faddeev— Popov determinant 



The curvature integral (I3.25P is linearly divergent. Prom the identity J = exp(— j AxA)^ see 
Eq. (|3.13l) . valid to one-loop order, it is evident that the divergences in x{k) can be traced 
back to the Faddeev-Popov determinant. Eq. (|3.13|) suggests that the counter terms required 
to renormalize the Faddeev-Popov determinant, or more precisely its logarithm, have to be of 
the form ~ / AA . Since the log of the Faddeev-Popov determinant can be regarded as part 
of the "action", we renormalize the Faddeev-Popov determinant as 



J ^ J ■ A J = exp 



TrlnG-^ + C^(A) / (rxAl{x)Al{x 



or by using the representation Eq. (|3.13l) . we obtain 



J ■ A J = exp 



A{x-C^{K))A 



(3.32) 



(3.33) 



Obviously, the counter term C-^(A) has to be chosen to eliminate the ultraviolet divergent part 
of the curvature x{k)- Thus the renormalization condition reads in momentum space 



X{k) - C^(A) = finite . 



(3.34) 



As usual, there is some freedom in choosing the finite constants of the right hand side of 
Eq. (I3.34p . In principle, we could just eliminate the ultraviolet divergent part of xi^) by 
appropriately choosing the counter term C^{K). However, it is more convenient to choose 
C-^{K) to be the curvature at some renormalization scale /i, resulting in the renormalization 
condition 



Cx(A) = x(/i) 
and in the finite renormalized curvature 

X{k) = x{k) - x(/^) • 



(3.35) 



(3.36) 



It is easy to check that this quantity is indeed ultraviolet finite and obviously it satisfies the 
condition 



fi) = . (3.37) 
By adopting the renormalization condition Eq. (I3.35p . the renormalized quantity x{k) in Eq. 



(I3.36P depends on the so far arbitrary scale /x. By choosing the renormalization condition 
Eq. (|3.37|1 this renormalization scale becomes a parameter of our "model", since it defines the 
infrared content of the curvature x{k) kept in the renormalization process. For instance. 
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3. Variational solution for the vacuum state 



choosing /x = 0, the whole infrared divergent part of the curvature is chopped off. The 
numerical calculations below use a finite fi > 0. 

In Ref. [12 Ij . a different renormalization condition was chosen, keeping from the ultraviolet 
divergent quantity x(//) the finite part x'(/^)- This amounts to setting 

C^iA) = x(/i) - Xip) , (3.38) 
which results in the renormalized curvature 

X{k) = m+X'{f^) ■ (3.39) 

Only the divergent part of the curvature is subtracted, keeping fully its finite part. Here, we 
have an extra parameter x'(m) of the theory and the renormalization scale n is not related 
to the zero of the renormalized quantity x(^) in Eq. ()3.39|1 . The latter method of removing 
divergences may be looked upon as an alternative, the former (I3.36P will be used below. Both 
methods can be accounted for by the specific choice of the counter term AJ' to the Faddeev- 
Popov determinant. 



Counter terms to the Hamiltonian 



In order to eliminate the divergences in the gap equation, it turns out the counter terms 

AH = Co{A) j d^x Al{x)Al{x) + iCi{A) j d^ x Al{x)Ill{x) (3.40) 

should be added to the Hamiltonian. Here the coefficients Co(A) and Ci(A) depend on the 
momentum cutoff A and have to be adjusted so that the UV singularities in the gap equation 
disappear. Note that these coefficients multiply ultralocal operators which are singular in 
quantum field theory. For the wave functional at hand, the expectation value of the counter 
term Eq. (|3.40|1 is given by 

AE = Co{A)^tii6'"' j d^xuj-\x,x)-Ci{A) J d^x (^^(x)Q^(x)) . 



(3.41) 



With the identity 

{Atix)Qt{x)) 



d^x'oj-\x,x') {uj{x',x) - x{x',x)) 5° 



shown in Ref. [83], Eq. (|3.4ip can be written in momentum space as 



AE = {Nl - 1) (27r)3r(0) 



Co (A) / d^k 



uj{k) 



Ci / 



uj{k) - x{k) 



io{k) 



Taking the variation of this expression with respect to uj{k), we obtain 



5AE 

Suj{k) 



-Co(A)- 



Ci(A) 



(3.42) 



(3.43) 



(3.44) 



'uj^k) 'iv-^{k)_ 

Note that x{k) depends only on the ghost propagator but not on the gluon energy w, at least 
as long as x{k) is not yet the self-consistent solution. Adding the counter terms from Eq. 
(|3.43|) to the energy in the gap equation (I3.26p . one finds 

u:\k) = k^ + x\k) + /o + Uk) - 2Co(A) - 2Ci{A)x{k) . (3.45) 
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Renormalization of the gap equation 

We now turn to the renormalization of the gap equation (|3.45|) which already includes the 
counter terms. After the renormalization of the Faddeev-Popov determinant, the curvature 
x{k) can be replaced by the renormalized one Eq. (|3.36p . Note that the Coulomb integral 
Iuj{k) does not depend on any constant part of the curvature, see Eq. (I3.28bp . so that x{k) 
could have been replaced right away by the finite quantity Replacing x{k) by x{k) and 

using the relation 

Uk)=lL^Kk) + 2miL'Kk) , (3.46) 
the gap equation in Eq. (|3.45l) becomes 

u^\k) - x\k) =k' + ll + li^\k) - 2Co(A) + 2xik) [l^'Hk) - Ci(A)) . (3.47) 

The integrals li"'~^'^^(/c) are the linearly (n = 1) and quadratically {n = 2) UV divergent parts 
of Iuj{k) in Eq. (I3.28bp . see also Ref. [83]. The integral in Eq. (|3.28ap is also quadratically 
divergent but independent of the external momentum. We can therefore eliminate all UV- 
divergences by choosing Co(A) ~ A^ and Ci(A) ~ A. In principle, the coefficients of the 
counter terms Co(A) and Ci(A) have to be chosen to eliminate the UV divergent parts of the 
quantities appearing in the gap equation. This means that we should choose these infinite 
coefficients as 

(l'. + lL'\k)) -2Co(A) = (3.48) 

\ / UV divergent part 

lL'\k) -Ci(A) = (3.49) 

UV-divergent part 

Note that the UV divergent parts of l!^~^''^\k) are by dimensional arguments independent of 
the external momentum k. Technically it is more convenient to choose the following alternative 
renormalization conditions. As usual we have the freedom in choosing the renormalization con- 
ditions up to finite constants. Given the fact that is independent of the external momentum 
and the differences 

A/W(A:,z.)=/W(A;)_jH(^) (3.50) 
are UV-finite, we can eliminate all UV-divergences by choosing the renormalization conditions 



lO + 42)(A: = z.)-2Co(A) = (3.51a) 
li'\k = u)-C{A) = 0, (3.51b) 

where v is an arbitrary renormalization scale, which could be chosen to be the same scale 
^ of the renormalization of the curvature, but given the fact that with the renormalization 
prescription (I3.35p . the scale fj. becomes a physical parameter, the two renormalization param- 
eters u and ^ need not necessarily be the same. With the renormalization conditions (|3.5ip 
the renormalized (finite!) gap equation reads 

u;\k)-x\k) = + A/(2)(fc,z.) + 2x(A:)A/W(A:,z.) . (3.52) 
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Assuming that the integrals ' (A;) are infrared finite and furthermore that the renormal- 
ized curvature xik) is infrared divergent, the infrared hmit of the renormahzed gap equation 
is given by 

lim (uik) - m) = ^iL'Hk = 0,u). (3.53) 
fc— >o 

The 't Hooft loop, which yields a perimeter law in the confinement phase, can be explicitly 
calculated using its representation found in Ref. |122| . In the Gaussian vacuum, it is sensitive 
to the quantity (|3.53p and it is shown in Refs. [1231 1121] that only the choice u = gives the 
perimeter law. With this choice the renormalized gap equation becomes 

u;\k)-x\k) = e + Al!,^\k,0) + 2xik)^lL'Hk,0) . (3.54) 



Renormalization of the ghost and Coulomb form factor 



The Dyson-Schwinger equation for the ghost propagator is given by Eq. (|3.29p and bears a 
logarithmic divergence in the ultraviolet part of the integral. By a simple subtraction at /x^, 
this divergence is removed. Hence, 



1 



1 



d{k) d{fid) 



9 



dH[i- {i ■ kf 



d{l - k) 
{£ - kYoj{i) 



{k ^ fid) 



(3.55) 



is a finite equation. In principle, one can choose different from /i. In view of the horizon 
condition, it is convenient to set fid = while keeping a finite fi (in particular for the sake of 
the curvature). 



The DSE for the Coulomb form factor is given by Eq. (|3.30|) and can be subtracted at /xj to 
yield the finite expression 



f{k) = fifif) + 



d^£{i-{i- k) 



2^ d\l-k)f{l-k) 

{£ - kyu{e) 



{k ^ fif) 



(3.56) 



We choose here fxj = fx. The coupled DSEs then contain the following undetermined param- 
eters: the renormalization scale fi of the renormalization of the curvature and the renormal- 
ization constant f{fi) of the Coulomb form factor. The parameter d{fjLd) may (or may not) be 
taken care of by implementing the horizon condition. 



3.4 Full numerical solutions 

The renormalized and thus finite DSEs (13.541 13.551 I3.56|) are cast in a form tractable for a 
numerical study. Solving such a coupled set of integral equations is pursued by iteration. 
An educated guess for initial form factors is plugged into the integrals and the results are 
recorded by the Chebychev approximation to be processed further. For asymptotic values of 
momentum, it is instructive for the numerics to make use of the analytical results, cf. the 
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methods used in Refs. |124l 11251 11061 [T26] . The set of nodes cover a finite momentum region. 
In order to extrapolate this to the whole infinite momentum range, the general algebraic forms 
obtained analytically serve as ansatze. However, contrary to what has been previously done 
in Dyson-Schwinger studies, the parameters of these asymptotic forms are still determined 
numerically. E.g., in the infrared it was shown analytically that power laws provide a solution. 
The numerics will have the asymptotic form of the power laws (I2.50p as an input but determine 
the various exponents and coefficients as parameters by nonlinear least-squares fitting. Thus, 
this numerical method provides a check on the analytical results of Table 12.11 rather than 
imposing the latter. 



Previously, a numerical solution of the Coulomb gauge DSEs was found in Ref. The 
infrared behavior agreed approximately with the exponent k = 0.398 of Table 12.11 A heavy 
quark potential was obtained that confines but does not rise linearly, as compared to the 
Wilson loop on the lattice. Therefore, a string tension cannot be extracted from these results. 
After analytical investigations of the infrared behavior of the Coulomb gauge Green functions 
in Ref. |103j . it became clear that yet another solution, with k = ^ should exist. Eventually, 
improved numerical methods were able to indeed confirm k = and reveal a strictly linearly 
rising potential in Ref. [86]. These results will be exhibited in this section. 



The technical details of the numerical calculation are described in Ref. [86]. It shall suffice here 
to mention the essence of the set-up. For convenience, the gauge coupling g is set to unity. The 
functionjfl u>{k) and x{k) are both described in the far infrared by a power law A/k"' where 
A and a are extracted by fitting the solution. It is possible to gain more information on the 
intermediate momentum regime by setting up an infrared expansion of the function 

uj{k) - xik) = c + cik'^ , /c^0,7>0. (3.57) 

and again determining its parameters by fitting. While both uj{k) and xik) are infrared 
enhanced, we have c < oo. Recalling that the 't Hooft loop requires for confinement u = in 



Eq. (I3.53p . we realize that c = 0. 

In the ultraviolet, we use the asymptotic behavior of the form factors found in Ref. ^3] and 
make the corresponding ansatze for fc — > oo, 

uj{k) = k (3.58a) 

xik) ~ k/^\n{k/m^) (3.58b) 

d{k) ~ l/^J\n{k/md) (3.58c) 

f{k) ~ l/^ln(A:/m/) (3.58d) 

extracting the coefficients as well as the different scale parameters m^, m^, nif using least- 
squares fitting. An ultraviolet behavior different from Eq. (I3.58P will be proposed in chapter [H 

All numerical plots displayed in this section were calculated by D. Epple (HSJ- As seen from the 
left panel in Fig. 13.11 both ijj{k) and x(^) are enhanced like 1/A; in the infrared, as predicted in 



^We drop the bar on x{k), still meaning x(fc) in Eq. I|3.36p . 
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0.001 0.01 0.1 1 10 100 1000 0.001 0.01 0.1 1 10 100 1000 



Figure 3.1: Left: The gluon energy uj{k) and the modulus of the scalar 
curvature xik)- Right: The difference uj{k) — x{^)- 

the stochastic vacuum. The right panel of Fig. 13.11 shows the function u}{k) — x(^) vanishing 
for A; ^ 0, in accord with our choice of the renormalization condition c = 0. Fig. 13.21 (left 
panel) shows the infrared enhanced ghost form factor d{k). Thus, the functions d{k), Lj{k) 
and x(fe) are all enhanced as 1/k in the infrared. On the right panel of Figure [321 the heavy 
quark Coulomb potential Vc{r), as given by Eq. (|2.74|) . is shown. An exactly linearly rising 
behavior (within an estimated error of less than one percent) is found, again in agreement 
with the infrared analysis of chapter [2l The linearly rising potential allows us to fit our scale 
from the string tension. Lattice calculations |127l 1751 1128| show, however, that the Coulomb 
string tension ac is about a factor of 1.5... 3 larger than the string tension aw extracted 
from the Wilson loop, in agreement with the analytic result |111| that the Coulomb string 
tension is an upper bound to the Wilson loop string tension. Setting ac = aw therefore 
yields a small quantitative error. However, the qualitative features of the results shown in 
Figs. 13.11 and 13.21 are reliable on the one hand, and they successfully describe the physical 
nonperturbative phenomenon of confinement on the other. Quarks are confined by a linearly 
rising potential Vc{r), and gluons are confined since the dispersion relation uj{k) diverges for 
k ^ 0, prohibiting the propagation of gluons over large distances. 

The numerical confirmation of the infrared behavior predicted by the analytical calculations 
in the stochastic vacuum, see chapter O deserves the following two comments. Firstly, among 
the full structure of the gap equation (I3.54p , the infrared leading contributions were obviously 
already accounted for by the stochastic vacuum. Fig. 13.11 clearly shows that the curvature 
xik), by which the gluon propagator in the stochastic vacuum is expressed over the entire 
momentum range, is the infrared asymptotic function of the (inverse) gluon propagator U!{k). 
In retrospect, these findings may be understood as a confirmation of the Gribov-Zwanziger 
scenario. The second comment is more a mathematical one. Given the full structure of the 
DSEs, how can an asymptotic solution be obtained analytically? We make power law ansatze 
for the infrared behavior and, knowing that in the UV they are not valid, we use them for 
integrating the propagators over the entire momentum range. This method was first used in 
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Figure 3.2: Left: The ghost form factor d{k). Right: The heavy quark 
Coulomb potential Vc{r). 



Landau gauge DSE studies in Ref. |107| and later led to the now accepted infrared exponents 
by the calculations in Refs. [97l I102| . From the experience with the numerical evaluation of 
integrals, this procedure seems acceptable, but a justification would be desirable. In Ref. [54] 
some arguments are given why infrared asymptotics may be computed in the way mentioned 
above. A sketch of a proof was given in the context of the "infrared integral approximation" in 
Ref. [103| . It is also worth mentioning the useful reference |129| in the context of asymptotic 
expansions. Here, we give a short and rather obvious explanation why the method is correct. 
It follows from the Riemann-Lebesgue lemma that convolution integrals C{k) of the type 

C{k) = J dqDi{q)D2{q-k) = J dq J dx J ^i(x)^2(2/) e-^'^'-'^^"^)^ 

dxDi{x)D2i-x)e-''''' (3.59) 



asymptotically approach the function C°-'^{k) obtained by replacing in the integrand the in- 
frared asymptotic functions Df'^{k) = linik^Q Di{k) and Df^{x) = linix^oo Di{x) [i = 1,2), 

C^'ik) = lim C{k) = [ dxbl\x)b^'{-x)e-'^'' 
k^o J 

= j dqD1%q)D^'{q-k) . (3.60) 

This argument can be applied to the (renormalized) DSEs. Replacing the scalar products from 
the Lorentz structure by identities such as 2£ • A; = A;^ + — (£ — A;)^, all (one-loop) DSEs turn 
into a sum of convolution integrals of the type C{k) given in Eq. (|3.59p . 

Recently, there have been investigations of Coulomb gauge Green functions on the lattice 
|131l 11281 [70] . In Fig. 13.31 the lattice data for the ghost propagator and the inverse gluon 
propagator are shown in comparison to the continuum solutions. The qualitative infrared 
enhancement of the ghost propagator is reproduced by the lattice calculations. At the same 
time, the gluon form factor ijj{k) seems to be infrared enhanced as well. It is difficult to extract 
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Figure 3.3: Ghost propagator jJ30y and inverse gluon propagator JZOy 
from lattice calculations, in comparison to the continuum results in the 
Hamiltonian approach. 
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reliable data for the infrared sector on the lattice. Yet, the lattice data for the Coulomb gauge 
gluon propagator of Ref. |131| gives k « 0.49, in agreement with the continuum result. In the 
ultraviolet the striking feature of Fig. 13.31 is that the lattice data for the function uj{k) rise 
stronger than linearly. This might actually be an indication that the ultraviolet properties of 
the continuum solutions need to be improved. We will come back to that in section [431 

Finally, let us mention that in 2 + 1 dimensional Coulomb gauge YM theory, numerical con- 
tinuum results from the Coulomb gauge Hamiltonian approach were recently obtained in Ref. 
|132| . Similarly to the d = 3 case, it is found that all form factors are infrared enhanced, 
with infrared exponents that approximately agree with k = -g, the value given in Table EH] for 
d = 2. 

3.5 Vertices 

The energy density that was minimized in the previous section to give rise to the vacuum Green 
functions was shown to be rather insensitive to the choice of the wave functional, according 
to the discussion in section 13.11 The three-gluon vertex we are turning to in this section will 
show a strong dependence on the wave functional, in particular on the exponent A of the 
Faddeev-Popov determinant. Knowing the infrared behavior of the three-gluon vertex, it will 
be possible to further investigate the ghost-gluon vertex in the infrared. 

Three-gluon vertex 

As pointed out above, it is only the Faddeev-Popov determinant that influences the infrared 
behavior of Yang-Mills theory. The solution obtained for the infrared exponents of the propa- 
gators was found to be independent of the three-gluon vertex, in particular, since it contributes 
to neither the ghost self-energy term nor the ghost loop contribution to the gluon self-energy. 
In the Coulomb gauge, we have seen that any of the vacua ^'a[^] given by Eq. (|3.ip minimizes 
the energy w.r.t. A, evaluated to one-loop order in the DSE (two-loop order in the diagrams for 
the energy). The question is how the three-gluon vertex changes in the infrared for different 
values of A without resorting to the one-loop approximation. 

The full three-gluon vertex is defined as 

{r,,n)tki^,y,z) = {At{x)A'^{y)AUz)) . (3.61) 
In the particular case of A = 1/2 it is found that 

j VAAiAjAk e-^^'^^ = (3.62) 

in a symmetric integration domain. Leaving aside a discussion whether the Gribov region 
is symmetric about A = 0, we here extend the integration domain to the whole gauge-fixed 
configuration space, for practical purposes. Hence, the three-gluon vertex vanishes for ■01/2 
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Figure 3.4: The (complete) DSE for the three-gluon vertex derived from 
the generating functional given in Eq. i3.2\) . 

[83| . Now consider the case A 7^ 1/2. To one-loop order in the gap equation we may employ 
the form (I3.14P for the generating functional Zx[j]. The three-gluon vertex then vanishes for 
any A since the weighting of path integration remains in a Gaussian form. 

On the other hand, without the use of the one-loop approximation, A 7^ 1/2 will give a non- 
zero three-gluon vertex, in contrast to Eq. (|3.62p . as will be shown. Thus, the three-gluon 
vertex shows great sensitivity to the choice of the vacuum wave functional tpx, a behavior not 
exhibited to one-loop order. Making the choice A = permits the standard representation of 
the Faddeev-Popov determinant by ghosts. One can then derive |103| the Dyson-Schwinger 
equation for the Coulomb gauge three-gluon vertex. It is depicted diagrammatically in Fig. 13.41 
Its complete form comprises a diagram with the unknown two-ghost-two-gluon vertex which 
is truncated here. For the following calculation, the ghost-gluon vertex is used at tree-level. 
The DSE for the proper three-gluon vertex then reads |103| 

rijk{ki,k2,h) = N,J d''£DG{i)DG{£ + h)DG{e-k2)r°i{£)r°jie-k2)rl{e + h) , (3.63) 

where the outgoing momenta fcj obey the conservation law 

ki + k2 + k3 = . (3.64) 

The vertex given by Eq. (|3.63l) is projected onto the tensor subspace spanned by the tensor 
components of the tree-level vertex. Due to Bose symmetry, the coefficient functions of these 
six components are all the same, but their signs alternate as the vertex without the color 
structure is antisymmetric under gluon exchange. One finds 

Tijk{ki,k2,k3) = - i{k2)i5jkF{kl,kl,kl) + i{k3)i6jk F{kl, kj, kj) 

+ i{ki),5ikF{kl kl kl) - i{k^)j5ik F{kl kl kl) (3.65) 
- i{ki)k6ijF{kj, kl fc|) + i{k2)kSij F{kl kl kf) . 

Equating Eq. (I3.63P with (|3.65p and contracting with these six tensors, yields a set of six linear 
equations for F, the solution of which reads 
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Figure 3.5: Calculation by M. Leder, published in J J Qgf Dressing func- 
tion of the Coulomb gauge proper three-gluon vertex at the symmetric 
point. The dashed curve shows in contrast the perturbative case where 
the propagators in the loop are tree-level, i.e. k = 0. 



—N f 

{{kf + ki ■ k2){-2J2 -J4 + 3J5) + {k^ + ki ■ k2){2Ji - 3J3 + Je)) (3.66) 

where 

Ji:={ki-i)f J2 ■.={k2 ■ e)e^ J3 ■.={k2 ■ e){ki ■ e) 

J4 :=fc|£2 J5 ~(^k2 ■ ef Je :=(A;i • k2)f . 

The integral (I3.66P depends only on the ghost propagator in this truncation, and despite the 
infrared enhancement of the latter it is convergent. The numerical calculation of the form factor 
F at the symmetric point, where A;^ = A:| = fe| =: A;^, shows a strong infrared enhancement, see 
Fig. 13.51 For the ghost propagator we used the numerical results of Ref. [83j where k = 0.425. 
A fit to the data in Fig. l3.5l yields an infrared power law, F{k'^) ~ {k'^)^^-'^^ . It is clear from the 
discussion on page [72] that this asymptotic power law can be obtained analytically. Replacing 
the integrand in Eq. (I3.66P by its infrared behavior and shifting I X£, one finds from the 
homogeneity that ^(A:^) ~ (/j2)d/2-2-3K_ pj^g 

ging in the value k = 0.425 gives an infrared 
exponent of —1.775 which agrees (within errors) with the numerical result in Fig. 13.51 At large 
momenta the vertex vanishes which complies with asymptotic freedom since in the Gaussian 
vacuum there is no tree-level vertex. 

The above results for the infrared behavior of the Coulomb gauge three-gluon vertex can be 
generalized to any value of k and any dimension d. Therefore, we can also make statements 
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d 






3 
3 

4 


0.398 

1 

2 

0.595 


1 


(^2)1.775 
1 





Table 3.1: Infrared power laws of the function F{k'^). 

about the Landau gauge where d represents the dimension of Euclidean spacetime. Noting 
that d/2 — 2 — 3k = — ac, see Eq. (|2.54|) . we find 

F(/c2) ~ -— ^ . (3.68) 

With the analytical results for k in Coulomb {d = 3) as well as in Landau {d = 4) gauge, 
the exponents in Eq. (|3.68p have the numerical values shown in Table 13.11 The Landau gauge 
result agrees exactly with Ref. [133| . 

Another interesting kinematic point is where one of the gluon momenta, say ki, is set to zero 
while the others remain finite. A technical remark: we switch here to the off-shell gauge 
condition. Thus, no ambiguities are found when trying to evaluate transverse projectors at 
zero momentum. Another technical problem is the following. Trying to calculate the vertex 
with one momentum vanishing from Eq. (|3.65|) by setting ki = 0, the projections onto the 
tensor components fail because the determinant of the coefficient matrix that defines the 
tensor expansion vanishes in this caself] It is advisory to impose /ci = in the DSE (|3.63p . 

Ti,k{0, k, -k) = -ig'^N, j d^i £,(£ - k)j£kDl{i)DG{i - k) . (3.69) 

One can then realize that this integral exists. It can be expanded into a tensor basis con- 
structed by the only scale k and Lorentz invariant tensors, i.e. {ki5jk,kj6ki,ki.6ij}. However, 
the only component that survives the transverse projections of the gluon legs ofTij^ with finite 
momenta, is obviously kidjk- Thus, we can write 

ri,fc(0, k, -k) = -ig^N.kiSjk kmtnpik) j d'^eiminhDlii)DG{i - k) + . . . (3.70) 

^In this context, one might note that the infrared limit of any tensor integral is non-trivial. Given an integral 

we can construct a tensor basis from the external scales {fc'*''} and Lorentz invariant tensors. According to the 
Passarino-Veltman formalism |112| . the above integral can then be expanded in this basis, which is nothing but 
solving a set of linear equations for the expansion coefficients in this basis. If one sets up a tensor expansion 
for finite {fe'*^} and then tries to perform the infrared limit of a single external momentum, say fe^*-* — + 0, the 
coefficient matrix becomes singular, and the tensor expansion is not well-defined. Instead, one can set 
from the beginning (if the integral exists here) and construct the tensor basis spanning a vector space which is 
of a lower dimension than beforehand. The expansion coefficients are then well-defined. 
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where the elhpsis represents irrelevant tensor components which shall be discarded henceforth. 
In the asymptotic infrared, we use the power laws (|2.5Q|) for the propagators and obtain from 
Eq. (irfn]! 



TijkiO, k, -k) = -ig^B^N.hSjkj^^z^ , k^O (3.71) 



with 



/u2\aG-az 

^ + 2k, 1 + k) - H3(2 + 2k,1 + K)/k^) 



1 r(f-2/.)r(f-/.)r(2-f + 3k) 

2(47r)'^/2(d - 1) T{d - 3K)r(l + K)r(2 + 2k) 



(3.72) 



In view of the strong infrared divergence of the three-gluon vertex, one has to check the ghost 
dominance in the propagator DSEs. The infrared power law of the three-gluon vertex (I3.68p . 
expresses that the vertex dressing replaces the infrared exponent of a gluon by that of a ghost 
propagator, for any dimension d. The infrared hierarchy of terms in the gluon DSE remains 
untouched, since even with the dressing of the three-gluon vertex, terms involving it remain 
subleading to the ghost-loop in the infrared. E.g., the gluon loop, which has an infrared 
exponent of d/2 — 2 — 2aA with a tree-level three-gluon vertex, attains an infrared power 
law with the exponent d/2 — 2 — ua — ac if the vertex is dressed. Clearly, this term is still 
subleading w.r.t. the ghost loop which bears an infrared exponent of d/2 — 2 — 2aG- In the 
Landau gauge, the infrared hierarchy was checked systematically by a skeleton expansion in 
Ref. [T33] . 



Ghost-gluon vertex 



The ghost-gluon vertex is investigated here with focus on the (off-shell) Landau gauge. The 
reason is that in this gauge the corresponding DSE is available from Ref. [93]. In the Coulomb 
gauge, with the choice A = for the wave functional \1'a[^], the derivation of the DSE should 
be straightforward, albeit tedious, and technically equivalent to the three-gluon vertex. 

Since neither the DSE studies [HS] nor the lattice calculations [Ml [Ml llOOj show any infrared 
divergences, the dressing function of the ghost-gluon vertex must be some finite function. To 
investigate the consequences of a finite dressing function of the ghost-gluon vertex, let us 
assume, for simplicity, that it is given by a finite constant, 

T^{k;q,p)=CoTl{q) , (3.73) 

The constant Cq can be determined by self-consistently solving the DSE for the ghost-gluon 
vertex in the infrared gluon limit, see section \2M. where T^{q) = igq^j. is the tree-level ghost- 
gluon vertex. Then, the infrared analysis of the propagators can be performed in the same 
way as above. The ghost self-energy and the ghost loop are both multiplied by the constant 
Cq. After evaluation of the integrals, this constant appears on both sides of the equation 
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CqIg = CqIa, cf. Eq. (|2.6ip . and thus trivially cancels. Therefore, a constant dressing of the 
ghost-gluon vertex is completely irrelevant for the infrared behavior of the propagators. 

The question that arises is if a non-constant dressing of the ghost-gluon vertex might result in 
a change for the determining equation (|2.6ip of k. The investigations in Ref. [93] showed that 
after one iteration step of the ghost-gluon vertex DSE, the vertex remains approximately tree- 
level over the whole momentum range, i.e. Cq ~ 1. Also, the results in Ref. [93] confirmed that 
for vanishing incoming ghost momentum p, the ghost- vertex becomes tree-level in the Landau 
gaug^, see Eq. (|2.41l) . Recalling the discussion in section [231 this is expected to hold true 
in the Coulomb gauge. If we discard the irrelevant component of the vertex along the gluon 
momentum k, Tfj_{k;q,p) also becomes tree-level for vanishing outgoing ghost momentum q 
[STllISll, i.e. 

hmT^{k;q,p) = limT^{k;q,p) = Vliq) . (3.74) 



However, the infrared limits of the ghost and gluon momenta are generally not interchangeable 
In particular, zero gluon momentum yields a dressing that is different from unity, although 
quite close to it, as we will see. The following relation shall define C, 

CTliq) = lim (limT^{k;q,p)] + lim Uxn.T ^(k-q,p) \ = rO(g) , g^O. (3.75) 



Does C / 1 or a non-constant C affect the value for the infrared exponent k? To see this, it is 
not necessary to get involved in a numerical calculation but we can argue qualitatively instead. 
Consider any loop integral that involves the ghost-gluon vertex. Wherever it may appear in the 
loop diagram, the ghost-gluon vertex is always attached to ghost propagators. The integrand 
will be strongly enhanced for those loop momenta where the ghost propagator diverges, i.e. for 
J* — > 0. Since the gluon propagator, on the other hand, is finite for all momenta, any infrared 
singularities in the integrand of the loop integral can actually be due to the ghost propagator 
only. The ghost-gluon vertex does not introduce any additional singularities, since its dressing 
function is finite. For the value of the integral, the singularities in the integrand will give the 
dominant contribution. The only value of the dressing function of the ghost-gluon vertex that 
is relevant to the integral, is then the one where any of the ghost momenta vanish. According 
to Eq. (|3.74l) . the vertex is tree-level in these limits. We can therefore infer that in any loop 
integral the tree-level ghost-gluon vertex will yield the correct result. This circumstance can 
thus be traced back to the horizon condition and the transversality of the gluon propagator. 
In Ref. [9^, infrared divergent dressing functions of the ghost-gluon vertex were employed and 
it was shown that k does not vary much. 

Nevertheless, the constant C, defined by Eq. (|3.75|) is not entirely meaningless since the in- 
troduction of a running coupling, see chapter 3] below, makes use of it. One can actually 
analytically calculate C by means of the DSE for the ghost-gluon vertex [93], see Fig. 13.61 

'^This agrees with the corresponding Slavnov-Taylor identity in the Landau gauge. 
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+ / V + 



Figure 3.6: Tie (truncated) Dyson-Schwinger equation for the ghost- 
gluon vertex. 



T,{k;q,p) = rO(g) +rf G^)(A;;q,p) +r(f^^)(fc;(?,p) . (3.76) 

(GGA) 

Here, is a graph with two full ghost and one full gluon propagator in the loop, 

Dg{1 + k)r^(k; e, i + k)DG{i)qata^{e -q){i + k)pl^ , (3.77) 

(GAA) 

and has two gluon and one ghost propagator in the loop, but involves a proper reduced 

three- gluon vertex F^i^p, 

rf^^)(A:;(Z,p) = la^iTliq)DGii-q)Tpii + k;q-i,p) 

Dppie + k)r^,p{k; e, -l - k)D,^{£) . (3.78) 



Comparing Eq. (I3.76P and Eq. (|2.39|1 . the four-point function T^^{k,£;q,p) used in section [231 
can be expanded to yield the two graphs Y^/^'^'^^ and T^i^'^'^\ as well as a proper four-point 
function omitted here. For further details, see Ref. 1931. 



Since T^{k;q,p) exists in the limit /c — > [Ml EHJ [99] , we set k = in the integrands which 
greatly simplifies the tensor structure of Eqs. (|3.77ll and (I3.78p . Furthermore, the proper ghost- 
gluon vertices that appear in the loop integrals are rendered tree-level, as discussed above. We 
then get 

rf«^)(0;<7,g) = tg^C^ ja^ll^D^p{l-q)qpl^Dl{e), (3.79) 



rf^^)(0;g,g) = I <fi q^qpD^,{l)Dpp{£)DG{l - q)T ^,p{Q;£,-l) . (3.80) 



and 



Naively, we would expect from ghost dominance in the infrared that the contribution (I3.80p is 
subdominant since it incorporates only one and not two ghost propagators, like (I3.79p . Using a 
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d 


K 


C 


3 


0.398 


1.089 


3 


1 

2 


1 


4 


0.595 


1.108 



Table 3.2: Solution for the infrared dressing C of the ghost-gluon vertex. 



tree-level three-gluon vertex, we can calculate both integrals for g — > in the infrared integral 
approximation and indeed find that the graph (|3.80p becomes negligible. If the dressed three- 
gluon vertex (|3.7ip is included, the graph (|3.80|) has a substantial contribution to this limit 
of the ghost-gluon vertex. The calculation is then somewhat more involved and deferred to 
appendix O, but one can extract from it the value of C solving the infrared limit of the DSE 
(|3.76|1 self-consistently. For the various solutions of k found in chapter 12.61 the numerical 
values of C are shown in Table 13.2^1 

It is quite remarkable that in the Coulomb gauge with the solution k = ^, the two non-trivial 
graphs that appear in the DSE for the ghost-gluon vertex show an exact mutual cancellation 
in the infrared gluon limit. 




so that the corrections to tree-level vanish and C = 1. Therefore, the inter changeability of 
limits is recovered in this case only and the ghost-gluon vertex becomes tree-level in all infrared 
limitslfl For other values of k and d, the ghost-gluon vertex is not regular in the infrared. It 
must be granted that the omission of the four-point function in the vertex DSE spoils the 
argument. A skeleton expansion of the omitted object would be an interesting project. The 
original derivation of the Landau gauge solution for k in Ref. [97j relied on the regularity of 
the ghost-gluon vertex in the infrared. Their result k = 0.595 for d = 4 does not produce a 
regular ghost-gluon vertex though, the infrared limits are not interchangeable. As long as the 
omitted four-point function has no effect, this would contradict the working hypothesis of Ref. 

m- 



*Note that the results l|3.2p are independent of Nc- The color trace that occurs in the loop diagrams of Eq. 
I|3.76p yields a factor of Nc/2, see Eqs. I|3.77p and l|3.78p . but it cancels with the propagator coefficient term 

g^AB^ = l/(/GiVc). 

^An extension to ghost-antighost symmetric Landau gauge is straightforward and is seen not to alter the 
results for C. 
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3.6 Coulomb form factor 



The success of finding a linear heavy quark potential Vc{r — > oo) ~ r in the Coulomb gauge 
Hamiltonian approach relies on the approximation made for the expectation value {F), see 
section [2771 The Coulomb form factor f{k) was defined by Eq. (|3.24l) . 

and here we try and assess the factorization (F) = k'^DQ{k), i.e. f{k) = 1. To this end, we 
employ the operator identity 

F = ^i9G) (3.83) 

which follows easily from gdG^^ /dg = G^^ — Gq^. We refer to the interesting relation (I3.83p . 
first written down in Ref. [135| . as the Marrero- Swift relation. Taking the expectation value 
of (|3.83p in the Gaussian vacuum and neglecting contributions of 

Plugging the ghost DSE (I2.48P into Eq. (|3.84l) . an integral equation for f{k) is derived, 

f{k) = -|J A _ N.k" (l - {i . kf) DA{i){gDG{i - k))^ (3.85) 

= l + g^NcjdH(l-{i-kf^DAii)Dl{e-k)f{e-k){e-kf (3.86) 

which agrees with the DSE (|3.30|1 for f{k) displayed above. With a power law ansatz for f{k) 
in the infrared, 

/(A:) = ^^, k^O, af>0, (3.87) 

we now aim at a solution of the Coulomb form factor DSE (I3.86p . For A; ^ 0, one finds 

(/fc2)~a/ = {k^)-^fg^AB^NJf{K,af) (3.88) 

where 

{d - 1)k r jaf) r (2k) r {d/2 -2k- 2a f) 

«/j ^^^yn r {d/2 - 2k) r {d/2 - a/ + 1) r (a/ + 2k + 1) ■ ^ ' 



From Eq. (|3.88|1 we can see that 

g^AB'^NJf{K, Of) = 1 (3.90) 

has to be fulfilled. The above relation from the Coulomb form factor DSE can now be plugged 
into the coefficient rule (|2.6Q|) from the ghost DSE to find 

Ig{k) = If{K,aF). (3.91) 
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We now specify the spatial dimension d. For d = 3, Eq. (|3.91l) yields 

_ (-1 + 2K) cos(7r (a/ + 2K))r(4-2a/)r(-2K)r(l + 2a/ + 4/«) sin(7ra/) 
~ vr (-1 + a/) (3 + 2k) {-l + 2af + Ak) T{2 + 2k) ' 



(3.92) 



The numerical solution to this equation gives k as a function of a/, as shown in Fig. 13.71 One 
can see immediately that for any value of a/, the ghost exponent k yields 

K<\. (3.93) 

Therefore, recalling the results for n in Table ITTl there exists no value of a/ for which all three 
DSEs are satisfied. 

It is instructive to focus on the case where 

a/ + 2fi: = l (3.94) 

since this leads directly to a linearly rising potential for static quarks. Plugging this constraint 
into Eq. JSJS]), we find 

^ _ -2 K (3 + 2 k) cos(2 TT k) r(-l - 4 k) r(l + 2 k) 

lo.yo) 



(-1 + 2k) T{-1-2k) 

which has the numerical solution k = 0.245227. Surprisingly, this result is exactly agreed upon 
by lattice calculations [128| . where k = 0.245(5) was found. However, one should notice that 
the lattice calculations carried out so far in Coulomb gauge and also in Landau gauge use too 
small lattices to give reliable results in the infrared. 

In 2 + 1 dimensional Coulomb gauge, the calculations are equivalent. From a simultaneous 
treatment of the ghost DSE and the gluon DSE we find from the stochastic vacuum results in 
section [23] a unique solution for the ghost exponent |102| . 

K = \. (3.96) 

5 

Note that with angular approximation the result is k = 1/4. In a recent publication [132| . the 
value (|3.96p for k was confirmed numerically. On the other hand, if we consider only the ghost 
DSE and the equation for the Coulomb form factor for d = 2, one gets by enforcement of the 
condition Eq. (ITnT]) : 



T{af)T{l - aj -2k)T{k)T{2 + k) _^ 
r(2 - a/) r(-K)^ r(l + a/ + 2 k) 



(3.97) 



The numerical solution is shown in Fig. 13.81 If we require the potential to be linearly rising, the 
solution has to obey a/ + 2k = ^. Eq. ()3.97|1 then leads to the numerical value of k = 0.138. 
Lattice results do not agree with this value |l36j , although they do for 3 + 1 dimensions. 

It is interesting to note that, leaving the DSE for uj aside, a restriction on the infrared exponents 
(|3.93|1 arises from self-consistency of the equations for d and / only. This restriction also 
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Figure 3.7: The ghost exponent Ki djS dj function of Of is shown by a solid 
hne for 3 + 1 dinaensions. A dashed hne indicates the values for which a 
linear potential will arise. 





3. Variational solution for the vacuum state 



concerns the infrared behavior of uj via the sum rule (|2.54p . Consider a fixed wave functional 
with a kernel uj that yields k > | by its infrared behavior. In principle, it should be possible to 
calculate d and / for such a wave functional. However, due to the restriction (|3.93p . no solution 
can be found. We therefore conclude that the Eq. (|3.93p must be due to the approximations 
made in the DSEs for d and /. The relaxation of the horizon condition allows for a self- 
consistent solution of all DSEs where f{k) and d{k) are infrared finite. This will 
be shown in section [3?7l However, we consider the infrared enhancement, of d{k) in particular, 
as the driving mechanism in the Gribov-Zwanziger confinement scenario. Therefore, we infer 
from the above discussion that the DSE for /(fc), as it stands in Eq. (13.8611 . along with the 
resulting restriction (I3.93|l on k, must be ignored. 

Let us reconsider the Marrero-Swift relation (|3.83p with focus on the neglect of In the 
infrared limit, the ghost propagator becomes independent of the gauge coupling as seen in 
the stochastic vacuum, 



d_ ^ d fVAG[gA]DetG-^[gA] 
dg^ ' dg j VAV)eiG-^[gA] 

d jVA' G[A']V)eiG-^[A'] 



fVA' Bet G-^[A'] 



, (3.98) 



where in Eq. (|3.98p the notation reads G ^[gA] = —d^ — gAd. Using the above relation and 
taking a derivative d/dg of the ghost DSE (|2.48l) yields 

= -N,k^ j dH (l - {i ■ fe)2) Dcii - A;)^ {g^DAii)) / . (3.99) 

Neglecting one must run into contradictions. In the expectation value of the Marrero- 

Swift relation (I3.83|l . the neglected terms can be carried along for the stochastic vacuum to 
give 

a_ d jDAgG[A]J[A] 

dg^^ ' dg JVAJ[A] 

JVAJ[A]f^{gG) JVAJ[A]^gG j j^aJ[A] gG f VAJiA]^ 



JVAJ[A] I'DAJIA] !VAJ[A] fVAJ[A] 

(3.100) 



There are two extra terms that were omitted in Eq. (|3.84|) . having the mathematical structure 
of a covariance. This covariance also occurs in the Gaussian vacuum. Since with Eq. (I3.98P 
we have ^ {gG) = (G), the covariance cancels the Coulomb propagator {F) on the r.h.s. of 
Eq. (|3.100|) in favor of a ghost propagator (G), 

(^> + (^sG)-teG>(^) = (G). (3.101) 



The derivation of the Coulomb form factor DSE (I3.86P would therefore fail if the approximation 

dtp 
dg 



— ~ was relaxed. This does not mean that the DSE (|3.86|1 is wrong within the chosen 
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Figure 3.9: Left: The subcritical ghost form factor d{k). Right: The 
Coulomb form factor f{k). Calculated by D. Epple in Ref. jj20y . 



approximation scheme. Applying Wick's theorem to derive an integral equation for f{k) in 
the one-loop approximation also leads to Eq. (|3.86p . However, the above considerations do 
show that, for the asymptotic infrared, the DSE (|3.86p for f{k) describes a higher-order effect. 
On the level of approximations used, we find from Eqs. (|3.84p and (|3.98|1 that 

Plugging this into the Coulomb form factor DSE (I3.30p yields 

= l + if) ,^0 ,3.103) 

which explicitly contradicts the ghost DSE ()3.29|) . note the plus sign. This contradiction is the 
analogon of Eq. HTm for ^ = 0. 

In conclusion, the quality of the approximation scheme must be improved in order to satisfy 
the DSE for the Coulomb form factor. Vertex corrections, which effectively increment the loop 
order, would be an interesting starting point in that direction. 



3.7 Subcritical solutions 

Although it was argued above that in the chosen approximation scheme, the DSE (|3.86p for the 
Coulomb form factor f{k) should be ignored, we here relax the horizon condition to still find 
a solution to all DSEs. The renormalization parameter d^^ := d~^{^d = 0) distinguishes the 
subcritical solutions for d,^^ > and the critical ones for ^ = 0. The subcritical form factors 
d{k) and /(fc), the gap function uj{k), and the curvature x(^) are all infrared finite functions, 
as the numerical calculation shows |120| . In Fig. 13.91 the functions d{k) and f{k) are shown for 
do ^ = 0.1 and d^^ = 0.027. The functions a;(fc) and x(^) are shown in a common plot in Fig. 
13.101 for several values of d^^ . These functions solve their respective DSEs self-consistently. 



87 



3. Variational solution for the vacuum state 



10000 
1000 
100 
10 

1 

0.1 

0.01 

0.001 

1e-04 
1e 



-04 0.001 0.01 




1000 



100 1000 10000 




0.01 



1000 



Figure 3.10: Left: The gluon energy u;{k) and the modulus of the scalar 
curvature, \x{k)\, in the subcritical case. Note that xi^) changes sign at 
k = 1. Internal units are used. Right: The heavy quark potential Vc{r) 
for the subcritical solutions and a linear fit for the intermediate regime. 
Calculated by D. Epple in Ref \TWI . 



The DSE (|3.86p for f{k) is solved as well for subcritical . Lowering the value of dg , the 
solutions are seen to break down. The smallest possible value for which all DSEs are solved 
self-consistently with positive d{k) and f{k) was found to be dp ^ k, 0.02. In the numerical 
calculations, the other renormalization parameters were chosen such that x{^^ = 1) = and 
/(/X = 1) = 1. With the subcritical solutions so obtained, one may calculate the heavy quark 
potential Vc{r), as shown in Fig. 13.101 For intermediate quark separation r, the potential 
Vc{r) is approximately linear; fitting the slope to the lattice string tension, r is found to be in 
the range of hadronic phenomenology. However, for r — > oo, the heavy quark potential becomes 
questionable since permanent confinement is lost for the subcritical solutions. In addition, the 
running coupling a{k) vanishes in the infrared, as seen in the forthcoming chapter. Altogether, 
the implications of the subcritical solutions for the infrared physics are not compatible with 
Coulomb confinement. This indicates that the horizon condition must not be abandoned. 
Despite the inconsistency in the DSE for /(fc), which is due to higher-order effects, we favor 
the solutions with critical dg ^ = 0, i.e. with the horizon condition implemented where the 
functions a;(/c), x{k) and d{k) are infrared enhanced. 
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4 Running coupling and ultraviolet behavior 

A quantitatively correct treatment of QCD, with the possibility to describe phenomenological 
effects, needs to take into account accurately its scale Aqcd- In the absence of flavor charges, 
Nf = 0, this is the only scale of the theory. The value of Aqcd can be related to the 
perturbative expansion of the (Gell-Mann-Low) beta function, 

Pia) = ^ = -/3o5' - - 029' + 0{g^) (4.1) 

that describes the dependence of the renormalized gauge coupling g on the renormalization 
scale /X. Solving the differential equation (|4.1|) with neglect of 0{g^), the gauge coupling is 
seen to diverge at the scale 

This breakdown of perturbation theory occurs at a finite Aqcd, if 5^(m) > 0: giving rise to 
dimensional transmutation of the dimensionless gauge coupling to the dimensionful energy 
scale Aqcd- The 0(5^) coefficient of the beta function M.ljl . /3o, is a gauge invariant quantity 
that is calculated in perturbation theory to be 

Higher order coefficients of the beta function were obtained in cumbersome calculations |137^ 
I138| and serve for an improved matching of Aqcd to experiment. 

High-energy experiments are able to provide us with the scale Aqcd immanent in QCdQ A 
physical process that involves measurable quantities, e.g. the mass of the Z boson, needs to be 
described by theoretical calculations which are then compared to the experimental data. The 
corresponding scattering matrix is computed by means of perturbative QCD, preferably in a 
covariant gauge, in a certain chosen renormalization scheme. As a result, the running coupling 

aik) = ^ (4.4) 

can be extracted for the momentum k of the given process. For the mass Mz = 91.19 GeV of 
the Z boson, it is found that a{Mz) = 0.1187 in the MS scheme and by the extrapolation of 
the perturbative a{k) into the infraredj^ 



/3oln(^7A^CD) 



In(^VA^c^) 



^2 ln(/.VA^cD) P^ln\f,yAlc^) 



2 



■2 



0t 



(4.5) 



^Estimates can also be obtained directly by lattice simulations |139) . 

^Various conventions exist for the definition of the /3 function and its coefficients of the perturbative ex- 
pansion. The barred symbols in Eq. I|4.5p are related to those of Eq. I|4.ip by /3o ~ ff^yi , Pi ~ (ffjT and 



P2 - 
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4. Running coupling and ultraviolet behavior 



one finds Ajg^ = 217 MeV [14Q| . They do, however, depend on the renormalization prescrip- 
tion, hence the subscript on A.jjg, as an artefact of perturbation theory. A change in the 
renormalization prescription can in fact be compensated for by a finite renormahzation group 
transformation. It is the truncation of the perturbative series that makes the running couphng 
scheme dependent |141| . 

The variational solutions presented in the previous chapter only loosely account for the ultra- 
violet properties of the theory and it will be shown that the /? function turns out to be slightly 
incorrect. The first few coefficients of the beta function, in particular /5o, record inaccuracies of 
the solutions and any observable is immensely sensitive to such errors, as can be seen directly 
in Eq. (14.21) . If the nonperturbative variational solutions are to give quantitatively sensible 
results, they have to reproduce /3o with high accuracy. In this chapter, we first recall the crucial 
properties of the YM coupling constant and then define a nonperturbative running coupling 
to come by some information on how the Green functions are to behave in the ultraviolet. 

4.1 Anti-screening of color charges 

An effect specific to the non-abelian nature of YM theory is the anti- screening of color charges. 
The interaction of two color charges on a short distance can be treated by perturbation theory. 
The perturbative corrections increase the magnitude of the (renormalized) color charges, i.e. 
the potential is "anti-screened", contrary to the abelian theory. Anti-screening is a crucial 
property of YM theory and intimately connected to asymptotic freedom. 

One may find the potential between two color charges by setting up a basis of states 

in the absence of external charges, then turning on pext, that is p = pdyn — > ^dyn + Pext, and 

calculating the correction to the potential energy Ec to second order in PextP 



The screening vacuum polarization is of course still present as the manifestly negative contri- 
bution E^ in (|4.6p . Its calculation in a perturbative expansion in the bare gauge coupling 
gB to order is familiar from the vacuum polarization in the abelian theory [^. To let 

two point-like charges, p^^^{x) as in (|2.72p . approach each other from an infinite separation 
to r, the energy Ec\oo — Ec\r = / d^kVc{k) e^'^''^ is needed. The momentum space potential 

^The Faddeev-Popov determinants J are omitted since they do not contribute in the lowest-order pertur- 
bation in gB. 




(4.6) 
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4.1. Anti-screening of color charges 




Figure 4.1: Perturbative screening and anti-screening contributions to 
the color Coulomb potential. The first graph represents the usual 
Coulomb interaction, the second graph shows the genuinely nonabelian 
anti-screening contribution and the last graph corresponds to the usual 
screening term. These graphs are not Feynman diagrams. 



yc{k) = VQ^\k) + V^^ \k) comprises the screening contribution 

where the cut-off A regularizes the UV divergence and the quadratic Casimir yields here for 
SU{Nc) the value C2 = A'c- The screening effect is, however, over-compensated by the con- 
tribution V^'^\k). To see that, expand the non-abelian Coulomb operator around its abelian 
part, 

F'''{x,y) = (Go + 2GogBAdGo + 3GogBAdGogBAdGo + 0{g%)y\x,y) , (4.8) 

and then take the expectation value, see e.g. Ref. tl42j . where the third term in (|4.8p strength- 
ens the potential. 



Note that in the abelian theory only the first term in the perturbative expansion (|4.8p ex- 
ists, and thus vl^~^\k) = g'^/k'^. The anti-screening correction in Eq. (|4.9p is genuinely non- 
abelian. In Fig. 14.11 the neatly separated screening and anti-screening contributions are shown 
diagrammatically. The net effect of the 0(5^) contributions increases the potential between 
static SU{Nc) color charges, 

as first demonstrated by Khriplovich for Nc = 2 [58]. The net anti-screening correction to 
the static quark potential can be used to define the renormalized coupling g{fJ.), with the 
ambiguity in subtracting the divergence inherent in the renormalization scale From the 
so-defined g{fi), the (3 function may be calculated via Eq. (|4.1|) . reproducing its Nobel-prize 
awarded perturbative expansion (with the interpretation of asymptotic freedom |143l 1144] ). as 
shown in Refs. [62l[Tl5] . 

If we wish to extract correctly from the color Coulomb potential the /3 function for the vari- 
ational solutions of chapter [31 the form factors are to behave appropriately in the ultraviolet. 
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4. Running coupling and ultraviolet behavior 



What exactly "appropriately" means will become clear in the course of this chapter. For 
starters, consider the Coulomb potential 

vt'\k) = '^d\k)m (4.11) 

as the part of the energy in the variational state that depends on the separation of the charges, 
see Eq. (|2.74p . The self- energies, which only regulate infrared divergences, are omitted. To 
lowest order in the bare gauge coupling gB^ d{k) = f{k) = 1. The 0{gj^) term of Eq. (|4.1ip 
can be found by expanding the DSEs for d (|3.29|1 and for / (jS.SOp up to 0{g\), using the 
lowest-order expression uj{k) = /c, 

d{k) = l + g|^ jd'l{l-{i-kf) ^^^^^^^ (4.12) 







■ kf) 


ld^i(i- 


{!■ 





m = ^ + 9l^J^^{^-(^-kf)j^l^ (4.13) 

Note that d{k) and f{k) have the identical perturbative 0{g'^) expansion. Plugging these into 
Eq. (|4.11|) . yields three times the same ©(g^) term, 

Vt'\k) = f (l + 3 • 5lf /^^ (l - {i • kf) ) + 0{g%) , (4.14) 



cf. the factor of 3 in the expansion of the non-abelian Coulomb operator (|4.8p . Regularized by 
a UV cut-off A, the integral in Eq. (|4.14D to 0{g%) gives 



in agreement with Eq. (14. 9|) . Note, however, that only the anti-screening contribution VQ^\k) 
is thus accounted for. Determining the state l^*) by the variational principle in the absence 
of color charges and deriving the color Coulomb potential Vc{r) as a perturbation around 
the vacuum, necessitates that the screening contribution \k) is carried along in the full 
second-order perturbative expression as in Eq. (|4.6p . 



The renormalized coupling g is defined so as to render the expression for Vc{k) finite, order 
by order in perturbation theory. Decomposing an infinite quantity into a finite and an infinite 
part always leaves an ambiguity that is controlled by the renormalization scale /i. Using 

gB{A) = Zg{A,fi)g{f,) (4.16) 

as the definition of the renormalized, ;U-dependent g{fJ-), the expression V^'"''\k) in Eq. Km 
turns out finite up to 0{g'^) if the renormalization constant is chosen acl 



1^12iV, A2 
2 (47r)2 3 //2 



Zg{A,^,) = 1 - - . (4.17) 



""Swapping qb for g does not matter to the order considered. 



92 



4.2. NONPERTURBATIVE RUNNING COUPLING 



The P function, as defined by Eq. ()4.1I) . can now be calculated from Eqs. ()4.16|1 and ()4.17|) to 
lowest order, 



This result does not agree with the canonical value of /Jq, see Eq. (|4.3p . but it is clear why: 
only the contributions of longitudinal gluons were considered to order 0{g^), the screening 
contribution from transversal gluons was omitted. For the same reason, the perturbative 
renormalization group calculation of the color Coulomb potential [146] relates Vc{k) to the 
running coupling g'^{k) by k'^Vc{k) = xog'^{k) with xq = jf . In the next sections, we aim to 
use a nonperturbative definition of the running coupling. 



4.2 Nonperturbative running coupling 



Since the variational Coulomb gauge solution from chapter [3] is nonperturbative, we seek here 
a nonperturbative definition of the running coupling in order to match a{Mz) = 0.1187 to our 
solutions. In the Coulomb gauge, there are two possibilities for its definition, either via the color 
Coulomb potential, or by means of the ghost-gluon vertex [lOlj . The nonrenormalization of 
the ghost-gluon vertex has proven above to be a useful common property of both the Coulomb 
and Landau gauge. Here, it is first discussed how the Landau gauge ghost-gluon vertex may 
define a nonperturbative running coupling and afterwards we turn to the Coulomb gauge. 



Landau gauge running coupling from ghost-gluon vertex 

The local off-shell formulation of Landau gauge YM theory allows for multiplicative renor- 
malization of the fields by the introduction of renormalization constants Z(A, /x). Some of 
these are related to one another by the fact that Zg, the renormalization constant of the gauge 
coupling, can be defined by any of the interaction terms in the YM Lagrangian, namely the 
ghost-gluon vertex, the quark-gluon vertex, the three- and the four-gluon vertex. The equality 
of the gauge coupling for all these interactions is expressed by the Slavnov-Taylor identity 

3/2 ... ^ ~ 

where Z\ = ZgZ^ is the renormalization constant of the three-gluon vertex, Zi = ZgZ^Z-^ 
the one of the ghost-gluon vertex, while Z3 and Z3 renormalize gluon and ghost fields, 

Ab = zl^^A , CB = Zll'^c , CB = Zl^'^c . (4.20) 

An unambiguous definition of the Z's requires renormalization prescriptions. We define the 
theory at the scale /i by setting 

DA{k = /i; ^) = ^ , Dcik = is;fi) = ^. (4.21) 



93 



4. Running coupling and ultraviolet behavior 



The nonrenormalization of the ghost-gluon vertex in Landau gauge infers that Z\ is a finite 
number, independent of ii. Therefore the quantity 

g\{k)DA,B{k-k)Dlj,(k;k) = Zfg\fi)DAik;fi)Dlik;ti) (4.22) 

is finite, due to Zi < oo, and ^-independent, since the l.h.s. in ^-independent. One may set 
Zi = 1 but any other finite number will also do in an appropriate renormalization prescrip- 
tion, recall section 12.51 It is left arbitrary here. In the perturbative momentum subtraction 
scheme, the values Ca = Cq = 1 define the renormalized quantities on the r.h.s. of Eq. (|4.22p . 
Nonperturbative Green functions might not allow for such a simple subtraction scheme. We 
therefore generalize it [53] by setting 

CaCI = 1 . (4.23) 

Other schemes can be accomplished by performing a finite renormalization, see below. 

In Eq. (|4.22l) . one may choose freely. Leaving /x once unspecified and setting it once to fj, = k 
leads with (|4.21l) to a nonperturbative definition of the running coupling [53j . 

g\k) = k^g\p)DA{k; fi)Dl{k; fi) . (4.24) 

While the factors on the r.h.s. are //-dependent, the running coupling g^{k) on the l.h.s. is 
renormalization scale- and scheme independent. Note that the finite value Zi drops out. 

1/2 1/2 

Applying a finite renormalization such as ^ ^ Zg A, c ^ Zg c in order to discuss the case 
CaCq / 1, the running coupling can be expressed by the transformed quantities, 

g\k) = Zl k'^g'{ii)DA(k-ii)Dl(k-^) . (4.25) 

where we have set Zfz^z^ Zf. Note that these simple scalings are not renormalization 
group transformations in the usual sense since they do not depend on the scale /i. It will be 
shown explicitly below that the (asymptotic) behavior of the running coupling is ignorant to 
such scalings. 

A different point of view is frequently found in the literature, see e.g. Ref. |124| . One may 
always absorb finite transformations by a transformation in g such that Zf = ZgZ^Z^ is left 
unchanged (usually chosen to be unity). The running coupling (|4.24l) then stays form-invariant. 
Here, we prefer to leave g unchanged and allow for a change in Zi and use Eq. (I4.25P as the 
running coupling. 

Coulomb gauge running coupling from ghost-gluon vertex 

A running coupling in Coulomb gauge may be defined in analogy to the Landau gauge from 
the ghost-gluon vertex. The multiplicative renormalization, however, is not as well-understood 
in the Coulomb gauge. From within a class of interpolating gauges where multiplicative renor- 
malization was discussed in Ref. [147] . the Coulomb gauge limit may be taken, breaking the 
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4.2. NONPERTURBATIVE RUNNING COUPLING 



path for a definition similar to Eq. (|4.25p in the Landau gauge |101| . We adopt this definition 
here. 

Using the notation d{k) for the ghost form factor (I3.23P and a;(fc) for the equal-time gluon 
propagator (13. 5p . both renormalized multiplicatively, 

dB{k; A) = Z3(A, fi)d(k; fi) , conik, A) = Z^iA, fi)io{k; fi) (4.26) 
d{k = fi; = Cd , uj{k = ^) = C^k , (4.27) 

with C(7^C| = 1, we find from the ghost-gluon vertex analogously to Eq. (|4.22p . 

gliA)dlik;A)^u;^Hk;A) = Zfg\ti)d\k; f,)^u;-\k; f,) . (4.28) 

Allowing for a change in the renormalization prescription, C~^C'^ ^ 1, one may either absorb 
it by a change in (7, or (as above) redefine Z\. Accordingly to (|4.25ll . the latter choice leads 
directly to the running coupling 

g^{k) = Zlkg^{li)d^{k,p.)^-\k,ii) (4.29) 

Note that k is here k = \h\ and all interactions are instantaneous!! The running coupling 
(|4.29p was given in Ref. |1Q1| with a prefactor of |, chosen as to recover the Landau gauge 
infrared limit from within the interpolating gauge. According to the discussion above, this 
factor needs to be interpreted as = |. 



Coulomb gauge running coupling from the color Coulomb potential 

Another renormalization group invariant that qualifies to define a nonperturbative running 
coupling is the instantaneous color Coulomb potential Vq^)- It was shown in Ref. [72] that the 
product gA^ is a renormalization group invariant and in Ref. |146| how the quantity (17^ ^0^0) 
actually relates to Vc'(fc). We may accordingly define the running coupling in Coulomb gauge 
also by 

g\k)=^^g^{ix)k^Vc{k;ii), k = \k\. (4.30) 

The factor of accounts for the fact that Vcik) only comprises the anti-screening contribution 
and over-emphasizes this effect. While in the ultraviolet the definition ()4.30|1 is expected to 
agree with the definition (I4.29P of the running coupling from the ghost-gluon vertex, the 
infrared behavior turns out to disagree drastically. In view of the confining long-range part 
of Vc'(fe), the running coupling in Eq. (14.30(1 diverges for fc — > 0, a clear instance of infrared 
slavery. On the other hand, the definition (I4.29P will yield a running coupling that freezes in 
the infrared, see below. 

'For equal times, the tree-level gluon propagator is D\{]i) = see Eq. I|3.3ip . Therefore, the definition 
l|4.29p is larger than in Landau gauge l|4.25p by a factor of 2. 
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4. Running coupling and ultraviolet behavior 



4.3 Ultraviolet behavior of Green functions 



Given relations that express the running couphng in terms of the Green functions nonpertur- 
batively, one may extract information on these Green functions in the UV domain, where the 
running couphng is known from perturbation theory. Solving the 0{g^) expansion of the /3 
function (|4.ip as a differential equation for g'^{k) yields for large /c^ = k^k^ 

a'C") = ■ (4-31) 

The ultraviolet behavior of the nonperturbative propagators has to be such that the non- 
perturbative running couplings defined in section W?2\ approach the function (I4.3ip . Further 
information can be extracted by solving the ghost Dyson-Schwinger equation self-consistently 
with an ansatz for the ultraviolet behavior. This ansatz requires an educated guess. 

Despite asymptotic freedom, the Green functions do not asymptotically become tree-level in 
the ultraviolet. Weinberg's theorem states that there can be logarithmic corrections |148| . and 
an ansatz for the gluon and ghost propagators for /c » /i might therefore read 



1 

A;2 In'' 



Dcik) ~ TTTTYTk^, (4-32b) 



The exponents of the logarithms are specified by the renormalization group in relation to the 
anomalous dimensions of the gluon and the ghost field |143j . 

7a(5) = ^ = + (4.33a) 

7a(5) = ^ = 1%9' + 0{g') (4.33b) 



7 = 




13 




~ 22 


5 = 


7g 


9 




/So 


~ 44 



In the Landau gauge, this relation is given by [143j 

13 

(4.34a) 
(4.34b) 

Oftentimes, the logarithmic exponents 7 and 5 in (|4.32l) are referred to as anomalous dimensions 
of the gluon and the ghost, although strictly speaking they are only the lowest order coefficients 
of the anomalous dimensions 7^1 and 7g defined in Eq. (|4.33p and divided by Pq. 

We now try to reproduce the well-established results (I4.34p from our nonperturbative study of 
Dyson-Schwinger equations in Landau gauge. In the same procedure, we may then go on to 
determine the UV behavior of Coulomb gauge propagators. 
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4.3. Ultraviolet behavior of Green functions 



Landau gauge anomalous dimensions 

Consider the (d = 4) Landau gauge ghost DSE multipHcatively renormahzed in the scheme 
(I4.2ip . where ^^DaI/w) = Ca, fJ-'^DcifJ') = Cq and the ghost-gluon vertex is evaluated at some 
momentum configuration to define a value of Zi, 

D^\k) = k'^Zs - g^N^Zlk^ j d^i (l - {! ■ kf'^ DA{i)DGii - k) . (4.35) 

A renormalization prescription that does not agree with (I4.23P is accounted for by a change in 
Zl. In order to study the ghost DSE (I4.35P in the asymptotic ultraviolet limit, the following 
ansatze are made: 

A;2(l + %ln^(^)) kHvC^i^,) 



A:2(l + %ln^(4)) A;2ln^(^) 



These satisfy the renormalization prescriptions (|4.2ip f or k = fi and behave like the proposed 
ultraviolet behavior (|4.32l) for A: ^ oo, with coefficients A and B. Rewrite the ghost DSE 
(OBll by 

1 



k^Dcik) 



Zs-g'N,ZtABlG{k) , (4.37) 



where the integral Tdk) can then be asymptotically evaluated for /c ^ /i with a dimensionally 
regularized UV divergencelf] 



^ oo) = / ct^-^H ^ } \^ = I 



y+5 I 



oo 



; drism^ r]{l + Oie)) I ^ ^ — r~7V ) x := eln ( ^ 

(27r)3 7 1 {n'^n^^^(^) 

oo 

3 1 1 



327r2 



(l + 0(6))(/x2)-^-^i I dxx-^-'e-^ 



172 ) 

2 



647r2^ ^ V £1-7-5 1-7- J 



(47r)2 4e'^-J-S (47^)2 4(1 _ ^ _ 



^One can integrate tor £ < k and £ > k and expand in ^/fc and fc/£, respectively. Only one term is 
dominant which turns out to be the one found by the "angular approximation". We restrict the calculation to 
0<1 — 7 — (5<1. Furthermore, the expansion of the incomplete gamma function 



J A ^ 



(-l)"r" 



n=0 -!(^ + ") 



is used. 
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Terms that vanish for e — > were discarded. Matching this result to the l.h.s. of the ghost 
DSE (|4.37p leads to a sum rule of the anomalous dimensions, 

7 + 25 = 1 , (4.39) 

as well as the relation 

The infinite part of the renormalization constant is identified by 

^3 = ^^^ + finite (4.41) 

where setting the finite part to zero puts QCD into a nonperturbative phaselll 

We may now plug the asymptotic forms (|4.36l) into the running coupling (|4.25l) to obtain with 
Eq. (lOnH 

,.(.^„,.5,y_l^.__l^ (4.42) 



Note that the finite value of Zi specific to the renormalization scheme drops out. The above UV 
limit expression for the nonperturbative g'^{k) can be compared to the perturbative expression 
(;i:3TD to yield 

- ' ' (4.43) 



(47r)2 4/3o 44 ' 
in agreement with the perturbative result for 6 in Eq. (I4.34p . 



Coulomb gauge anomalous dimensions 



Since the above UV analysis of the nonperturbative ghost Dyson-Schwinger equation success- 
fully yielded the correct anomalous dimensions for the Landau gauge, we use the same method 
to extract the anomalous dimensions for the gluon and ghost fields in the Coulomb gauge. 



The multiplicatively renormalized ghost DSE reads in terms of d{k) and io{k), 

d{£ - k) 



d-\k) 



g^N.Zl dH kf 



and the ansatze 



UJ 



-Hk) 



a 



-1 



k 



1 + S^ln7(g 



d{k) 



2uj{e){e-ky 



It was discussed in Ref. [78] that the horizon condition does not allow Z3 
contradicts perturbation theory. 



(4.44) 



(4.45) 



1 for g ^ and hence 
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give for k = iJ, the values and C^, resp., whereas the k —>■ oo behavior is determined by 
the coefficients A and B and the exponents 7, 5, similarly to the Landau gauge. Setting the 
renormalization constant to 

~ £(5 + 1)1 

Z3 = ^ ^ , (4.46) 

enforces the horizon condition (cf. Eg. I4.4ip . and one may calculate the integral in Eq. (|4.45ll 

for A; ^ 00 by 



-1-^ I 4) = Z,-g^Njf I d^~-^l (1 - (i.kf) ^ 



where we have matched the exponents by means of the sum rule 

7 + 25 = 1 , (4.48) 
and the coefficients are found to to obey 

With these results, the Coulomb gauge running coupling (|4.29l) gives for the ultraviolet limit 
9\k ^ 00) = Zlg-^^ = J (4.50) 

and in comparison to its known gauge invariant behavior (I4.3ip this yields 

1 4iV,, 



(47r)2 3/3o ■ 



(4.51) 



With the canonical value (|4.3p for (3o and by virtue of the sum rule (|4.48|) . the logarithmic 
exponents of ghost and gluon propagators can thus be found to yield 

7 = ^, * = ^. (4.52) 

The lowest order anomalous dimensions (cf. Eq. I4.34P 7^ = /3o7 = and = (3o5 = 
agree with the purely perturbative calculation [56l 11491 172] . Therefore, a full numerical 
nonperturbative solution is expected to yield the ultraviolet behavior specified by (|4.45|) with 
the logarithmic exponents (|4.52|) . 
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4. Running coupling and ultraviolet behavior 



Recalling the numerical solution from chapter [HI the values for 7 and 6, see Eq. (|3.58l) . were 

7num = , 6num = ^ • (4.53) 

Note that the sum rule (|4.48p is satisfied by the above values. In the calculation that led to 
(|4.52|) . we required /So to be the canonical value. We may turn the argument around, assume 
the numerical solutions (I4.53P and calculate /3o from the UV behavior of the resulting running 
coupling, cf. Eq. (|4.51|) . This procedure yields 

" (4vr)2 35„,„ - (4vr)2 3 " ^^'^^^ 

The disagreement of /Sg"™ with (3q is small, but in view of the exponential dependence of 
the QCD scale on /3o, a slightly erroneous result for Pq will make a sensible assignment of 
^QCD hopeless. The failure to find the correct Po is clearly due to the approximations made. 
The Gaussian wave functionals used eliminate the gluon loop in the gap equation since there 
exists no tree-level three-gluon vertex. From purely perturbative calculations it is known that 
the gluon loop is essential for the gluon self-energy to preserve gauge invariance and produce 
the correct Pq. A nonperturbative calculation therefore requires a truncation that involves 
the gluon loop as well. However, DSE studies in Landau gauge showed that the dressing of 
the three-gluon vertex in the gluon loop must be highly non-trivial [53l I124| . and to achieve 
results with the canonical Pq value a momentum dependent renormalization constant Zi is 
proposed [5l]- In the Hamiltonian approach with a wave functional, the extraction of the 
right Pq was pursued from the Coulomb interaction using perturbative UV tails for the form 
factors [1501 1151] . It is not clear at present how to manipulate the wave functional such that 
the nonperturbative UV behavior (I4.45P with correct values for 6 and 7 is reproduced from 
analyzing the gap equation in the Coulomb gauge Hamiltonian approach. 

The agreement of the result (I4.54p with Refs. [1521 [83] is coincidental. There, after factoring 
out qb from the ghost form factor, the P function was defined as the (logarithmic) derivative 
of the tree-level term in the ghost equation, being g^^{A). However, this term is actually Z3, 
see Eq. (14.441) . What was proposed as the P function was actually 

P=j^lnZs = 2-fd (4.55) 
am/i 

i.e. twice the anomalous dimension of the ghost form factor. With the specific value of 6 = 1/2, 
these objects are indeed identical, since 5 = 

An alternative way to calculate the P function from the nonperturbative variational solutions in 
the Coulomb gauge is to extract information from the heavy quark potential Vc{k). Since the 
nonperturbative running coupling can also be defined via Vc{k), see Eq. (I4.30p . the Coulomb 
form factor f{k) must yield 

Vcik) = ~ ^ . (4.56) 
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4.4. Infrared fixed point 



Making an ansatz for the Coulomb form factor f{k) in the ultraviolet, 

C 



f{k) = , (4.57) 

therefore provides a sum rule for the logarithmic exponents b and 99, 

2(5 + (/? = 1 . (4.58) 

Furthermore, the UV ansatz (|4.57|) for the form factor /(/c) can be plugged into its DSE (|3.30ll 
to yield after UV analysijfl 

9'ivJ?^B^^A . 1 . (4.59) 

By comparison of the above relation to the identity (|4.49p from the UV analysis of the ghost 
DSE, we infer that f = 5 and with the sum rules ()4.58|1 and (|4.48p this yields 

6 = ip = j=^. (4.60) 

One may now use Eq. (|4.51|) from the ghost-gluon vertex to calculate (3o as a function of 6, 

The occurrence of the anti-screening contribution is familiar from the consideration of 

the instantaneous potential, see section 14.11 Apparently, the non-instantaneous part needs to 
be taken into account, possibly in the manner proposed in Ref. |146| . 

Apart from the technical difficulties to arrive at the correct value for /?o automatically, the 
above considerations showed that it is possible to extract a certain ultraviolet behavior of the 
Coulomb gauge Green functions, having imposed the known value for Pq. It may therefore be 
expected that, once approximations can be relaxed, the ultraviolet behavior according to Eqs. 
and Km . i.e. 

Lo{k)^kln'/''k, d{k)^-^—, (4.62) 

in ' k 

will be found numerically. Lattice calculations |128| claimed a counter-intuitive UV power 
law UJ ~ A;^/^. The data displayed in Fig. 13.31 also indicate an UV enhancement of uj{k). 
However, a stronger than linear rising of uj{k) might as well be an indication of the logarithmic 
correction in (|4.62l) . The lattice calculations in [70] agree with such a statement. This should 
be investigated further. 



4.4 Infrared fixed point 

In this section we investigate the infrared behavior of the Coulomb gauge running coupling as 
defined by Eq. (|4.29p . Due to the nonrenormalization property of the ghost-gluon vertex, the 

*In the same manner as for the infrared analysis in section 13.61 the asymptotic evaluation of the DSE for 
f{k) gives no additional condition on the (power or logarithmic) exponents. This is due to the fact the product 
i^i^'^l^ behaves like ^ in the infrared and like -^j^ in the ultraviolet, as a consequence of the ghost-gluon vertex' 
nonrenormalization. 
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0.01 
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Figure 4.2: The running coupling a{k) as defined by Eq. Ii4.29\) with the 
variational solutions presented in chapter \^ 



running coupling attains a finite value in the infrared. The corresponding (3 function then has 
a fixed point in the infrared. This can be readily seen by noting that the sum rule (I2.54p for 
infrared exponents yields 



1 - 2(2aG + 1 + a^) 







(4.63) 



in the Coulomb gauge where d = 3. Hence, with the infrared power laws (|2.50p with coefficients 
2 A for uj~^{k) and B for d{k), we find for the infrared limit of the (equal-time) running coupling 



a{k) 



in Eq. K29\\ 



a. 



lim aik) = Z^q 



2tt NJcin) ' 



(4.64) 



where we have made use of Eq. (|2.60l) which relates the infrared coefficients A and B to the 
function Ig{h) defined by Eq. (I2.59hp l^ For the preferred solution with n 
IciK) to find 



^, we evaluate 



a, 



.1 

" 2 

cld=3 



2tt 



(4.65) 



In Fig. 14.21 the running coupling as given in Eq. (|4.29|) is shown with the full numerical solutions 
for d{k) and u;{k) given in chapter [3l The numerical value of Uc is in excellent agreement with 
Eq. g35|)0 



®A factor of needs to be reintroduced here, it was set to unity in the infrared analysis of chapter [2] 
^"in Ref. [86], the definition of the running coupling from Ref. [lOlj with a factor of | in Eq. I|4.29p was used 
without adjustment of Z\. Therefore, the value of ac was found to be higher by that factor, i.e. 167r/(3A^c). 
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The other critical solutions for various spatial dimensions d and the respective solutions for K{d) 
yield different values for the infrared limit of the running coupling. To study the dependence 
on the dimension d, we use here the Landau gauge running coupling a{k) := with g'^{k) 
defined in Eq. (|4.25p . The infrared limit yields 

die ■= lim a(k) = - — ^ , , , (4.66) 

cf. Eq. (|4.64|) . For the three solution branches Ka{d), Kb{d) and Kc{d) shown in Fig. 12.51 the 
value of ac is shown in Fig. 14.31 as a function of the dimension d. It is seen that for large d, 
add) rises exponentially. For d ^ 2, which corresponds to the Landau gauge 1 + 1 dimensional 
YM theory, the solutions Ka{d) and Kc{d) yield add = 2) = whereas for Ubid) the infrared 
limit etc approaches a finite value 



28/5^2(6)^(13 



ac\:j ■ N, = -^—^§^--fLLN, « 0.501iVe (4.67) 



The root Kb{d) does not show any qualitative change for d = 2 where all Landau gauge degrees 
of freedom are frozen. The solutions Ka{d) and Kc{d), on the other hand, do show a qualitative 
change at d = 2. Granted the approximations made, the latter solutions with Oc = for 
d = 2 indicate what one may expect: the 1 + 1 dimensional Landau gauge calls for a separate 
treatment. 

Returning to the d = 3 Coulomb gauge, it is reassuring to see that the running coupling in Fig. 
14.21 is a monotonic function. Otherwise the (3 function would have spurious zeros (a problem 
that occurs in the Landau gauge [106]). The infrared finite value of the running coupling infers 
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|3(a) 




a 



Figure 4.4: Sketch of the (5 function. The sohd hne quahtatively describes 
the behavior of the running couphng as shown in Fig. \4.2l The dashed 
hne shows the case where there is a small bump in the running coupling. 
With a dotted line, the (3 function is shown for the subcritical solution. 

that the (3 function has an infrared fixed point for |^ = In Fig. 14.41 several cases of 

the (5 function are sketched. The subcritical solutions (section [377|l . as opposed to the critical 
solutions, do not have a sum rule of infrared exponents since all form factors are infrared finite. 
The running coupling a{k) should have an ultraviolet behavior similar to the one shown in 
Fig. 14.21 since the form factor are qualitatively the same in the UV. However, in the infrared 
one finds 

a(fc ^ 0) ~ A; , (4.68) 

that is the running coupling from the ghost-gluon vertex vanishes in the infrared for the 
subcritical solution. This behavior is unacceptable for a theory that is strongly coupled in 
the infrared. In contrary to the 1 + 1 dimensional theory, the 3+1 dimensional case requires 
a non- vanishing coupling strength for the infrared degrees of freedom. Moreover, if we hold 
onto asymptotic freedom, the additional requirement (|4.68p yields a multi-valued /3 function. 
The interpretation of a{k) as a running coupling is then questionable since the zeros of the (3 
function are not necessarily fixed points jl53]rn Aside from the discussion in section E31 this 
indicates that the subcritical solutions are unphysical. 

Vertex corrections 

There exists the possibility of nonperturbative corrections to the running coupling. This might 
be interesting even for the extraction of Aqcd by comparison to experiment [140j . It was 
suggested in Ref. [98] to multiply the running coupling with a contraction of the ghost-gluon 

^^The infrared divergent running coupling l|4.30p from the heavy quark potential would not exhibit this 
infrared fixed point. 

^^The running coupling l|4.30p extracted from the Coulomb potential would yield an infrared finite limit for 
the subcritical solution and might serve as a more sensible definition in this case. 
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vertex in a certain momentum configuration. The infrared limit of the thus corrected running 
couphng is denoted by 



a 



(C) 



Hm a(k)C^ 



(4.69) 



The value C, defined by Eq. (|3.75|) . describes the dressing of the vertex in the infrared limit 
where first the gluon momentum is approached to zero and afterwards the ghost momentum 
is taken to be zero. With the help of the result obtained in table 13.21 for the number C, the 
corrections to Oc can be be specified. They are shown in table 14. l1 for the various solutions of 
K{d) in question, in comparison to the uncorrected values Oc- Note that the solution k = | for 
d = 3 is particular in the sense that it is not corrected by the vertex contribution since C = 1. 



d 




AT 

acNc 




3 


0.398 


4.488 


5.32 


3 


1 

2 


2tt 


2tt 


4 


0.595 


8.915 


10.94 



Table 4.1: Infrared Unite values of the running coupling for the various 
infrared solutions, with and without vertex corrections. For d = 3, we 
used the definition Uc in Eq. ij4.64)) whereas for d = A, Oc from Eq. ij4.66|) 
is used and reproduces the result in Ref. J97y without vertex correction. 



In concluding this chapter, we give a brief outlook on future investigations with the nonper- 
turbative running coupling in the Coulomb gauge. With an improved ultraviolet behavior of 
the Green functions, see section 14.31 it will be possible to adjust the scale in the Hamilto- 
nian approach to YM theory. This may be pursued by comparison to experimental scattering 
amplitudes from deep inelastic scattering. The nonperturbative feature of the variational so- 
lutions makes it then possible to quantify other interesting objects numerically, such as, e.g., 
the Coulomb string tension, or correlators of the high-temperature phase of QCD. 
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5 Inclusion of external charges 



In the Coulomb gauge, the Green functions were investigated in the infrared Hmit in chapter 
[2 and for the entire momentum range with variational methods in chapter El These vacuum 
calculations include only the interaction of gluons and ghosts, the coupling to the quark sector 
is suppressed by setting p^xti^) = 0- The translationally invariant Green functions so obtained 
were discussed in the previous chapters, and it was shown that the heavy quark potential is 
obtained by a perturbation in p^^^{x). Although the crucial feature of linear confinement is 
already given by the Green functions of the pure glue theory, the back-reaction of the external 
charges on the wave functional is lacking. One may expect that by extending the ansatz for 
the wave functional in the presence of a quark-antiquark pair, for instance, the energy can 
be lowered further by the formation of a flux tube connecting these external charges. In this 
section, various approaches are discussed by which the back-reaction of the external charges 
on the gauge field sector can be studied. 

5.1 Gluons in a quasi-particle representation 

The vacuum wave functional used in chapter [3] is of Gaussian type. In analogy to the quantum 
harmonic oscillator where the Gaussian ground state is annihilated by a linear combination of 
the coordinate and momentum operators, the variational YM vacuum state |0) suggests the 
introduction of an annihilation operator that yields a^(x) |0) = 0. Its hermitian conjugate 
a^(x)^ may then be considered the creation operator of a gluonic excitation at point x. By 
means of these gluonic excitations, an orthonormal basis of particle excitations can be con- 
structed |154| . In a normal order expansion of the Yang-Mills Hamiltonian, the term uja^a 
occurs so that uj{k) represents part of the energy of a gluonic excitation with momentum k. 
Since in the infrared, the dispersion relation departs substantially from the behavior uj{k) = k 
of a free particle, these particle excitations are regarded as quasi-particles, incorporating non- 
perturbative effects. 

Before we define the quasi-particle basis, recall that the choice A = ^ for the wave functional 



formally eliminates the Faddeev-Popov determinant J from expectation values. The Coulomb 
gauge scalar product ( | ) that comprises J after gauge fixing can be related to a scalar product 
( I with a fiat measure. 



in Eq. (HU 




(5.1) 




(5.2) 
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if we identify 

$[A] = (A|0) = v^^[A], d[A,U] = yrjlA]0[A,U]-^^. (5.3) 

With these manipulations, expectation values can be calculated straightforwardly by means 
of Wick's theorem. Note that the expectation values considered are not taken in the auxiliary 
Gaussian state ^[A] but actually in the state (|5.ip . The Gaussian "^[A] is suitable to define 
an annihilation operator 

atix) = ^jd^y {a^^ix, y)A^^{y) + iakm{x, yWmiv)) (5-4) 

that yieldfQ 

atix)\0)=0. (5.5) 
Above, the function aij{x,y) is defined as the inverse square root of the variational kernel cu, 

aij{x,y) = tij{x)a{x,y) , / d^ya{x,y)a{y, z) =uj'^{x,z) , (5.6) 



and has the same dimension as the field operator A. Adopting the canonical commutation 
relationJ§ 

[^^^(x),nf (y)]=it°,''(x,y), (5.7) 
the hermitian conjugate of Eq. (15. 4p obeys 

[<(x),a5(y)t] = t-''(x,2/), (5.8) 
and creates gluonic excitations by acting on the vacuum, 

n 

{x^'\ak, ik} ) = AA„ n <:(^^'V |0 > (5.9) 

k=l 

which can be shown to be orthogonal. Normalization of the states (|5.9p by an appropriate 
choice of J\fn niay be issued, if needed, in a finite volume V = f d^x. 



5.2 Coherent states 

With the set of gluonic excitations formulated in the previous section, we have a basis of the 
Hilbert space at our disposal in which a state can be constructed taking the presence of external 

^Strictly speaking, only the auxiliary vacuum |0 ) and not the true vacuum |0) is annihilated by the operator 
a. Employing the replacement In J = — J Ax A in the one- loop approximation, see Eq. I|3.13p . a definition of 
the operator a by means of that annihilates |0) could be pursued. 

^Here, we write the _L explicitly in order to emphasize that all fields are transverse. The commutation 
relations l|1.89p that involve the field dependent operator Tij are not considered here. 
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color charges into account. This state, is here chosen as a coherent superposition of the 
excited states (|5.9p . where a localized function controls the excitation of gluons in between the 
external color charges. This function, Z^{x), may be determined by the variational principle, 
minimizing the expectation value oi H' = H + Hext with the external charges pext contained 
in i^ext- 

Coherent states were first introduced by Schrodinger [155] in 1926 and further established by 
Glauber |156| in 1963 who laid the keystone for developments in quantum optics and thus 
gained the 2005 physics Nobel prize (along with Hall and Hantsch). It will be instructive here 
to briefly discuss the simple case of the three-dimensional harmonic oscillator and thus expose 
the main features of the coherent state. The Hamiltonian reads 



p2 1 

H = — + -muP'x^ = uja\ak + Eq , (5.10) 



where 

Ofe = + zcjpfc) and a\ = ^ - iapi^ (5.11) 

are linear combinations oixk andpfc that diagonalize H for cj~^ = ^Jrm2. Due to = iJjj, 

we have [oj, a^] = Jjj and the orthonormal states |n) = nfc("'fcO~^''^('^I)"'° |0) can be shown to 
yield the equidistant spectrum H \n) = En\n) with £"„ = ioiux + n^^ + + |). The vacuum 
state |0) for n = has minimal energy E'o = | and is annihilated by (|5.1ip . o |0) = 0. The 
solution to the latter differential equation yields the vacuum wave function in coordinate space, 

V^o(^) := (^|0) = -4^e-i-'/'^' . (5.12) 

We now define the coherent state 1^;) by acting with a unitary operator V on the vacuum state 
|0), 

\z) = U{z) |0) , U{z) = , zfc G C . (5.13) 

The coherent states so defined are normalized if |0) is, but they are not orthogonal. The 
continuous index z makes the set of coherent states over-complete in the given Hilbert space 
with a countable basis. Further properties of the coherent state are listed below: 

a) eigenstate of annihilation operator 

Ofe \z) = Zk \z) (5.14a) 



b) closure relatiorI§ 



1 = J- I cfz \z) (z\ (5.14b) 



''The integral is taken over the complex planes of all components of z. By inserting the closure relation into 
1 = (0| 0), the factor is cancelled by the integration of the polar angles in the complex planes. 
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c) minimal uncertainty 



d) Poissonian distribution 



{Axk){Apk) = ^ (no sum) (5.14c) 



|(n|z)|2 = 17^6-1^^1' (5.14d) 

k 

e) shifted vacuum wave function 

{x\z) {x - V2aRez\0) (5.14e) 

To see the property e) first note that the operator U{z) in Eq. (|5.13p can be expressed by 
coordinate and momentum operators x and p using Eq. (|5.11|) . 

having abbreviated 

Ik := V2a-Hm{zk) , Rk := V2aRe{zk) . (5.16) 
The well-known Campbell-Baker-Hausdorff formula 

In (e^ e^) = X + y + y] + ^[X, [X, Y]] - ^[Y, [X, Y]] + . . . (5.17) 

which infers e^"*"^ = e"^ e~2 '"^''^l for a c-number, helps to rewrite the coherent state 

(I5.13P and interpret it by a shifted vacuum wave function, 

tP.ix) = {x\ z) = e-^'^^^^ e'^^""^ e-^^^k i^oix) 

= e-l'^^^^ e'^^''^ i;o{x - R) . (5.18) 

Here, we have used {x\ ipk \y) = d^6^^\x — y) and assumed analyticity of il)o{x). That is, the 
wave function of the coherent state is, up to a x-dependent phase, the vacuum wave function, 
evaluated at a; — i?. This proves Eq. (|5.14el) . Equivalently, using {p\ ixk \q) = —d^6^^\p — q), 
we find for the momentum space wave function 

i^zip) = {p\ z) = e^'^"^^ Q-iRkPk _ /) , (5_i9) 

Expectation values in the coherent state can thus be easily related to vacuum expectation 
values, 



{z\0{x)\z) = j d^x^l{x - R)0{x)i)o{x - R) 

d^x i}l{x)0{x + R)ijQ{x) 
{Q\0{x + R)\Q) (5.20) 
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For an operator 0{x,p) that depends on both x and p, an equivalent relation is derived. The 
most general form shall be referred to as the shift rule: 

{z\0{x,p)\z) = (0| e-*^^^'=0(x + i?,p)e^^'=^'= |0) 

= {0\O{x + R,p + I)\0) (5.21) 

The utility of the above identities can be presented in a simple example. Consider the three- 
dimensional harmonic oscillator (|5.10p and construct coherent states as defined above, Eq. 
(|5.13l) . Now think of the particles as charged and an external electric field £ being switched 
on, H ^ H' = H — q£ ■ X. A prime in this chapter shall always indicate that external fields 
are present. Let a; = + x_\_ be split into a vector rry parallel to £ and one perpendicular to 
it (accordingly for R). The energy E'{z) in the coherent state is readily calculated using the 
shift rule (lOTTl . 

E'{z) = {z\H' \z)=Eq+uj\z\'^ -q£-{<d\x + R\Q) 

= Eq + -ua I + -—R^ + -—[R\\ 2 -o 2- (^-22) 

2 2 cj^ 2 cj^ \ " muj'^ J 2mu;^ 

Minimizing E'{z) clearly gives / = i?_L = and R\\ = The wave function '^0(2;) is 

therefore shifted to '4'zix) = ijjo{x — -^ji)- This happens to be the exact vacuum wave function 
of -ff', as seen directly by quadratic completion. 

The ground state energy (actually the entire spectrum) is shifted to Eq = Eq — 

In Yang-Mills theory, a less trivial situation is encountered. The Hamiltonian H is not di- 
agonalized exactly, but by means of the variational principle a quasi-particle basis can be set 
up. A coherent state constructed in this basis is used to calculate the expectation value of H' , 
and minimize the latter. We now discuss a comparable quantum mechanic example. Let the 
anharmonic oscillator be defined by 

i/^^h = — + -muj'^x'^ + Xm'^uj^x^ , A G M+ (5.24) 
2m 2 

may be a helpful quantum mechanic example. We also fail to diagonalize this Hamiltonian, 
but a Gaussian ansatz as in Eq. (|5.12l) with cr as a variational parameter aids to find an 
approximate ground state. Its width a can be shown to yield for 



(T^ = — + , C = 1350A2 + 30\/2025A4 - 3A2 - 1 (5.25) 

30Amu; V / 



the minimal energjQ 



3 3 15 
Eo = K + -mwV^ + — Am^cjV^ . (5.26) 
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Figure 5.1: Variational solution for the ground state of the anharmonic 
oscillator. Left: Energy Eq as a function of the parameter X of the quartic 
potential. Right: Width a of the Gaussian ground state, as a function 
ofX. 



As can be seen in Fig. 15.11 the Gaussian ground state wave functional becomes narrower for 
increasing the coefficient A of the quartic potential, while the energy increases. Note that for 
A ^ 0, the results for the harmonic oscillator are recovered. For a comparison to field theory, 
we discuss the variational solution to the c/)^ theory in the Hamiltonian approach in appendix 

El 

Similarly to the variational approach to YM theory, we have found an approximate Gaussian 
ground state and can create quanta by the creation operator in Eq. (15. lip with the Gaussian 
width a given by Eq. (|5.25p . Now switch on an external electric field, 

HU. = ^anh. -q£-x (5.27) 

and calculate the energy in the coherent state. The result gives the ground state energy Eq, 
see Eq. (I5.26p . and a correction dependent on z, 

E'{z) = {z\ H'^^^ \z) 

/2 1 

= Eq + + \m^u)^a'^R'^ + -muj'^(l + l^\moj(T^)R^ - q8 ■ R (5.28) 

2m 2 

For z = 0, the energy E'{z) is equal to Eq, as expected. Minimizing E'{z) w.r.t. z yields some 
value 2;™'°. The result will depend on the parameter A of the quartic potential and the electric 
force F = —q£. Without giving the result explicitly, we can see from Eq. (|5.28l) that for all 
values of A > and F / there exists a = RJ'"" ^ such that E'iRJ'"^) < Eq, i.e. the 
energy is lowered in the coherent state. 

We now turn back to the original task, writing down a coherent state in the quasi-particle basis 
of Yang-Mills theory. The over-complete set of states, here given by Eq. (|5.9|) . incorporates 

lowest-order expansion for small A agrees with Ref. [TSYj for the one-dimensional case. Particular care 
is required concerning the Riemann sheets when taking the A — > limit. 
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already the nontrivial effects of the Hamiltonian H (|1.107p without external charges, similar 
to the treatment of the quartic quantum mechanic potential above. Using these states, the 
coherent state for Yang-Mills theory is introduced here by 

\Z) = U\0) (5.29) 

with 

C7 = exp II d^x [Zl{x)al{x)^ - Z,"(x)*4(x)) | , (5.30) 

where the field Z^{x) controls the excitations of gluons. By means of Eq. (|5.4|] . one can rewrite 
U in terms of field operators, 

^7 = exp 11 d'x {iIt{,x)Al{x) - iRl{x)Iil[x))^ , (5.31) 

where the abbreviations 

i?j"(x) = V2 j d^x'ajk{x,x')Re{Z^{x')) (5.32a) 

//(x) = V2 j d^x'aJ^ix,x')lm{Z^{x')) (5.32b) 

were used. The expectation value of the Hamiltonian with charges, H' = H + Hg^t, in the state 
(|5.29l) will be calculated in the next section. Optimistically, we may hope that in analogy to 
the quantum mechanical example, a nontrivial solution for a coherent state exists that lowers 
the energy. 

We provide here the definition of the coherent state energy, 

E[Z] := (Z|i?|Z)flat = / (U^y h{u^) . (5.33) 

It is important to clarify which state it refers to as a Coulomb gauge expectation value of H 
(not H). Carefully treating operator ordering, the energy (|5.33|) can be rewritten as a Coulomb 
gauge expectation value. 



E[Z] = j VA'^*U^HU'^ 

- (^*-) (;^^'^) {:^"7i) {"^"Ti) 



/ 



VAJ {U^y H {U^) 
= {U-^\H\U-^) (5.34) 

It is seen that by calculating E[Z] via the auxiliary expectation value (|5.33|1 with a flat measure, 
the Coulomb gauge expectation value of the Hamiltonian H is taken in the state 

\Z) = [/ |0) , (5.35) 
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cf. Eq. (|5.29p . The corresponding wave functional is obtained by acting with U on the vacuum 
wave functional ^[^4]. With U given by Eq. (|5.30p . a Campbell-Baker-Hausdorff formula is 
needed to find U = J-^/'^U J^/"^ . Erom Eq. (f5T7ll . it can be derived that for [X, [X,Y]] and 
[y, [X, y]] c-numbers, 

Using this relation, and the one-loop Gaussian approximation for the Faddeev-Popov deter- 
minant (jS.lSp . we find after some algebra 

jj = Q-\SRxR+lAxRjj _ (5.37) 



5.3 Energy of Yang-Mills theory in the coherent state 

In this section, the calculation of the energy E'[Z] in the coherent state is presented. First of 
all, let us clarify how this energy is defined. The coherent state (|5.35|) is used to calculate the 
expectation value of the Yang-Mills Hamiltonian (|1.107|) in the presence of external charges 
/°ext(^)- Without external charges, the Hamiltonian reads 

H = Hk + Hp + hP (5.38) 

with the definitions of the kinetic part, 

H^ = ll 71:4] nfc(^) m^) (5.39) 

the magnetic potential part, 

Hp = \J d''xF^^{x)F^^{x) (5.40) 

and the Coulomb part He of the Hamiltonian which involves a contribution coming from 
dynamical gluonic charges, 

hP = d^l:,y] _L_ p^^J^) j[A] F^\x, y) 4„(y) . (5.41) 

In the presence of external charges, the total charge distribution p^{x) = p'^y^{x) + Pe^ti^) 
adds extra terms to the Coulomb part He, so that the Hamiltonian under consideration reads 

H' = H + H^(}^ + H^^^ (5.42) 

with 

H^c^ = y/ ^'Nl [^fdynWlA] F^''{x,y)pUy)+pt,,{x)F'^''{x,y)p''^yM) (5-43) 
and 

Hg^ = ^1 d''[xyKAx)F^''ix,y)pl,{y) . (5.44) 
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The energy in the presence of external charges in the coherent state |^) is defined by 

E'[Z] = {Z\H' \Z) , (5.45) 

where this expectation value involves after Coulomb gauge fixing the Faddeev-Popov determi- 
nant J. For the calculation of E'[Z], it will be extremely helpful to employ the shift rule as 
it was written down in Eq. (|5.21|) for the quantum mechanic oscillator. Here, we can find the 
equivalent "shifting effect" of the operator U, given in Eq. (|5.30|) . on the operator H' (the use 
of the tilde was explained in Eq. (15. 3p ). 



E'[Z] = {Z\H'[A,U]\Zh^.t 

= (0|^'[^ + i?,n + /]|0)flat (5.46) 

In the calculations below, we always use the transformations (|5.3|) . ending up with the fiat 
scalar product ( | For brevity, the notation 

(O)^ = ( OIOIO )flat , {d)z = (^lOI^)flat (5.47) 

is used. Furthermore, recall the matrix notation in coordinate space from chapter [H It will 
sometimes be convenient in this chapter as well. The spatial dimension d is left unspecified. 

The energy E'[Z] can be expressed by the results from chapter [3] if we set Z = 0, 

E'[0]=E + Vc. (5.48) 

Here, E is the vacuum energy displayed below Eq. (|3.19p and Vc is the Coulomb potential 
(|2.73|) in the vacuum state. In Eq. (|5.48|) . there are no contributions coming from 

2 

= y ( (^iQ) *Fpext + Pe.tF{iAQ) )^ 

= (5.49) 

We will see that in the coherent state the mixed term H^^^ does have a contribution to the 
energy. Let us emphasize that it is instructive to stick to the coordinate space representation of 
the Green functions. The transformation to momentum space is helpful in the vacuum where 
thus an overall infinite factor (27r)'^(5'^(0) can be extracted from the energy, due to translational 
invariance. The coherent state, on the other hand, accounts for localized external color charges 
and translational invariance is lost. Therefore, an overall volume factor cannot be expected 
and the Fourier transformation provides no simplification whatsoever. 

Let us turn to the calculation of the kinetic energy in the coherent state. The kinetic part Hk 
in Eq. (|5.39l) comprises the operator 

n^(x)[An] = n^(x)-l^^^ (5.50) 
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that may be understood as a function of A and 11. Applying the shift rule (|5.46p . we obtain 



Ek[Z] 



Hk[A,U[A,n]]) ^ = {Hk[A + R,U[A + R,n + I]] 
VA Il1{x)[A + R,Il + I]^[A]'^ 



(5.51) 



With the explicit form (|5.3p of the Gaussian vacuum wave functional '^[A] = (^|0), we use 

I[k[A + R,Ii]^[A] =iQk[A-R]^[A\ (5.52) 

to define 

l5\nJ[A + R] 



Ql{x)[A-R] = I d''x'uf^{x,x')A'^^{x') + ^- 



5Al{x) 



(5.53) 



with the obvious notation iofj{x,y) = 5"-^tij{x)uj{x,y). The kinetic energy in the coherent 
state thus yields schematically 



1 



Ek[Z] = -Ti[uj{AA)lo + uj{A 



5A 



+ 



1 / 5\nJ 5\nJ 

Ia Ja' 



d''xl\x) (5.54) 



which we calculate, following Ref. [83], with a quadratic completion. To this end, define the 
curvature Xz 



1 / 6^J[A + R] 
'2 \SA^{x)dA^^{y) ^ 



= +- Tr (^G[A + R]T';'''{x)G[A + R]T''.^''{y) 

With the definition of the ghost propagator Dq with a coherent background field, 

D^'^'{x,y) = {G[A])fix,y), 
the curvature (I5.55P can be expressed by 



ixz)^H^,y) = -TY Dirf"(x)Dir°'^(y) . 



(5.55) 



(5.56) 



(5.57) 



This implicitly establishes a definition of the proper ghost-gluon vertex rf '°(x) in the coherent 
state; the latter will in actual calculations be rendered tree-level, as before. Now, the second 
term in Eq. (I5.54P can be written using a partial integration as 



6\nJ[A + R] 
SA>^jiy) 



d''ziu;-'mx,z) 



6^J[A + R] \ 



5Al{z)5A]{y) / ^ 
d'z{u;-'mx,z){xz)tj{z,y) 



(5.58) 



In the one-loop approximation, one can realize by expansion of the Faddeev-Popov operators 
that the third term in Eq. (|5.54|) . (^^^^^^), is simply obtained by writing the first two 
terms as a modulo square. With Q = Q[A; R] as in Eq. (|5.53p . one thus finds 

[Qtix)Q',iy)) =1 [d'[uv][u;-Xzm^,u)^~Hn,v)[u;-Xz]iiv,y) (5.59) 
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and is ready to write down the result for the kinetic energy, 

Eu[Z] = \j d^lxyzWu; - Xzt,{^.y)^-\y.^)[^ - Xz]'ii^,^) + \j d'^xllixmx) . (5.60) 

By setting the coherent field to zero, Z = 0, the vacuum result (I3.20p is recovered. 

We now proceed to the Coulomb part He of the Hamiltonian where Eq. (|5.41l) amounts to the 
Coulomb interaction of gluonic charges, 

pSy^{x)[A,U] = Atix)Ulix)[A,U] (5.61) 

Taking the expectation value of such operators, relation (|5.52|) can be used. Noting that 
E^^^Z] = {H§\a,U])^ is a modulo square, the imaginary terms can be explicitly eliminated 
to find 

E^P[Z] = ^ ld'[xy]{{A + R)f{x)liix)F^^{x,y)[A + R]{A + R)f{y)lf{y))^ 

+ y / d''[xy]{iA + R)t\x)Qi{x)F'^^{x,y)[A + R]{A + R)fiy)Q%y))^ (5.62) 

with Q = Q[A; R] as given in Eq. (I5.53p . To 2-loop order in the energy, there are no contractions 
with the operator F. On the assumption that Dq is symmetric in color and coordinate space 
separately (just like Dg), we have 

{Qkix)[A;R])^ = ^TrD^rlix) = . (5.63) 
The non-vanishing contractions in Eq. (j5.62p then are 



(AQFAQ), (AQFAQ), RI{F)RI , (AIFAl) , (RQFRQ) 



A careful treatment of all indices giveq^ 



eP[Z] = -y y d''[xy]{f3)'"'Fl'''{x,y){f''f'' 

\u;-'y'{x,y) [ d''[nv][u;-xM^,u)uj-\n,v)[u;-x,]lj{v,y) 



+ { tfl{x)5\x,y) - j d'wu;-\x,w){x,)!^{w,y) 

tf;{x)5\x,y) - j d''wu:-\x,w){xj^{w,y) ] \ (5.64) 
+^ j d'[xy]Rt\x)li{x)Fr{x,y)R^'^{y)lf{y) 
I d''[xy]u;r.\x,y)tT {ii{x)Fz{x,y)i,{y) 

h'^lxyuvWu: - x,mx.u)u:-\u,v)[u: - x,]tjiv,y)Rfix)Fl%x,y)Rf{y) 



'The vacuum result lf3?22|l can be recovered for ^ = by using tr (r^f ") = -Nc{N^ - 1 
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5.3. Energy of Yang-Mills theory in the coherent state 



where we assigned 



Ft{x,y):={F'^\x,y)[A])^. (5.65) 



The expectation value of the dynamical charge Pdyn = ~ A • 11 in the coherent state is in view 
of Eq. (ran 

(pdV(x)[A,n]>^ = -Rt{x)li{x) . (5.66) 

Hence, the mixed terms ~ Pdyn-^/'ext in the Coulomb energy yield without contractions of the 
Coulomb operator F 

E^^\Z] = -g'l d%y]Rf{x)ll{x)F-,^{x,y)pl^,{y) . (5.67) 

The Coulomb energy of solely external charges obviously yields 

E^c\z] = ^1 d'[xy]p:^,{x)Ft{x,y)pt,{y) . (5.68) 

The magnetic potential energy in the coherent state is yet to be calculated. With the Hamilto- 
nian given by Eq. (|5.40l) . we find that the vacuum energy £'p[0] as given in Eq. (I3.2ip is found 
in addition to terms involving the coherent field, 

Ep[Z] = Ep[Q] 

+ 19"^ f ^"'^ f'^"^Ri{x)R]{x)K^{x)R'^^{x) 



+ g I d^xr'^ {diR%x)) R\{x)R]{x) 
+ \j d'^xRtix) {-d')Rnx) 



+ \g'N^-^^u;-H0,0) J d''x{Rt{x)f (5.69) 

All but the last term, call it £^p'^[Z], are quite trivially obtained using the shift rule (|5.47p . 
Since it is divergent, let us look at it more carefully. It comes from the quartic part of Hp after 
applying the shift rule and contracting two gauge fields A, leaving two coherent fields R, 

E^'''[Z] = j d^x{(yl,(x)4-(x)/?i(x)^,(x)>^ + (i?i(x)i,(x)i,(x)i?,(x)>^ 

+ {Ai{x)Rj{x)Ai{x)Rj{x)) ^ + {Ri{x)Aj{x)R^{x)Aj{x)) ^ 
+ {A,{x)Aj{x)Ri{x)Rj{x))^ + {Ri{x)Rj{x)Ai{x)Aj{x))^^ (5.70) 
Using the symmetry of the vacuum gluon propagator, we get 

Ef-'IZ] = \g^ j d^x {2(A,(x)4(x)i?,(x)^,(x)>^-2(A,(x)i2,(x)i?,(x)^,(x))^} 

= ^g^ jd'^x[ {tij{x)u;'Hx,y))\y^JrR,ix)Riix) 

- {tu{x)LO-\x,y))\y^JvRj{x)Rj{x)] . (5.71) 
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Using a Fourier transform for uj ^(x, y) and tr (f"f^) = -A^c'5"^ we recover the proposed 
expression in the last hne of Eq. (|5.69|) : 



d^x I [^j af'pU,{p)u;~\p)j R';{x)R\{x) - d^ptii{p)u;-\p)j R''j{x)R]{x) 
= I d^pu-Hp) I d'^x Rtix) (^5,,id - 1) - <5., ^) Rfix) 



]g'N,^-^-^u-H0,0) I d''x{Rt{x)f , (5.72) 



4 

where the identity 

d''pf{p')UM = 6.,^ I d'pfip^) (5.73) 

has been employed. We note here, however, that ti;~^(0, 0) = / d'^puj~^{p) is not necessarily a 
well-defined quantity. Actually, if we consider that the gluon will essentially behave like a free 
particle for high momentum, uj{p) p, there is an ultraviolet divergence of degree d — 1, i.e. 
quadratic for d = 3. We therefore have an undefined expression in the coherent state magnetic 
potential (|5.69l) which calls for renormalization. 



Renormalization 

The renormalization of the vacuum Green functions as discussed in chapter [3l suffices to render 
the shift AE'[Z] := E'[Z] — E'[0] of the coherent state energy from the vacuum energy finite, 
as will be shown now. The ultraviolet divergence in the magnetic potential energy (I5.69P comes 
as no surprise, for the magnetic potential operator (I5.40p is a highly local object. We now try 
to single out the divergent object from the magnetic energy (|5.69|) . To this end, we first realize 
that it is contained in the quadratic term in the (terminating) Taylor series of Ep[Z] about 
Z = 0, 

Ep[Z] = {Hp[A + R])^ 

Such a term with a second derivative inside a Gaussian expectation value can be rewritten as, 
see Eq. (IE.6|) in the appendix, 

2i 



= - [ d^[xy]Rtix)( ^ \ R\{y) 



-2 I d^[xyuv]Rnx)u;Z{x, u) ^^^f^ u;%{v, y)R]{y) (5.75) 
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5.3. Energy of Yang-Mills theory in the coherent state 



Note that E^^ ^ [Z] is not given by only in Eq. (|5.72p but also includes the other term 

in Eq. (I5.69|l that is quadratic in R. 

In view of the gap equation, we know that uj was chosen to minimize the vacuum energy, 
5{H)i^/6uj = 0. Therefore, the divergent expression 5{Hp)^/5ijj in Eq. (I5.75P can be expressed 
in terms of the kinetic vacuum energy {Hk)uj and the Coulomb vacuum energy {H^^)^^, 

Ef^[Z] = +2 / d'[xyuv]Rt{x)u;t:^{x,u) ^^f\^^''^^- u:%{v,y)R]{y) . (5.76) 



The above expression is now finite if we replace the bare form factors by the renormalized ones. 
To arrive at the finite expression for Ef\z\, one needs to evaluate 



iu;-')t:i,ix,u){u;-')i){v,y) (5.77a) 



i/) = ^{i^-'mx,u) j d^z{u:~y^{y,z)xil{z,v) (5.77b) 
5 (Q'L(u)Qi{v)) 1 

^ ' ~dit,,-l\aa'i^ i\,.a'cil „.\ / jd „ ,U , —\\hh' i „ , „J\,h'di„i 



-- J dV(a;-^)^? (x,:r')x?^(^',u) J {u:-%,{y,y')x%{y' ,v) (5.77c) 

The implicit w dependence of the form factors x(^) and of the vacuum Coulomb propagator, 

F^\x,y) := {F[A]r^ {x,y) , (5.78) 

is omitted in the computation of Eq. (I5.76P : it would bring about two-loop terms in the 
equations of motion. With the relations (|5.77p . Eq. (I5.76P can be shown to give 



Ef\Z] = ^Jd''[xyz]Rf{x){u;{x,yMy,z)-xix,y)x{y,z))Rnz) 
+ ^iV,y" d''[xy]FUx,y) 

d'^[uv]Rt{x)[uj{x, u) - xix, u)]tij {u)uj-\u, v) [uj{v, y) - x(w, y)]R'-{y) 
+ {tij{x)uj-^{x,y)) J d'^[uv]R^{u)[uj{u,x)u;{y,v) - xiu,x)xiy,v)]R]{v) 
-2 J d''uRnx)x{y,u)R'^{u)tji{y) (^6''{y,x)- j d''vu;-\y,v)x{v,x) 

'(5.79) 

The magnetic potential energy in the coherent state thus finally yields 

E,[Z]-EM = j d''xr'r''R'i{x)R]{x)Ri{x)R%x) 

+g j d^xr^" {m'^ix)) R\{x)R''j{x) + E^f\Z] (5.80) 
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where Ep[0] is given by Eq. (l3:2T]l and Ep '[Z] by Eq. (ISTTflll . The shift AEp = Ep[Z]-Ep[0] is 
indeed finite for the variational solution of the coherent field (any divergence will be suppressed, 
for Z = is an option) . 



Solution for the imaginary part of the coherent field 

The energy E'[Z] in the coherent state sums up the magnetic energy (|5.80|) . the kinetic energy 
(|5.60|) . and the Coulomb energy contributions (|5.65p . (I5.67P and (|5.68|) . The real and imaginary 
parts of the coherent field Z^{x) are to be determined by the variational principle, 

6E'[Z] _ SE'[Z]_ 

Since the Yang-Mills Hamiltonian has cubic and quartic terms in the fields, we cannot alge- 
braically solve these equations. However, the energy expression is quadratic in the imaginary 
part of Z. The relevant terms can be written in matrix notation as 

^« = l/M/ - bl . (5.82) 

This defines the matrix M as well as b by 

Mt^ix,y) = tf^{x,y) 

-g'Rfix)Fl'^ix,y)Rfiy) - ^{f'^y'u;;^\x,y)Fz'%x,y){f'Y' (5.83) 
blix) = j d''ypl^,{y)Ff{y,x)Rt{x) (5.84) 

The quadratic form (|5.82p attains its minimal value E^^^ = —^bA^^b for 

il{x) = j d%{M-Xi^,y)b'j{y) ■ (5.85) 

In order to write down the coherent state energy E'[Z] minimal for the imaginary part of 
the coherent field, replace the fields I^{x) by the expression (|5.85l) . This yields the complete 
energy expression 

E'[Z] = ^|d^[xyz][a;-xj?f(x,y)c.-i(y,z)[a;-xj$?(z,x) 
-y / d''[xy]{faf'Fl'''ix,y)if''f'' 

(w'^)J(a;,y) j d'^[uv][u; - Xz]'ltix,u)uJ-\u,v)[uj - Xz]kjiv,y) 
+ (tf/(x)5'^(x,y) - I d''wu;-\x,w){xz)!fiw,y) 
t^^{x)5''{x,y)- [ d''wu;-\x,w){xz)'ji^,y) 
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4 

-.2 



d'^[xyuv][uj - Xz]tk{x,u)uj Hu,v)[uj - Xz]kj{v,y) (Ri{x)Fz{x,y)Rj{y) 
+ / / d^[xyuv]pl^,{x)Ft{x,u)R\'^{u) {M-%' iu,v)R^%)Fi%v,y)pl^dy) 



ab 



1 

1 



+\g' I d''xr^r''m{x)R]{x)Rt{x)R%x) 



+ g j d'^xr"' [m'^jix)) R\{x)R'^{x) 

+ \ j d'^[xyz]R^{x) {uj{x, y)uj{y, z) - x(x, y)x{y, z)) R^{z) 
+ ^N,J d''[xy]FUx,y) 

d'^[uv]R^{x)[uj{x,u) - x{x,u)]tij{u)uj-^{u,v)[uj{v,y) - x{v,y)]Rj{y) 
+ {Uj{x)uj^\x, y)) j d'^[uv]Rt{u) [uj{u, x)uj{y, v) - xiu, x)x{y, v)]Rjiv) 
-2 / d''uRt{x)x{y,u)R'^iu)tjiiy) (6''{y,x)- [ d^v vo~\y,v)x{v,x) 



(5.86) 



with the abbreviation (M {u,v) given as the inverse of the matrix M in Eq. (|5.83|) 



Estimate for translationally invariant form factors 



The expression (|5.86|) for the energy in the coherent state is quite large and it might be 
useful to make crude approximations to get an idea about its main characteristics. In Eq. 



(I5.86p . the kernel uj{x,y) is fixed by the variational calculation in the absence of external 
charges, see chapter O The other Green functions, (Dq)"-^ {x , y) , F^{x,y) and (Xz)ij(^)y)) 
are dependent on the location of the charges and thus have lost their translational invariance. 
Moreover, neither the Lorentz nor the color structure can be expected to be trivial. Therefore, 
a minimization of the energy E'[Z] is a very costly calculation, not even vaguely comparable 
to the challenge in solving the integral equations in the vacuum. 

A possible simplification is achieved by rendering the Green functions insensitive to the pres- 
ence of quarks. Setting Dq = Dq, Xz — X and Fz = F^^, we can simplify the expression (|5.86ll 
for E'[Z] quite significantly, 

E'[Z] = Ek[0]+EP[0]+E^^\0]+Ep[R] 

+ N,^ j d''[xyuv]FUx,y)Rtix)[u; - x]ix,u)u;T^\u,v)[io - x]iv,y)R]iy) 

+ / / d'^[xyH/>Lt(^)^c.(x,n)iif (n)(M-i)^;(tx,t;)/?f • (5-87) 
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The imaginary part of Z is already chosen so that E'\Z\ is minimal, this gives the last term 
quadratic in F^. The term in the second line of Eq. (|5.87p is exactly cancelled by the term 
in Ey\B\ that came from the expectation value (RQFRQ)^, so that the energy shift AE' = 
E'[Z] - E'[0] yields 



AE'[Z] = ^g"" j d''xr'r''Rt{^)R]{x)Rt{x)Rj{x) 
+ g j d''xr'''{m-^{x))R\{x)R'^{x) 



1 ' .d. 



+ 2 j d''[xyz]R^{x){uj{x,y)uj{y,z)-x{x,y)x{y,z))Rt{z) 
+ —Nc / d'^[xyuv]F^{x,y)uj-^{x,y)Rl{u)[uj{u,x)uj{y,v) - xiu,x)xiy,v)]R'^{v) 
-N, I d'^[xyuv]F^{x,y)Rtix)x{y,u)R'^{u)tji{y) 
5^{y,x)- j d'^vuj~^{y,v)x{v,x) 



2 



+<7^y d'*[xyHy9^,t(^)i^c.(^,^^)^rW(M-^)^|(n,t;)i?f • (5-88) 

The prominent terms are the very first and the very last ones, resembling the energy of the 
anharmonic oscillator (|5.28p in an external electric fieldlfl Approximating M « 1, 

AE'[Z] = d''xR';{x)Rt{x)Ri^{x)R^^{x) 

+/ / d''[xyu]pl^,{x)F^{x,u)Rf{u)Rf{u)F^{u,y)pl^^{y) , (5.89) 



and furthermore replacing the matrix RkRk by its mean diagonal elements, 

Rk{x)Rk{x) ^ t ^''^k{x)Rk{x) ^ _^_N^^^^^ ^ cPix) := Rt{x)Rt{x) , (5.90) 

we get 

^E'[Z] = j d'xcfix) - j d'^xc^ix) (I d%FUx,y)pt^dy)y ■ (5-91) 

The expression (I5.9ip attains for (j){x) = 2g'^(^J d'^y F^{x,y)p'^^^{y))^ its minimal value 

AE'[Z] = J d'x d'yFUx,y)pUy)y^ ■ (5-92) 

By dimensional analysis we find in the infrared limit of quark separation r, noting that 
Fiv ~ a^r'^^~^'^~'^ with Coulomb string tension a^, we get in general dimensions AE'[Z] ~ 



_^4^16k+ii 4d gj^^ hence for d = 3 and k 



1 



AE'[Z] ~ -ay , (5.93) 



^Studying the infrared limit, the third, fourth and fifth terms can be neglected since Lu{k) = for fc ^ 0. 
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a completely senseless result. First of all, confinement is lost since Eq. (|5.93|1 overwhelms the 
confining potential for large r. But even worse, the energy E'[Z] is negative which contradicts 
the positive definiteness of the YM Hamiltonian. 

What has been done wrong? With various ansatze, such as for SU (2) 

Kix) = r'^^dm^ix) , (5.94) 

the energy in the form (|5.87l) was minimized and the erroneous large r behavior (|5.93l) was 
confirmed. Therefore, the approximations done after Eq. (|5.87p cannot be blamed, the awkward 
result (|5.93p must come from the use of vacuum Green functions. In a sensible approach, the 
Green functions must take into account the presence of external charges. 

The proper way to go about this has to involve solving a set of Dyson-Schwinger equations 
without translational invariance. One of these equations follows from the variational principle, 
6E'[Z]/6Z^{x) = 0. The ghost propagator Dq can be determined by the DSE 

{D^r'=G-'[R]-J:z (5.95) 

where the self-energy T,z is obtained by using the propagator Dq in the definition of S in Eq. 
(|2.3ip . A relation for the inverse ghost propagator, such as Eq. ()5.95|1 . might not be so useful. 
Instead, the equation 

D^ = Dg + J (fx (g[A] tY{x)RI{x) G[A + R]^^ (5.96) 

seems more applicable for solving the set of integral equations. 

It was not attempted to solve the DSEs in the variational coherent state. Immense computer 
power would be necessary and such a project is at best problematical. In the upcoming section, 
the effect of the coherent field on the confining property of the heavy quark potential will be 
discussed with alternative methods. 



5.4 The gluon chain 

Having learned from first approximations that the Green function need to be sensitive to the 
presence of external charges in order to arrive at a sensible result, we try to make estimates for 
such Green functions in this section. The quest for the correct wave functional has been long- 
standing. After Dirac's description of the abelian case |158| . the generalization to YM theories 
proved to be difficult |159| . There have been various attempts of formulating flux tube models 
(see e.g. Ref. |160| ). calculations in the Hamiltonian approach to QCD [161j . and also lattice 
calculations |162j that elaborate on the subject of the ground state in the presence of heavy 
quarks. One of the most intuitive notions might be the gluon chain model [163111641 1165| . and 
we will come back to it below. 

An advantage of the variational approach is that any ansatz for the wave functional will do. 
A too wild guess will be rejected if it leads to a rise in the energy. If the ansatz is somewhat 
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sensible, however, and a decrease in the energy can be detected, this wave functional is to be 
favored over the vacuum wave functional. A description of the back reaction of the presence of 
external charges on the vacuum wave functional is then found. With this reasoning, we here try 
to find a coherent field Z^{x) that leads to a lowering of the infrared Coulomb propagator. At 

(2) 

very high quark separation, the Coulomb energy [Z] can be expected to give a dominant 
contribution. We will now focus on the calculation of the respective expectation value in the 
coherent field. 

Note that the Coulomb propagator Fz corresponds to the vacuum counterpart F^^ with the 
coherent field i? as a background field, 



Fz = {F[A + R])^ = I VAMA]\^ew J <fxRl{x) 



F[A] . (5.97) 



5Al{x) 

We find a series in the coherent field where each term is treated using the identity 

^G[A] _ (^[^]r°''*(x)G[^] . (5.98) 



5Al{x 



Acting on the Coulomb operator F = GGq G, the variational derivative generates a sum of 
terms, 

dF[A] _ ^r.i-nO.a. .1 , .ipO,a 



SAtix) 



F[^]rr (x)G[A] + G[A]Tl'''ix)F[A] . (5.99) 



One can check that the combinatorial factors cancel the 1/n! from the expansion of the expo- 
nential function in Eq. (|5.97p . and one obtain^ll 

Fz = {F + FMG + GMF + FMGMG + GMFMG + GMGMF + ...)^ 
= ((1 + GM + GMGM + . . . ) F (1 + MG + MGMG + ...))^ 
« Y^Y.^DcMy^FUMDcy- (5.100) 

m n 

with the abbreviation 

Tab I 



M'^\x,y) = j d^zRl{z){Tl\z)Y\x,y) , (5.101) 

where the tree-level ghost-gluon vertex r^'"(x) was defined in Eq. (|2.20l) . In the last line of 
Eq. (jS.lOOp . the approximation of factorizing the Coulomb expectation value was employed 
(for a discussion, see section IXHll . We have thus expressed the Coulomb propagator Fz in the 
coherent field by vacuum expectation values, i.e. the vacuum ghost propagator Dq and the 
vacuum Coulomb propagator F^^ (|5.78p . and the coherent field itself via M in Eq. (|5.10ip . 

The Coulomb potential Vq (r) in the coherent field is defined for a given charge distribution 

Pextix) by 

VSir) := y / d'[xy]pt,,ix)Ffix,y)pl,iy) , (5.102) 



''A proof is given by showing {Rkj^)^G = n\{GM)"G by induction, which is trivial, and plugging it into 
G[A + R]^J2n 7diRkj^)"G[A] which leads with F^GGo^G directly to the second line in Eq. l[ETOO|l . 
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+ 



Figure 5.2: The Coulomb propagator in the coherent Geld, denoted by 
the double line on the l.h.s., and its expansion Ii5.100\) . The solid lines 
represent vacuum Coulomb propagators, the dashed lines vacuum ghost 
propagators, the dots are tree-level ghost-gluon vertices, and the wiggly 
lines indicate that gluons are excited by the coherent held. 

cf. the vacuum Coulomb potential Vc in Eq. (I2.73p . With two point-like heavy charges located 
at the origin and at r, 

pt,,{x)=6'^^6'{x)-6'{x-r)) , (5.103) 

one may think of the expectation value in Eq. (|5.102|) as shown in Fig. 15.21 As the charges 
are separated by a large distance r, it can be energetically more favorable to excite a gluon 
in between the charges. This notion basically corresponds to the gluon chain model proposed 
in Refs. [1631 11641 1165| . These gluons possibly screen the potential and thus account for a 
lowering of the Coulomb string tension ac to the string tension aw from the gauge invariant 
Wilson loop. Moreover, with excited gluons along the flux tube, the thickness of the string 
(lacking for the one-gluon exchange) may be generated. The coherent state is suitable to check 
whether this mechanism is realized in the framework of the variational approach. 

The guiding properties for choosing an ansatz for the coherent field are transversality and 
symmetry. Firstly, let us exploit transversality. The coherent field Z^{x), here the real part 
R'j^{x), is always projected transversally, i.e. it occurs only in the form tjk{x)R1{x). Without 
loss of generality, we can therefore set 

Rl{x) = e^jdih'^ix) (5.104) 

and determine the field h1{x) instead. 

Secondly, we can assume that the symmetry of the charge distribution is mirrored in the field 
R'j^{x). While in the absence of charges, the Green functions are translationally invariant, 
a charge distribution Pext(^) remove this invariance, yet leaving us with some weaker 
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symmetry. Assume for instance, one of the charges is located at the origin of our coordinate 
system, and the second charge is smeared on the surface of a concentric sphere with radius rjf] 



pl,,{x) = 5^H\x) - 5^'^ 6{\x\ - r) (5.105) 



1 

The exchange energy is the same as for the charge distribution Eq. (I5.103|) . and the choice 
(|5.105p may therefore equally be used for the purpose of these investigations. We refer to 
the setting (|5.105p as the sphere picture and Eq. (|5.1Q3p as the string picture. For the charge 
contribution of the sphere picture (|5.105p . the coherent field Z'^{x) can be expected to be point 
symmetric, 

Z^'pa;) = DZ^'ix) , (5.106) 

that is, for any rotation D G S0{3) about an axis that contains the origin a; = 0, it behaves 
according to Eq. (|5.106|) . Expanding Z"(a;) in spherical coordinates, this means that the 
coefficients merely depend on the distance \x\ from the origin. The field h'*(x) defined in 
Eq. (|5.104p inherits this property from Z^{x) via Eq. (|5.32p and suggests an expansion into 
spherical coordinates. On the other hand, if we use the string picture, see Eq. (I5.103|) . only 
a cylindrical symmetry may be expected, i.e. the rotations D are restricted to those about 
the symmetry axis along fq. It is then useful to decompose h^{x) in the cylindrical basis 
= {ep{x),e^{x),e^} with || fq, 

h^ix) = h^pip, z)ep{x) + h%{p, z)e^{x) + h^ip, z)e, . (5.107) 



Now consider the gluon chain integral (15. 100(1 rl Plugging in the definition of the field /i^(x) 
into the matrices M, any integral involving h'j:{x) is of the fornix 

[DGMDG]{z,y) = geijk J d'^xdfDG{z,x)hj{x)d^DG{x,y) = 1 (5.108) 

The above object corresponds to one "chain link" and the entire gluon chain is a convolution 
of chain links we intend to re-sum. 

We may choose the first integrations as the one where the charge sitting at the origin is attached 
to the ghost line. This corresponds to setting jz = in Eq. (I5.108p . The distance r between 
the charges being considered very large, the infrared asymptotic forms of the vacuum ghost 
propagator may be used. From the momentum space expression (see table [27111 . we find the 
dimensionally regularized (d = 3 + 2e) coordinate space ghost propagator 

VG{x,y)- J (Tk^^^^^^^e "(4^)^/2 r(l + k) V 4 ) ' ^^"^"^^ 



The bold type notation for vectors is reinstated for the remainder of this chapter. 

transformation to momentum space would be, on the one hand, convenient since the infrared asymptotic 
vacuum Green functions from chapter [2] can be used. However, a subtlety is accompanied by the Fourier 
transformation. The real part of the coherent field in coordinate space, ^(Z{x) + Z{x)*), transfers to ^iZ{k) + 
Z(—k)*) and only by knowledge of the parity properties of Z{x) can we tell whether the above momentum 
space quantity is real, purely imaginary, or complex. 

^"Recall that a caret denotes an object in the adjoint representation. Here: hj{x) — T°'h'j{x). 
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and taking the derivative as in Eq. ()5.108p gives for k = and d = 3, 

B Xk-yk 



dtDG{x,y) = ^dtm^ 



2tt^ {x — y 



|2 



(5.110) 



Plugging this into the chain hnk (15.10811 . we get 



B^ I'n \Y(x)-(x>iy) 
[DaMDa]{^.y) = -g— J d'x l,[J^_J' (5.111) 

The vector h(a;) is expanded in the cyhndrical basis QSa;, see Eq. (|5.107|1 . and we expand 
X = pxep{x) + ZxGz and y = py cos ipep{x) — pysm{pe^{x) + ZyBz accordingly. The outer 
product in Eq. (|5.11ip then reads 



XX y 



^ PyZx sin If 

PyZx cos If - PxZy I . (5.112) 

Y -pxPy sin if ^ 



The (/9-integration eliminates both the p- and the z-component of h"(a;) in Eq. (|5.11ip . To 
prove this, note that both hk = hk{p, z) \/k and the product x ■ y = p^Py cosip + z^Zy that 
occurs in the denominator are even functions in ip. The components hp and contribute with 
the outer product (I5.112|l a factor of simp, odd in 93, and hence must vanish. We are left with 
h^{p,z), the only function that will contribute to the corrections to the Coulomb propagator. 

In order to estimate the effect of the infrared (|a;| 00) behavior of h"(a;), we switch to the 
sphere picture. This is achieved by letting y || which infers py = 0, py = y, and cos 93 = 1. 
On the assumption that a non-vanishing h^^ component can be realized in the sphere picture, 
use the ansatz 

h^ix) ~ ^ . (5.113) 

By virtue of the homogeneity of the functions in the chain link integral (|5.11ip , one can show 
that for a/j < 1, the infrared strength of the Coulomb propagation is weakened by each chain 
link in the series in Fig. 15.21 This means that for r — > 00, the only linearly rising term would 
be the very first one of the expansion, the vacuum Coulomb propagator F^^. On the other 
hand, if a/i > 1, the infrared strength is enhanced by each chain link integral, meaning that 
for some minimal number of chain links, the integral would be divergent. Such a choice will 
be avoided by the variational principle. If (and only if) a/^ = 1, each chain link produces a 
constant factor and all terms of the expansion contribute to a linearly rising potential. The 
numerical value of this factor can be computed by means of the two-point integral (|B.26p in 
the appendix which applies to Eq. (|5.11ip by 

d'x , ^' , = 2^^ ■ (5.114) 

|ic|-*(a; — y)^ 

The string picture will require a more involved calculation but the statement aimed at here is 
only on a qualitative level. It shall only be noted that with the choice a/i = 1, the chain link 
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integrals i 

[DgMDgKKO, y) = LDg{0, y)dk (5.115) 

yield a calculable constant L, here not given explicitly. It is only important to remark that L 
is a matrix made of real numbers@ Integrating the gluon chain (|5.100p "from left to right", 
every chain link produces a factor of L in the infrared limit, 

VS{r ^ oo) 

= y j d^[xy]pt^,{x) (^l + L + LL + ...y'F^%x,y) + L + LL + . . .y'pLiv) 

(5.116) 

The excitation of a single gluon at first sight seems to lower the Coulomb interaction if L < 
(which depends on the self-consistent solution for Z^(x)). 

So far, the color structure has been suppressed. Recall that L = L"'T°' is a color matrix in 
the adjoint representation. The Coulomb energy Vq (r) in Eq. (I5.116|) is sensitive only to the 
diagonal element i^J^. The diagonal elements vanish for L. Therefore, the excitation of a 
single gluon does not contribute to the Coulomb energy. The gluons have to come in pairs. If 
we turn now to SU{2), two chain links are seen to produce a factor 




+ 



-2{{Ly + {L'f) 



which is manifestly negative. These contributions therefore lower the Coulomb interaction, 
due to the color structure. We will now proceed to re-sum the entire series (|5.116l) for SU [2). 

There are two geometric series in L occurring in the gluon chain (I5.116p . The value of the 
series can be computed by noting that with = Xla-^"-^"' 

(LL)"^ = LU''-rt2 , (5.117) 
which leads to the following result for odd and even powers of L, 

^l^yn+i ^ (_l2)"X , (L)'"+' = (-L2)"ll . (5.118) 
Hence, a geometric series in the matrix L can be expressed by a geometric series in the number 

oo oo 

2n+l , / f\2n+2 



c = i + l + ll + ... = i + Y,{lY'^ +Y.i^ 
1 ^ 1 

iTL^^+rTL 



n=0 n=0 

1 + -1—L + -l—lL . (5.119) 



^^The coordinate space ghost propagator can only be defined as a distribution. Therefore, we act with a 
derivative from the right on it which is always present in the gluon chain (see Fig. 15. 2p . 
^^It is not as straightforward to keep track of the complex phase in momentum space. 
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The series converges if 

< 1 , (5.120) 
and the elements of C are found in view of Eq. (|5.117p |^ 

= f d'^* + L"L* + L'^*) . (5.121) 



1 + L^ 

The Coulomb potential in the coherent field V^{r)^ see Eq. (|5.116p . depends on the matrix 
element (CC)^^. Defining a string tension az by the infrared limit of the Coulomb potential 
Vq (r) in the coherent field, 

VSir ^oo) = (CCf^acr =: azr , (5.122) 

we get from Eq. (|5.12ip thai0 

1 - L2 + (3 + L2) (L3)2 

fJz = ^ ^ — -Tc ■ 5.123 

(1 + L2)2 ^ ^ ^ 

The variables and are entangled by the convergence condition (|5.120p . With spherical 
coordinates in color space, ^ := |L| and =: |L| cos?], Eq. (I5.123P can be expressed in terms 
of independent variables < ^ < 1 and < r/ < vr. 



[1 + ^2-^2 



az=[l- sin^ r] \ \^^^ ] ac ■ (5.124) 
It is thus seen that 

< az < (Tc , (5.125) 

i.e. the Coulomb string tension is indeed lowered by the excitation of gluons. This very notion 
was proposed in the gluon chain model and it follows here from an ansatz for the wave functional 
of SU{2) YM theory with heavy quarks. The upper bound, az = ac is reached only if 77 = 0, 
i.e. |L| = L^. 

Depending on the coherent field Z^{x), the string tension az could in principle be lowered 
just above zero. In this case, L has no component along the color vector pext(a^), i-e. r/ = ^, 
and its length is maximal, i.e. |L| /" 1. Thus, we find az \ 0. However, one may expect that 
this leads to a rise in the magnetic potential energy. Adding more and more gluons into the 
Coulomb interaction must increase this part of the energy since it is ultralocal in the coherent 
field, see Eq. (|5.80l) . At a given mean number of excited gluons (the peak of the Poisson 
distribution), a saturation of the energy must occur. Realizing this effect computationally 
may be difficult, though. If the magnetic potential energy indeed saturated the number of 
gluons for large separations r, it would have to also rise linearly in r, otherwise it would be 
negligible in the asymptotic infrared. 



^■^This matrix is the inverse to 1 — L, as one may expect. 

^■^This result is sensitive to the approximations made above. At this stage, only qualitative statements can 
be made. 
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5.5 Squeezed states 

Although the estimate of the Coulomb propagator in the coherent state in the previous section 
indicated a lowering of the string tension, it is not clear whether a full self-consistent solution 
will exhibit the same feature. A possible alternative description of the back reaction of the 
gluonic sector on the external charges may be motivated, again, from the quantum mechanic 
harmonic oscillator. In three dimensions, the charged oscillator can be studied in an external 
uniform magnetic field B = iJey which couples to the angular momentum. To this end, the 
Hamiltonian ()5.10|) is supplemented by the substitution 

Pk^Pk-qAk, A = -^£CxB. (5.126) 

The new Hamiltonian H' = H + Hp + Hpi contains the Larmor frequency '^L = linearly 
in the paramagnetic term 

Hp=ulL\\, (5.127) 
and quadratically in the diamagnetic term 

Hd = ^mivlx'i . (5.128) 

The coupling of the angular momentum operator L = x x p to the external magnetic field B in 
Hp (|5.127p is reminiscent of the coupling of the dynamical charge p^^^{x) of YM theory to the 
external charge Pg^^ (x) in the Coulomb interaction PdynFpext ■ This similarity is acknowledged 
by noting that p"^^^ = — A"^ • 11* is an outer product of coordinate and momentum operators, 
albeit in color space, just like the angular momentum L. For the calculation of the energy in 
the coherent state, this has important consequences, as will be seen. 

Before we focus on the coherent state, let us give the exact ground state solution of the rather 
simple Hamiltonian H'. The primed quantities always include the effect of the external field. 
Thus, we can interpret 

^ / 2 2 I -2 2\ , 1 / /2 2 



H + Hd = -uj[ cj^pJ + a-^x\ + -J {a%\ + a'-^xi) (5.129) 



as one harmonic oscillator of frequency oj (and width cr^ = ;^) in the ey -direction, and two 
oscillators perpendicular to it with frequency oj' = (w^ -|- w|)^''^ (and width a'^ = ;^^). Using 
the annihilation operators 

ay = icjpii) , ak = ^{(T'~'^Xk + i(T'pk) , k = x,y, (5.130) 

where aj^ = OxGx -|- o-y^y, the Hamiltonian H' can be re-expressed, but the paramagnetic part 
Hp is still not diagonal. By means of a Bogoliubov transformation, 

a± = {ax ± iay) , (5.131) 
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the complete Hamiltonian is diagonal, 

H' = uj ^aluz ^ 2) ^ (^+"+ ~'~ + 1) ~ (^a\_aj^ — a^_a- + 1^ , (5.132) 

in the orthogonal basis of the states {a\.)"'+ (oL)"— (al)"'^ |0). The exact ground state of H' 
therefore has the energy 



E'q = -uj + +ljI> -uj = Eo, (5.133) 

i.e. it is higher than without the external field (w^ = 0). The interesting part is that the exact 
ground state iPq{x) = ip^ ||(a;)^g j_{x) is Gaussian, 

(TTfJ^) ' 

^oA^) = -^;T7^e-i</-", (5.135) 

but the width a' in the directions perpendicular to the external field B is changed, as compared 
to the width a of the Gaussian in the direction parallel to B. By switching on the external 
field B, the ground state wave function is squeezed by a factor smaller than 1, 



V 1 + ^ 



in perpendicular directions. 



The effect of squeezing cannot be expected from a coherent state which merely shifts the 
vacuum wave functional. We now construct a coherent state in the eigenbasis of H, and 
determine the parameters by minimizing E'{z) := (z| H' \z). This is equivalent to the 
approach in YM theory, see section [531 The minimal value E'{z) can then be compared to E'q 
in Eq. Kim . 

Define the coherent state by 



|0) , a = ((j-^a; + fcjp) (5.137) 



and calculate the coherent state energy E'{z) = E{z) + Epiz) + E£){z), 

E{z) = cu^(Imz)2 + (Rez)2 + (5.138a) 
Ep{z) = 2u;L(Rez_L X Imz_L) -ey (5.138b) 
Ed{2.) = ^mwl ((T2 + 2o-2(Rez_L)^) (5.138c) 

Note that the coupling of the angular momentum to the magnetic field brings about a term 
Ep{z) linear in Imz^. This also occurs in the energy expression (|5.86p of YM theory where 
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the couphng of pdyn to pext is responsible for such a term. By quadratic completion, it is easily 
found that for Imz = ^Rez x ey the energy E'{z) is minimal. In this case, 

= ^uj + l"^ +Lo{Rezf . (5.139) 

Minimizing this expression w.r.t. Rez, one finds that Rez = and thus Imz = 0. Therefore, 
the trivial solution |z) = |0) is energetically most favorable. In the first line of Eq. (|5.139p 
one might hope for a non-trivial solution due to the minus sign. Such a situation occurs also 
in YM theory. However, this negative term is exactly cancelled by the diamagnetic energy 
()5.138cp with cj^ = Since the coherent state with Z = still gives a higher energy than 
the exact ground state energy (|5.133|1 . 

E'{0) -E'o = u;(^l + ^^-^Jl + ^^ >0, (5.140) 

it must be inferred that the coherent state cannot mimic the squeezing of a Gaussian wave 
function. 

Of course, the comparison of the coherent state to the quantum mechanical harmonic oscillator 
and the Yang-Mills theory is by no means supposed to indicate a complete analogy. For 
instance, quartic and cubic terms may allow for a nontrivial shifted ground state — the coherent 
state — to be energetically more favorable. Also, the coupling of the "outer product" pdyn to 
the external field in YM theory is mediated by the Coulomb operator F[A] which is dependent 
on the field operator A itself. Such a situation is not encountered for the harmonic oscillator 
and we have seen in the previous section that the expectation value of F[A] in the coherent 
state may alone provide part of the desired effect. Despite these differences of the oscillator 
coupled to a magnetic field on the one hand and the gluonic vacuum with external charges on 
the other, it might be a useful idea to consider squeezed states as an alternative to coherent 
states. This is a possible future investigation. An introduction to squeezed states can be found 
in the literature on quantum optics [166| . 
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Summary and outlook 

The gauge principle, which is the main ingredient in the very definition of the gauge field sector, 
recurs in all aspects of Yang-Mills theory. Starting with the quantization procedure, discussed 
in chapter [H gauge invariance requires a careful treatment which is particularly difficult for 
nonabelian gauge groups. The generator of (time-independent) gauge transformations is the 
Gauss law operator. In the canonical quantization approach, it is not possible to promote the 
Gauss law to an operator identity. We have shown that by means of fixing the gauge completely 
on the classical level and imposing quantization conditions within the Dirac bracket formalism, 
the Gauss law does hold on the quantum level. The gauge-fixed Hamiltonian operator thus 
derived determines an unambiguous time evolution of the quantum system. It agrees with 
the Christ-Lee Hamiltonian [2^ where a projection on gauge invariant states is employed. 
This agreement is notably reassuring in the sense that there is no ambiguity in the choice 
of the quantization procedure. For Yang-Mills theory in the temporal Coulomb gauge, it 
makes no difference whether Dirac quantization {first quantize, then constrain) or constrained 
quantization (first constrain, then quantize) is used. In the light of the gravitational force, the 
method of quantization is to date a controversial issue. 

In the temporal Coulomb gauge-fixed configuration space of Yang-Mills theory, the uniqueness 
of the gauge fixing condition demands the restriction to a compact region, the fundamental 
modular region. In the Gribov-Zwanziger scenario, this infers the enhancement of boundary 
effects in that region. We have investigated in chapter [2] the consequences for the Green func- 
tions in the temporal Coulomb gauge in d spatial dimensions of (d-l- l)-dimensional Minkowski 
spacetime. The Landau gauge results for (d + l)-dimensional Euclidean spacetime can be 
obtained by shifting d ^ d + 1. In the temporal Coulomb gauge, the stochastic vacuum 
{"^[A] = 1) is sufficient to account for the infrared asymptotics of the theory. The horizon 
condition, which enhances the ghost propagator in the infrared, is along with the nonrenor- 
malization of the ghost-gluon vertex the most important feature for the infrared behavior of 
propagators in the temporal Coulomb or the Landau gauge. Using power law ansatze for these 
propagators, the set of Dyson-Schwinger equations was solved with a tree-level ghost-gluon 
vertex. A single infrared exponent k can be extracted to specify the infrared power laws of both 
the ghost and gluon propagators, for a given dimension d. In the 3 -|- 1 dimensional Coulomb 
gauge, the result k = ^ is shown to yield a heavy quark potential which rises linearly. Quark 
confinement is thus, qualitatively speaking, a consequence of the horizon condition. It was 
argued that this result, is not necessarily unique. First of all, there is one further solution for 
d = 3, K « 0.398, which leads to a Coulomb potential Vc{r) that rises less than linearly |83| . 
With the Coulomb potential being an upper bound to the gauge invariant quark potential 
Vw{r), the latter solution would indicate that Vwir) also rises less than linearly, contradicting 
the lattice results. In this sense, it is more likely that the solution k = | is realized. Moreover, 
the off-shell gauge condition introduces an additional parameter C into the infrared analysis, 
on which the solution for k generally depends, due to the approximation of the ghost-gluon 
vertex. It was found that all solutions but the one for d = 3 and k = ^ depend on the value of 
(. This infers that the latter solution is less sensitive to the approximations made. Its corre- 
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spending continuous branch of solutions Ka{d) = ^ — 1 as a function of the dimension d exists 
for all d, even for very high dimensions where the other solutions cease to exist (see Fig. 12.61) . 
The infrared power laws were furthermore extended by an ansatz that incorporates powers of 
logarithms in the infrared. It was found that values for these exponents exist for which the set 
of integral equations is also solved in the asymptotic infrared. The power law exponents, how- 
ever, remain unchanged. Hence, there is a degeneracy in the asymptotic power law solutions, 
also for d = 3 and k = |. Qualitatively, the logarithmic corrections do not change the power 
law behavior and thus also not the confining property of the Coulomb potential. 

A full solution for the propagators of Yang-Mills theory in the temporal Coulomb gauge 
is expected to show both, the infrared behavior that infers confinement and the ultraviolet 
behavior that agrees with perturbation theory. The variational solution to the Yang-Mills 
Schrodinger equation, discussed in chapter [3l provides such a solution approximatively. The 
Gaussian types of wave functionals were shown to reproduce exactly the infrared behavior 
(k = |) anticipated in the previous chapter. Thus, numerical results with a linearly confining 
potential were obtained [SB]. While the propagators are insensitive to the details of the wave 
functional, it was found that the three-gluon vertex is extremely dependent on it. For one 
wave functional considered, it is equally zero, whereas for another it shows a strong infrared 
enhancement. However, in the Dyson-Schwinger equations of the propagators the dressing of 
the three-gluon vertex was found to have no effect in the infrared |103| . The ghost-gluon vertex 
was studied in the infrared gluon limit from its corresponding Dyson-Schwinger equation. It 
was shown that generally the infrared limits of the vertex' momenta are not interchangeable, 
except for the solution d = 3 and k = ^. Again, the regularity of the vertex indicates that 
this solution may be favored. Chapter [3] is furthermore concerned with the assessment of the 
approximation of the Coulomb form factor. It was shown that the horizon condition needs to be 
relaxed in order to arrive at a solution that satisfies the Coulomb form factor DSE. The infrared 
analysis of the ghost DSE and the Coulomb form factor DSE proves that a simultaneous 
solution with infrared divergent form factors is impossible. The so-called "subcritical solution" 
[120] where the form factors are infrared finite solves all DSEs but does not provide the expected 
infrared behavior of the quark potential. If a critical solution exists, then the Coulomb form 
factor DSE, as it stands, will be violated. Higher-order effects must then be included. An 
approach that includes higher-order effects may be to calculate the gap equation in the form 
of Eq. (IE.14P up to two loops (avoiding a three-loop expression for the energy density). This 
may be tedious but possible. 

The ultraviolet tails of the variational solutions for the Coulomb gauge Green functions have the 
correct power law behavior but the anomalous dimensions, related to the powers of logarithms, 
turn out incorrect. This might be an effect of the simple Gaussian shape of the wave functional. 
In chapter m the ultraviolet gluon energy cj(A;) ~ kln^^^^ k and ghost form factor d{k) ~ 
l/ln^/^^/c were motivated. This followed from a nonperturbative definition of the running 
coupling by the ghost-gluon vertex. In a calculation carried out in both Landau and Coulomb 
gauge, the ghost DSE was analyzed in the asymptotic ultraviolet using appropriate ansatze. 
From the requirement that the nonperturbative running coupling coincide with the ultraviolet 
limit of perturbation theory, it can be deduced what the anomalous dimensions are to leading 
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order. The way to reproduce these anomalous dimensions must be to extend the Gaussian 
wave functional. At the very least, the gluon loop must be included. A non-trivial dressing of 
the three-gluon vertex, similar to Landau gauge studies, is most probably necessary. This is a 
possible future investigation. 

Moreover, chapter |4] dealt with the infrared limit of the running coupling which can be calcu- 
lated analytically. In all gauges that interpolate between the Landau and the Coulomb gauge, 
this value is the same, but changes discontinuously in the Coulomb gauge limit |103| . Vertex 
corrections to the infrared limit of the running coupling were also calculated, showing that the 
only solution that has no vertex correction is the one where d = 3 and k = ^. 

The final chapter [5] uses a quasi-particle representation of gluonic excitations to incorporate 
the back reaction of the presence of external charges on the gluon sector. To this end, coherent 
states are motivated from quantum mechanical examples. The calculation of the Yang-Mills 
energy in the coherent state produces a vast expression that depends on a localized background 
field. It was possible to remove divergences from this expression by usage of the gap equation. 
Even after the explicit variational solution for the imaginary part of the background field, the 
energy expression requires some approximation to be processed further. A first approximation, 
rendering the Green function insensitive to the presence of external charges, yielded senseless 
results. Subsequently, the Coulomb potential was estimated in the asymptotic infrared by 
an ad hoc ansatz for the background field. It was possible to expand, calculate, and resum 
this expectation value and indeed show that the string tension can be lowered by such an 
ansatz. The determination of the exact factor by which the string tension is lowered, requires 
further investigation. Another quantum mechanical example was put forward to motivate that 
squeezed states may serve equally well for the incorporation of external charges. 
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A Conventions and notation 

A.l Units, metric and group conventions 

Throughout this thesis, natural units are used, 

h = c=l. (A.l) 

Since he 200MeVfm-\ Eq. (lAHl infers that Ifm corresponds to 200 MeV. The Minkowski 
metric is here defined by the metric tensor 
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(A.2) 



Antihermitian generators of group transformation in SU (Nc) are normahzed as follows. For 
the fundamental representation, we choose 



tr (^T"r^j = -^6"'' (A.3) 
and for the adjoint representation 



It follows that 

Nr 



^^arpbrpc^ jade jheg jcgd c jabc 



A.2 Notation 



In momentum space integration, the phase space factor 27r is absorbed by the definition 

in analogy to h. In the case of multiple integrations, the abbreviation 

d[xyz] = dx dy dz (A. 7) 

is frequently used. 

Matrices in coordinate and intrinsic space, for instance G, have the components G"'^{x, y). The 
functional trace "Tr" is then defined by 

TrG = j d'^[xy]6'^{x,y)5^''G''\x,y) (A.8) 
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Functional determinants can be related to the functional traces (lA.Sjl by the identity Det G = 
exp Tr In G. 

The Fourier transformation of a function 0{x), operating on spacetime, defines by 

0{x) = [ d'^pdip)e'P-'= (A.9) 



the function 0{p) operating on momentum space. Although these are different functions, the 
tilde is often omitted. For a function of two variables, we may define a Fourier transform by 

0{x,y) = J d^q] e'P-''d{p,q)e-'''-y . (A.IO) 

In translationally invariant systems, 0{x,y) = 0{x — y), and 

d{p,q) = [ d'^[xy] e~'P-''0{x,y)e''^-y 

ddy e%-p)-?// [ d'^xe-'P^'^-y^ 0{x-y) 



= {2T,f5''{p-q)d{p), (A.ll) 

such that the propagator-like object 0{x—y) is directly related to 0{p) via Eq. (|A.9j) . Whereas 
in coordinate space, a matrix in coordinate and intrinsic space is denoted by a single symbol, 
for instance G, the Fourier transformed object is usually a function of a single momentum 
scale and has no color structure, cf. Eq. (|2.11l) . 

Higher n-point functions are defined equivalently. A vertex function with translational invari- 

ance, 0{x^ y, z) = 0{x — y,x — z), leads to 

d{k,q,p) = j d%yz] e-'''-''0{x,y,z)e-'i ye-'P-' 

= i2nf 5'^ {p + q + k)d{k,q) (A.12) 

and momentum is conserved by the delta function. Choosing the momentum routing sets the 
signs in the exponent of the definition of the Fourier transform. 

Functional derivatives can be shown to yield 

^ df^p e-'P-'' . (A.13) 



SO{x) J " SOip) 
Note that the sign in the exponent is opposite of that in Eq. (IA.9h 
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B Two-point integrals 

Here we sketch the evaluation of nonperturbative Euchdean loop integrals in asymptotic limits. 
According to the discussion on page[72], we can replace the integrands by their asymptotic forms 
and encounter the two-point integrals 

'^^ (^2)a((^ _ fc)2)/3 ^ «,/3eK,rnEN. (B.l) 

These can be shown to be homogeneous functions of the momentum k. By a scaling of the 
integration variable, £ \i, one readily finds that since the two-point integral can only 
depend on the scale, it should obey /3) ~ {k'^Y the exponent of the power law can 

be determined to he k = d/2 — a — j3 + m. After applying the usual trick of introducing 
Feynman parameters, 

1 1 

1 r r r^a-i 13-1 1 

' dx / ciy^(x + y-l)—^^— (B.2) 



CfCf J J ' ' {xCi+yC2r+P B{a,f3)' 

^ ^ 

where B(a,P) is the Euler beta function, we can shift the integration variable, i ^ i — yk. 
The integrand then depends on i'^ only and we can integrate it out. For m = 0, 1, 2, 3 we need 
the following standard integrals |2]: 

d'i ^ 1 r(n-d/2) fiy-^/' 

(£2 + A)" (47r)'^/2 r(n) ^ ' 

a'nd, ^ 1, 1 nn-d/2-i) (ly-'i^-' 

(£2 + A)" 2 ^^(47r)'^/2 r(n) \/\) ^ ' 

Integrals with an odd number of vectors i in the numerator vanish by symmetry. For our 
purposes, we have A = xyk'^. The two-point integrals can be straightforwardly computed 
using the identity 

1 1 

Jdxjdydix + y- l)x--V-' = B{a, (3) = (B-5) 



for the Euler beta function. One then finds the results 

1 r(d/2 - a)r((i/2 - /3)r(a + /? - (i/2) 2n 



-0 



^ , _ 1 r(d/2 - g + l)r(d/2 - 0)T{a + /3 - d/2) 2W/2-a-g+i .u^ 

^^^"'^^ - (4^0^ na)nmd-a-(5+i) ^^-^^^ 

„ 1 r(d/2 - g + 2)r(d/2 - P)T{a + (3- d/2) 2.d/2-a-0+2 

-2{a,p) T{a)T{P)T{d-a- P + 2) ^ ' 

1 1 r(ci/2 - g + l)r(rf/2 - /? + l)r(a + P-d/2-l) ..d/2-a-0+2 

^ 2 {47ry/^ T{a)T{p)T{d - g - /? + 2) ^ ' ^ ' 

^ . .^ 1 r(d/2 - g + 3)r(d/2 - /?)r(g + /3 - d/2) 2^d/2-a-ff+3 

"3l«,PJ T{a)T{P)T{d-a- p + 3) ^ ' 

3 1 r(d/2 - g + 2)r(rf/2 -P + l)r(a + p-d/2-l) 2^d/2-«-g+3 m 

^ 2 (4^)'i/2 T{a)T{p)T{d - g - /? + 2) ^ ' ^ ' 
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B. Two-point integrals 



The above formulae are valid for those values of a, /3, m only for which the integrals converge. 
At £ = A:, a pole is integrable as long as (3 < d/2. On the other hand, the infrared convergence 
at £ = depends on m: 

{d/2 + m/2 for even m 

d/2 + m/2 + l/2 for odd m ' ^ 

One can relate this inequality to the requirement that the arguments of the first gamma 
functions in the numerators of Eqs. (IB.6P be positive. The ultraviolet convergence is contained 
in the third gamma function of the numerators. For convergence, the relation 

I d/2 + m/2 for even m ^ 

y d/2 + m/2-l/2 for odd m ^ 

has to be satisfied. Obviously, odd values of m, compared to even values, work "in favor" of 
convergence both in the infrared and in the ultraviolet, due to the angular integration. In 
sums of UV divergent two-point integrals, the angular integration eliminates the divergence in 
some cases. For instance, the Brown-Pennington projection of the gluon propagator in Landau 
gauge [163, 

Ri=\k)D^,{k)^\nK (B.9) 

eliminates the unphysical quadratic divergences proportional to the metric tensor, by virtue 
of Rfi^{k) in Eq. (|2.82p . We demonstrate that this elimination can be due to the angular 
integration for the Brown-Pennington projection of the off-shell ghost loop (I2.81al) for d = 3 
and the ultraviolet ghost propagator (I3.58cp . 



2Nc 



i-ai.kY + ic-i)- 

A 1 



d{i)die - k) 



di — , / dx 



167r2 J y/^F) 



-1 



1 - C^' + (C - 1) ^ 



1 



{£ - kf 
kx 



{P + A;2 - 2ekx)y/ln{P + k^ -2ekx) 



- ^ fdi-^(^-(3-C) + 0(^ 
= (3 - C) O (i^) + finite (B.IO) 
For C = d = 3, the linear divergence of the ghost loop vanishes. 

Another cancellation of UV divergences is encountered when implementing the horizon condi- 
tion in the ghost DSE (I2.49p . The proposition claimed in section [2^ is that although the ghost 
self-energy S is UV divergent, the naive usage of the formulae (IB.6P yields the correct result 
for the finite difference d^^{k) — (i"^(0). The proof is here restricted to two-point integrals 
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B. Two-point integrals 



I{k) = Ho(a, P) with n := a+f5 such that I{k) does not converge in the UV. Let ^ > n > f — 1- 
With the Feynman parameter trick (|B.2j) . we get 



1 1 



Hk) 







ii-k)y 



1 1 



j dx j dy 6{x + y - l)x''~^y^-^Jix,y;k) 



(BAD 







where = 2tt'^^'^ /T{d/2) and the function 



J{x,y;k) = - I df 



2j (£2 + A)" 





(B.12) 



with A = xyk"^ is regulated by the UV cut-off A. Employ the mean- value theorem, 



1 



1 



nA 



to realize that J can be written as J = Jreg + Joo with 



—1^ o<e<i, 



(B.13) 



Jo 



i 1 (A2)f-^-^oo 



nA 



Jreg(A:) = I df 



2^2"! A^oo 



2 7"" (£2 + ^A)"+i 





< CXD 



(B.14) 
(B.15) 



for the values 2>''^>2~^'^^ exponent n = a + p. In order to determine J^eg = J — Jo, 
explicitly, write 



Jieg{k) 



A2 





A2 



(£2 + /^)n (^2)n 



'2^/)2\|-l 



d 1 1 



r(n) 



(it 



n-l 



g-(^2+A)t_g-f2i 



(B.16) 



With help of the identity 



(B.17) 
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B. Two-point integrals 



we can now compute Jreg: 



>reg- 

A2 





1 oo 



1 r(f)r(n-f + 1) 1 d 



^ ^ 



Id r(f)r(n-f) 

Plugging this result back into the definition of I{k) in Eq. (IB.llj) . one finds 



I{k) = Ioo+Ireg{k) (B.19) 
1 



(4vr)tr(f)(f-n) 



^oo = — ttt^tt:^ -(A')^~" (B.20) 

— n) 

^ (27r)'^r(f) S(a,/3)r(n) 2 4 '2 ^ 



1 T{d/2 - a)T{d/2 - P)T{a + p- d/2) „^2^d/2-a-/3 



(47r)'^/2 T{a)T{p)Tid - a - (3) 



The regular part of I{k) would have been found correctly by direct usage of Eq. (|B.6aj) . despite 
the UV divergence. 

Finally, another result of a two-point sort of integral is provided that comes into play in the 
gluon chain in chapter 15.41 It reads 

X(a,,3) := /.i^^ ■ (B.22) 

where in d dimensions \£ x k\ can be understood as the projection of i to the hypersurface 
perpendicular to k. The Feynman trick (IB.2j) yields after the shift i ^ i — yk (which does not 
affect the numerator in Eq. (IB.22j) ) 

1 1 

^("' /3) = jdxjdy6{x + y- 1) x--^"' j dH j^^^^, ■ (B.23) 



If Q is the enclosed angle of t and A;, |£ x A;| = iksinO, its integration gives 

TV 

I desm^-'d = B{^,^) , (B.24) 
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B. Two-point integrals 



whereas the remnant angular integrals simple yield surface of the unit sphere with d ^ d — 1, 
Thus, Eq. (lR23]l is calculated by 

1 00 







1 r(f + i-a)r(| + i-/?)r(f)r(a + /?-f 

(47r)'^/2 r(a)r(/3)r(f -i)r(d+l-a-/3) 

(B.26) 
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C. Infrared limit of the ghost-gluon vertex 



C Infrared limit of the ghost-gluon vertex 

The rather technical project of determining the infrared hmit of the ghost-gluon vertex to yield 
the solution shown in Table [3?2l is presented here. We use the truncated vertex DSE as shown 
in Fig. 13.61 To illustrate how to keep track of color group and vertex factors, first note that the 
triangle diagrams comprise a trace of generators, tr ^f'^'f^'fc^ ^ g^g y^^n g^g three vertex factors 
of ig, i.e. 




{igY tr ( f^f^f" 



^ 2 



(C.l) 



due to Eq. (jA.Sp . We now consider the infrared gluon limit of the first loop diagram (|3.79p . 
With the identity 

^2 {I . qf 



(■ai3{q)ta0 - q) 



(C.2) 



we can express this diagram in the infrared limit of the remnant ghost momentum q by two- 
point integrals, 



limrfG^)(0;g,,) = r(0)(g)C5^^AS^(H,(l + 2K,^-2. 



■3(2 + 2k, ^ - 2k) /q' 



rf(,)ciii, 



where 



d-l 



(47r)'^/2 d (1 + 2 k) T{d/2) 



(C.3) 



(C.4) 



and we have used Eq. (|2.60|) . The factor of C appears in Eq. (|C.3p since the infrared limit 
is taken as in Eq. (|3.75p . A factor of | stemming from the color trace (|C.ip is left explicit. 
Considering the second loop integral of the ghost-gluon vertex DSE, the result (|3.7ip for the 
three-gluon vertex in the infrared limit is plugged into the expression (|3.78p to find that 



liinrP^)(0;g,9) = {-2)ig' A' B^hq^^ J d^e i ■ q (l - (p ■ if 



1 L 



i-2)igq^ ^ ^ {Ei{d/2 -k,1 + k)- E^{d/2 - k + 1, 1 + K)/q^) 



Tfiq) (-2)^ 



1 hh 



4 I\ 



(C.5) 



143 



C. Infrared limit of the ghost-gluon vertex 



where I3 was defined in Eq. (|3.72|) and 

1 d-l 



(47r)'^/2 {d - 2K)T{d/2 + 1) ■ 



(C.6) 



The factor of (—2) in Eq. (IC.5I1 is the prefactor of the ghost triangle in the three-gluon vertex 
DSE, see Fig. 13.41 Altogether, the infrared gluon limit of the ghost-gluon vertex, defined by 
Eq. (137751) yields, 

and one can find C to obey 



C = 2^4^id-l)d^T{l + -)-Kri-)jTid-3K)n--^n-KyT{- + K)T{l + 2K) 
r(2 + 2k) - y^r(^ - 2K)3r(l + ^ + KfTi2 - ^ + 3k)|^ 
I (d-l) {d-2K) T{d - 3 K)T{-K)r{l + 2 Kf 

(^4i+«r(i + ^)r(^ - K)r(-K)r(^ + k) + 0Fr(^ - 2K)r(i + ^ + ^)) } 



(C. 



As can be checked, this leads to the numerical values of C given by table 13.21 for the various 
solutions of K. A cancellation of the two loop diagrams, see Eq. (I3.8ip . must yield C = 1 in 
Eq. (|C.7p . In this case 

^^hiG^ id -i){d- 2^)r((i - 3K)r(i - K)r(2K)r(2K + 2) 
hh r2(f -2K)r(f + K + i)r(2 + 3^v-f) 



and among the solution Ka{d) and that solve Iq = Ia, see Fig. 12.5^ Eq. IC.9I uniquely 

yields k = ^ with d = 3. 

One might object that the Bose symmetry of the ghost triangle contribution to the three-gluon 
vertex is broken, for only the upper vertex is dressed while the others are at tree-level (see Fig. 
13.41) . In the derivation of the DSE, the channel to start the calculation in is an arbitrary choice. 
Therefore, one might as well find the upper vertex dressed. If we repeat the calculation, the 
DSE for the ghost-gluon vertex now reads 



(C.IO) 
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C. Infrared limit of the ghost-gluon vertex 



cf. Eq. (|C.7p . and the solution can be found to be 

C = A{d-l){2K + lfT^(^ + l)T{d-'iK)T{^-K + l)T\-K)T\2K + l)/ 

|d(2K + i)r(^)r(K + i)T\'^- + K + i)r(3K + 2 - ^)r'{^ - 2k) 

+{d -l){d- 2k)T{^ + l)T{d - 3k)T{-k)T{2k + 1)T{2k + 2) 

(r(^ - 2K)r(K + i)r(^ + + 1) + 2r(^ + i)r(^ - ^)r(-K)r(2K + 2)) | 

(C.ll) 

Again, the limit d 3 and k — > | yields the result C = 1. Apparently, it makes no difference 
which one of the vertices is at tree-level. 
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D. HAMILTONIAN approach to (p^ THEORY 



D Hamiltonian approach to 0^ theory 

The (/)^ theory is a popular toy model to illustrate some aspect of a physical field theory in a 
simpler context. Unfortunately, the (j)^ theory in d = 4 dimensions suffers from the triviality 
problem which states that after renormalization, the coupling parameter is equally zero and a 
free field theory is considered |168j . This problem can be circumvented by choosing a negative 
coupling parameter which is at the prize of losing the boundedness of the spectrum. In this 
appendix, we will use the (j)'^ theory in order to show how the variational principle can be 
applied, in accordance to chapter [3l 

With the classical Lagrangian density 

C = \id^cl>)d^<P-^m^<p''-^<P'' (D.l) 



and the conjugate momenta 

TT = 

6d, 



^ = #1 = 5V (D.2) 



we can define the classical Hamiltonian 



H = ld'x (ivr^ + \ (V0)^ + imV^ + ^0^) . (D.3) 

The equal-time canonical quantization relations, [(p{x,t),'ir{y,t)] = i6^{x,y), lead to the repre- 
sentation 

idcpix) 



of the momentum operator, where x is a 3-vector and the time is fixed in the Schrodinger 
picture. Thus, the quantum theory is defined by the Hamiltonian (|D.3p in terms of the field 
operators. 

The time-independent Schrodinger equation H {ip) = E lip), along with the semi-definiteness 
of for A > infers that the vacuum wave functional iplcp] = {(j)\ip) can be determined by a 
variational principle 

Eo = (V'l H \%b) — > min. (D.5) 
As an ansatz, we choose a Gaussian wave functional 

= AAexp^-^y" cfx j (fy(l){x)uj{x,y)(j){y)^ = Af exp ^-^(j)^uJxy(l)y^ (D.6) 

where the short-hand notation keeps the coordinate dependence as an index and the integration 
is implicit by means of the sum convention. The kernel oj is determined by the variational 
principle (|D.5|) . 
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D. HAMILTONIAN approach to (f)^ THEORY 



Expectation values are computed by 

where the generating functional Z\j] is defined as 

= j V(j) eXp {-(i)x^^^y4)y + j. 



o 



Z[j] 



(D.7) 



j=0 



AA^ j exp I • exp {^jxUJ^yjy 



exp I -j^a;^.j^jy 



(Di 



Above, after a quadratic completion, the expression in the curly brackets is (1) = 1 by virtue 
of the normalization of the wave functional via N . 

The vacuum energy Eq is now calculated term by term, as a functional of w. First of all, the 
mass term yields 



e4 



Ju^uv jv 



m? S 6 
ui 7t' e 



xx' 



j=0 



The above term is basically the trace of the propagator 

{4>{x)<f>{y)) = l-uj~'^{x,y) . 



(D.9) 



(D.IO) 



In order to understand better the meaning of the quantity uj^^, use the explicit notation and 
relate it to its Fourier transform, 



d^[xy]io~\x,y)6^{x,y) = J d^[xy] J d^koj-\k) e'^'^^'-y^ S^{x,y) 

V [ d^kuj-^{k) 



(D.ll) 



where V = j d^x is the volume of coordinate space. The first variational derivative in the 
calculation (ID.9P brings down a source term jx> from the exponential so that the second 
variational derivative acts not only on the exponential but also on the latter source term jx'- 
Setting sources to zero at the end of the calculation generally eliminates some terms. Some 
combinatorics is necessary to deal with the quartic term in the energy, we just write down the 
result: 



^ ^ d'^x{{(i){x)f) = ^v I d^eLo-\e)uj-\i-k) 



4! 



(D.12) 
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D. HAMILTONIAN approach to (p^ THEORY 



The kinetic energy is found to give 



V(t) 



V / d^kuj{k) 



where the propagator (ID.lOp was used. 

Furthermore, the spatial derivative terms of the Hamiltonian (ID.3P yield explicitly 



d^xU{x)(-dl)cp{x)) 



d'[xy]6\x,y){-dl){ct>{y)4>{x)) 
d\xy]5\x,y){-dl) f d'k^u^~\k) e'^-^'^-y^ 



J d^kuj-\k)k^ 



Summing up all contributions, the energy is found to be 

E = V Jd^k (^^uj-\k)ik^ + m^) + ^ujik) + ^ j dH uj-\k)uj-^ {£ - k)^ 

and it is minimized w.r.t. the kernel uj, 

5E/V 1 







which leads to the "gap equation": 



(D.13) 



(D.14) 



(D.15) 



(D.16) 



(D.17) 



The term "gap equation" is sensible since even if the mass parameter m was zero to begin 
with, there would be a mass dynamically generated by the last term in Eq. (|D.17jl . Hence, it 
describes the lowest allowed energy level of the spectrum, i.e. a mass gap. That term stems 
from the quartic term in the Hamiltonian and is also present in the YM Hamiltonian. It is a 
divergent constant subtracted in our approach. In the (/>^ theory, Eq. (ID.17jl must be solved 
iteratively and renormalization becomes necessary. For further reading, see e.g. Ref. [169| . 
If such a term is of any significance for YM theory, then it will have an effect only on the 
intermediate regime of oj{k). The infrared dominant curvature term from gauge fixing is of 
course absent in the (f)'^ version (ID.lTjl of the gap equation. 
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E Identities for Gaussian expectation values 

The Gaussian expectation values 

{0)^ = j VA^*[A]0^[A] (E.l) 
with the normalized wave functional, 



can be shown to yield some useful identities. Take a variational derivative of Eq. (IE.2P w.r.t. 

ujf^{x,y) =tij{x)ujix,y)6''\ 

to find 

. (E.4) 



For expectation values (lE.lh of operators O = 0[A], a variational derivative thus gives 

^^^^^ = {At{x)A^{y))^{0)^ - {A^{x)A]{y)0)^ , (E.5) 



i.e. a covariance. Using Wick's theorem, a few manipulations yield the identity 

\ - -4 / d'^[uv]u:^^{x,u) J^f)- ^y ^Av,y) . (E.6) 



For O = Hp[A], this leads directly to Eq. ([535]). 

A little more effort is required for expectation values of operators O = H[A,Jl\ that involve 
the momentum operator 11. The equivalent of Eq. (IE.5j) is then 

^^J^ = {^nx)A^iy))jHl - \{ [At{x)A^{y),HUl , (E.7) 

whereby [ , ]^ denotes the anticommutator (+), or the commutator (— ), respectively. The 
anticommutator in Eq. (IE.7j) can be rewritten as 

HU = -^f^a0^^ + ^nx)HA%y) + (E.8) 



using [Ai, [Aj,H]-]- = — ^ with Al{x) = i6/6Ill{x). With knowledge of the wave func- 
tional, 



AUxMA] = -j dS{u-')f^{x,y)j^^nA] , (E.9) 
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E. Identities for Gaussian expectation values 



the last two terms in Eq. (IE.8h result in 

{At{x)HA%y) + A'^iy)HA^{x))^ 



HA'>{y)-A'^iy)H 



I 



6H 



A'Ay)) +{HUtAn,y) 



d''v(^HA';iy)AUv)) u;'^Ju,v)- / d'^v {A>;iy)A^^{v)H) y^Ju,v) 
6H 



+ {H)^t%{u,y) 



Anx)A'j{y),H ) +2{u;-')t^{x,y){H) 



d^n(u;-i)-(x,n) 



SH 



(E.IO) 



The anticommutator in the last line can be rewritten as a commutator, according to the 
following identity, 



which holds for Gaussian expectation values, such as (|E.ip . 



fE.lll 



Proof: Inside the expectation value, the field operator A can be expressed by a momentum 
operator H, see Eq. IE.9[ Thus, 



{[At{x),0]^)^ 



dM^-X(-^yHjjF7-,o^o-. 



M (y) M (y) 



i I d''y{u;'X{x,y){[U,0]^)^ 



Therefore, 



6H 



6A<^{u) 



d%{u-Xn{y,V 



5A-^{u)5Ai{v) 



(E.12) 



(E.13) 



Plugging Eq. (IE.13P into the expression (|E.10p and the latter with (IE. 81) into Eq. (|E.7p . some 
cancellations occur and we end up with 



1 



--J d'^H(^-^)-(x,n) 
Such an identity was also derived in Ref. |154| . 



5^H 



6A-^{u)dAi{v) 



In] 



v,y). (E.14) 
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